R を使用して、Pinheiro and Bates book (S および S-PLUS の混合効果モデル) の演習を解決しようとしています。nlme パッケージのデータ ピクセルを使用します。
演習では、lmList を使用して、犬ごとに 1 日の個別の 2 次モデルを当てはめます。lmList 近似の係数の個々の信頼区間をプロットします。切片と線形係数のみが Dog で大きく変化しているように見えることを確認します。
私は次のようにモデルを適合させました：
fm2Pixel.lis<-lmList(pixel~day+I(day^2)|Dog, Pixel)
しかし、この関数を使用して信頼区間をプロットしようとすると:
plot(intervals(fm2Pixel.lis))
R は mi にこのエラーを与えます:
Error in `[<-`(`*tmp*`, use, use, ii, value = lst[[ii]]) :
subscript out of bounds
誰かが間違っているものを修正する方法を知っていますか?