[R] Collecting results of a test with array

Bert Gunter gunter.berton at gene.com
Fri May 25 16:10:59 CEST 2012


?abind  ## in R package abind. Very handy for this sort of thing

Something like:

do.call(abind, yourlist)

You can also do it by hand along the lines Michael described: e.g.
something like

array(unlist(yourlist), dim = c(2,3, 12) )

-- Bert

On Fri, May 25, 2012 at 6:16 AM, R. Michael Weylandt
<michael.weylandt at gmail.com> wrote:
> My first thought is to unlist them and stick them in a vector and then
> to form that vector into an array of the desired shape.
>
> Something like this:
>
> x <- vector("list",12)
> for(i in 6*(1:12)){x[[i/6]] <- matrix(i:(i+5), ncol = 2) }
>
> print(x)
>
> x.out <- array(unlist(x), dim = c(2,3,12))
>
> print(x.out)
>
> which I believe is what you want. [If not, I think it will at least
> get you started]
>
> Hope this helps,
>
> Michael
>
> On Fri, May 25, 2012 at 5:21 AM, Francesca <francesca.pancotto at gmail.com> wrote:
>> Dear contributors
>> I have tried this experiment:
>>
>> x<-c()
>> for (i in 1:12){
>> x[i]<-list(cbind(x1[i],x2[i]))  #this is a list of 12 couples of time
>> series I am using to perform a test
>> } # that compares them 2 by 2
>> #
>> #################
>> #trace statistic
>> test<-data.frame()
>> cval<-array( , dim=c(2,3,12))
>> for (i in 2:12){
>> for (k in 1:2){
>> for (j in 1:3){
>> result[k,j,i]<- ((ca.jo(data.frame(x[i]),ecdet="none",type="trace",
>>      spec="longrun",K=2))@cval[k,j])
>> }}}
>>
>> I have a problem in collecting the results of a
>> test.
>> The function ca.jo creates an object with various attributes,
>> one of which is the "cval" that i can access through @cval.
>> The attribute cval is an object of dimension 2X3.
>> I am running recursively the test with ca.jo for 12
>> couples of time series, so I have an output of 12 matrices of 2X3
>> elements and I would like to create an object like an array
>> of dimension (2,3,12) which contains each matrix @cval
>> produced by ca.jo for the 12 subjects that i tested.
>>
>> Can anyone help me with that?
>> I hope my explanation of the problem is clear.
>> Thanks in advance for any help.
>>
>> --
>>
>> Francesca
>>
>> ----------------------------------
>> Francesca Pancotto, PhD
>> Università di Modena e Reggio Emilia
>> Viale A. Allegri, 9
>> 40121 Reggio Emilia
>> Office: +39 0522 523264
>> Web: http://www2.dse.unibo.it/francesca.pancotto/
>> ----------------------------------
>>
>>        [[alternative HTML version deleted]]
>>
>>
>> ______________________________________________
>> R-help at 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.
>>
>
> ______________________________________________
> R-help at 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.



-- 

Bert Gunter
Genentech Nonclinical Biostatistics

Internal Contact Info:
Phone: 467-7374
Website:
http://pharmadevelopment.roche.com/index/pdb/pdb-functional-groups/pdb-biostatistics/pdb-ncb-home.htm



More information about the R-help mailing list