r - 绘制 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 给出了这个错误：

Error in `[<-`(`*tmp*`, use, use, ii, value = lst[[ii]]) : 
  subscript out of bounds

有没有人知道如何解决任何问题？

score 1 · Accepted Answer

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

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

1 回答 1