Dear All,
Suppose I have a regular linear mixed model
Y|b~ X*beta+Z*b+error;
b~N(0, Psi)
and would like to fit it to a simulated data, where Y, X and Z are known.
I checked some reference but examples are given with the specified variable
names. Could some one help me out, please?
# R code to try
require(MASS)
require(lme4)
n <- 20
beta <- c( -0.5, 1)
psi <- matrix(c(1,0.2,0.2,1),2,2)
b <- mvrnorm(n=1, mu=c(0,0), Sigma=psi )
X <- cbind(rep(1,n),1:n)
Z <- X
Y <- X%*%beta+ Z%*%b+rnorm(n,0,1)
?lmer
Thank you.
Best wishes,
Jie
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
