354ページ6行目
lme1 <- lme(ex~cage*group,random=~agec|subject)
は
lme1 <- lme(log.ex~agec*group,rondom=~agec|subject)
では? このように修正すると,実行できて結果も本と一致します。 以下,summary(lme1)の出力中,exはlog.ex;cageはagecでしょう?
355ページの
lme2 <- update(lme1, corr=corCAR1())
は PowetPC G5, Mac OSX 10.4.8, R2.4.0の環境ではエラーになりました。
以下にエラーlme.formula(fixed = log.ex ~ agec * group, random = ~agec | subject, : nlminb problem, convergence error code = 1; message = iteration limit reached without convergence (9)