为 lmList 拟合中的系数绘制单独的置信区间

Question

我正在尝试使用 R 解决 Pinheiro 和 Bates 书中的练习（S 和 S-PLUS 中的混合效果模型）。它使用 nlme 包中的数据像素。

练习说：使用 lmList 为每只狗在一天中拟合一个单独的二次模型。为 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

有谁知道如何解决问题吗？

Answer 1

Dog == 9 只有两个观测值。这导致二次参数的估计值 NA 而 intervals 无法处理。如果排除此子集，它会起作用：

fm2Pixel.lis <- lmList(pixel ~ poly(day, 2, raw = TRUE) | Dog,
                       Pixel[Pixel$Dog != 9,])
plot(intervals(fm2Pixel.lis))

Plotting individual confidence intervals for the coefficients in the lmList fit