Hello there
someone who knows this?
im having some trouble with making a function that comes out with what
i want.
i have these vectors:
x1 <- factor(rep(1:2500,3*10))
x2 <- factor(rep(1:25000,3))
x3 <- factor(1:75000)
y <- rep(rnorm(2500,mean=0,sd=2),10*3)+
rep(rnorm(25000,mean=0,sd=sqrt(2)),3)+rnorm(75000,mean=0,sd=1)
its a simulation of data and i want to run them fx 1000 times, and in
every run make an output from this
lmer1 <- lmer(y~1+(1|x1)+(1|x2))
with the variance on x1 and x2 and residual (fx in 3 vectors or a
dataframe or something....
is this possible? is there someone who knows this?
thanks in advance
Rina
[[alternative HTML version deleted]]
______________________________________________
R-help@stat.math.ethz.ch 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.
