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]]

