Displaying 3 results from an estimated 3 matches for "simu2".
Did you mean:
simul
2003 Mar 05
3
reading in tab delimited data in a loop
Dear all,
I need to read in 4 sets of tab delimited data in a loop. The 4 data sets are called "simu1.dat", "simu2.dat" and so on. I know what I need on the righthand side of the read.table expression but I can't the left hand side of it to work (see the line in bold below). Can you kindly help? Many thanks.
simu1 <- matrix(0,30,3)
simu2 <- matrix(0,30,3)
simu3 <- matrix(0,30,3)
simu4 <-...
2010 Apr 13
0
Error using reshape method
...imu11<-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...
2001 Apr 09
4
fastest R platform
...;
obj <- mle.b.D.u(b, D.u, y, x, z); print("mle"); print(obj$D.u);
write(obj$D.u, file = "simu1.dat", append=T);
obj <- REML.b.D.u(b, D.u, y, x, z); print("REML");
print(obj$D.u);
write(obj$D.u, file = "simu2.dat", append=T);
}
}
###################################################
#global variable m, n, n.fixed, n.random
n <- 10;
m <- 25;
n.fixed <- 4;
n.random <- 2;
main();
print(date());
=====
Jason G. Liao
Department of Biometry and Epidemiology
Me...