Thank you very much for your suggestions. They both worked!
Actually I am trying to make an array of length 7, where each array
contains a different number of 2 by 2 matrices. The vector "n" defines the
different number of matrices across the arrays. Therefore, the parameter
"l" has to change in every array i. Of course I cannot put the vector "n"
directly into the dimensions of the array (under the name "results").
First, I wanted to see if the code gives rational results assuming that
each array has the same number of 2 by 2 matrices. Fortunately, the results
are rational (thanks to your help, as well). Now I am taking the problem to
the next level and I want to construct this array of different number of 2
by 2 matrices!
Any suggestions would be really very helpful!
On Tue, Oct 2, 2012 at 6:23 PM, Rui Barradas <ruipbarra...@sapo.pt> wrote:
> Hello,
> See if this is it.
>
>
> Nx <- rep(0,length(x))
> Ny <- rep(0,length(y))
> n <- (x+1)*(y+1)
> results <- array(0, dim=c(2,2,64,7))
> # l <- 1 # <--------------------------- This changed place
>
> for(i in 1:length(x)){
> Nx[i] <- length(1:(x[i]+1))
> Ny[i] <- length(1:(y[i]+1))
> l <- 1 # <----------------------- To here
>
> for(j in 1:(Nx[i])){
> for(k in 1:(Ny[i])){
>
>
> Hope this helps,
>
> Rui Barradas
> Em 02-10-2012 15:20, Loukia Spineli escreveu:
>
>> I want to make a multi-dimensional array. To be specific I want to make
>> the
>> following array
>>
>> results<-array(0,dim=c(2,2,64,**7))
>>
>> This is the code I have created but it gives no result due to the error
>> "subscript out of bound".
>>
>> x<-rep(7,7) # Missingness in intervention
>> y<-rep(7,7) # Missingness in control
>>
>> arraynames<-list(Group=c("**Success","Failure"),Outcome=c(**
>> "Intervention","Control"),**Trial=c("1":"7"))
>> mat.stat<-array(c(9,16,10,15,**66,12,44,23,102,88,66,104,277,**
>> 60,247,119,23,43,20,41,201,**162,122,
>> 263,14,41,4,41),dim=c(2,2,7),**dimnames=arraynames);mat.stat
>>
>> Nx<-rep(0,length(x))
>> Ny<-rep(0,length(y))
>> n<-(x+1)*(y+1)
>> results<-array(0,dim=c(2,2,64,**7))
>> l<-1
>> for(i in 1:length(x)){
>> Nx[i] <- length(1:(x[i]+1))
>> Ny[i] <- length(1:(y[i]+1))
>> for(j in 1:(Nx[i])){
>> for(k in 1:(Ny[i])){
>> results[,,l,i]<-mat.stat[,,i]+**matrix(c(c(0:x[i])[j],c(0:y[i]**
>> )[k],-c(0:x[i])[j],
>> -c(0:y[i])[k]),nrow=2,ncol=2,**byrow=T)
>> l<-l+1
>> }
>> }
>> }
>> results
>>
>> Any suggestion would be really welcome!
>> Thank you very much!
>>
>> Loukia
>>
>> [[alternative HTML version deleted]]
>>
>> ______________________________**________________
>> R-help@r-project.org mailing list
>> https://stat.ethz.ch/mailman/**listinfo/r-help<https://stat.ethz.ch/mailman/listinfo/r-help>
>> PLEASE do read the posting guide http://www.R-project.org/**
>> posting-guide.html <http://www.R-project.org/posting-guide.html>
>> and provide commented, minimal, self-contained, reproducible code.
>>
>
>
[[alternative HTML version deleted]]
______________________________________________
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.
