Hi,
Although my question is not directly linked to R functionality, I hope
you can forgive me for posing it here. I have been looking for the
answer for a long time (~ 4 weeks) and have not been able to find it. My
question is:
Suppose I have an m*n matrix, with a random (normally distributed)
What do you think about the following:
library(nlme)
set.seed(1)
n - 3;m - 4
s.e - 0
(X0 - array(rep(rnorm(n), each=m)+s.e*rnorm(m*n),
dim=c(m, n)))
s.e - 1
(X1 - array(rep(rnorm(n), each=m)+s.e*rnorm(m*n),
dim=c(m, n)))
X. - data.frame(Row=as.vector(row(X)),
