To be fair I don't really understand what the error is telling me so I
don't
know how to correct it.
I'm trying to reshape my data so that I do repeated measures ANOVA on it.
The warning message: "Error in d[, timevar] <- times[1L] : subscript out
of
bounds"
The program:
library(QRMlib)
library(Hmisc)
ss<-5
mu<- c(0,2,3,12)
mu2<- c(0,2,10,20)
dims<- 4
co<- c(.3,.4,.5,.4,.3,.2)
replis<- 3
stdev<-c(1,2,5,10)
#1st set
corrmat <- matrix(1, dims, dims)
corrmat[upper.tri(corrmat)]<-co
selectlow <- lower.tri(corrmat)
corrmat[selectlow] <- t(corrmat)[selectlow]
sigma<- diag(stdev) %*% corrmat %*% diag(stdev)
simu1<-rmnorm(n=ss, Sigma=sigma, mu=mu, d=dims, rho=co)
simu11<-cbind(simu1,data.grp=1)
#2nd set
corrmat2 <- matrix(1, dims, dims)
corrmat2[upper.tri(corrmat2)]<-co
selectlow2 <- lower.tri(corrmat2)
corrmat2[selectlow2] <- t(corrmat2)[selectlow2]
sigma2<- diag(stdev) %*% corrmat2 %*% diag(stdev)
simu2<-rmnorm(n=ss, Sigma=sigma2, mu=mu2, d=dims, rho=co)
simu22<-cbind( simu2,data.grp=2)
data.all<-rbind(simu11,simu22)
new.data.all<-rbind(data.all)
data.all<-rbind(data.all)
data.frame(data.all)
Y<-cbind(data.all[,1], data.all[,2], data.all[,3], data.all[,4])
rmlong.data<-reshape(data.all, varying=list(names(data.all)),
direction="long" ) # Part that's giving me trouble.
Thank you for your help!
--
View this message in context:
http://n4.nabble.com/Error-using-reshape-method-tp1837911p1837911.html
Sent from the R help mailing list archive at Nabble.com.
