[R-sig-ME] Does lmer have it's own random number generator?
    Mike Lawrence 
    Mike.Lawrence at dal.ca
       
    Mon Jan  3 03:51:51 CET 2011
    
    
  
Hi folks,
Does lmer have it's own internally seeded random number generator? I
just came across an example where I get different results across
repeated identically specified lmer fits even when I use set.seed():
library(ez)
data(ANT)
#repeat the below a few times, the result
# on the last line should vary occasionally,
# sometimes TRUE, sometimes FALSE,
# and sometimes there will be a warning
# after the fit
set.seed(1)
fit = lmer(
	formula = rt ~ (1|subnum) +
		cue*flank*group*block +
		cue*flank*group*I(block^2) +
		cue*flank*group*I(block^3)
	, data = ANT[ANT$error==0,]
)
any(is.finite(attr(fit,'deviance')))
    
    
More information about the R-sig-mixed-models
mailing list