[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():


#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

fit = lmer(
	formula = rt ~ (1|subnum) +
		cue*flank*group*block +
		cue*flank*group*I(block^2) +
	, data = ANT[ANT$error==0,]

More information about the R-sig-mixed-models mailing list