Hi again...
I will get an array with 5 matrix 2x2, with dimension 2x2x5
This code below works fine... And I am applying the same procedure...
The problem might be when I do the permutation of the array, but I checked
the dimension and its all right...
array1 <- array(1:30,dim=c(3,2,5))
array2 <- array(1:20,dim=c(2,2,5))
result <- array(dim=c(dim(array1)[1], dim(array2)[2], dim(array1)[3]))
for (i in 1: dim(array1)[3]){
result[,,i] <- array1[,,i]%*%array2[,,i]
}
Thanks...
Ben Bolker wrote:
>
>
>
> MarcioRibeiro wrote:
>>
>> Hi listers,
>> I am having some trouble in a matrix multiplication...
>> I have already checked some posts, but I didn't find my problem...
>> I have the following code...
>> But I am not getting the right multiplication...
>> I checked the dimension and they are fine...
>>
>> id_y <- array(1:10,dim=c(2,1,5))
>> id_yt<-aperm(id_y,c(2,1,3))
>> m_id<-array(dim=c(dim(id_y)[1],dim(id_y)[1],dim(id_y)[3]))
>> for (i in 1:dim(id_y)[3]){
>> m_id[,,i]<-id_y[,,i]%*%id_yt[,,i]
>> }
>> m_id
>>
>>
>
> What did you expect to happen?
>
> Ben Bolker
>
>
--
View this message in context:
http://www.nabble.com/Matrix-multiplication---code-problem-tp22835003p22836042.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
