同じく，farawayさんの文章から。mcmcsampという関数があるのですね。参考。しかし，farawayさんの文章を読むかぎり，lmerはなかなか扱いが難しそう...

farawayさんの文章を読んでいて，simulateという関数があるということを知りました。これは便利。library(lme4)内にはlmer用のsimulateもあるよう。 > xx > yy > plot(yy ~ xx, , cex = 2) > m1 > summary(m1) Call: lm(formula = yy ~ xx)Residuals: Min 1Q…