[R] Collecting results of a test with array

[R. Michael Weylandt]

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.
>