nlme中双()函数的麻烦

Question

我很难让pairs()函数在nlme中工作。以S和S中的Pinhiero和Bates混合效应模型为例.

模型本身运行得很好

fm1Theo.lis <- nlsList(conc ~ SSfol(Dose, Time, lKe, lKa, lCl), data = Theoph)

但这对夫妇..。

pairs(fm1Theo.lis, id = 0.1)

...returns这个错误

Error in as.data.frame.default(x) : 
  cannot coerce class "c("nlsList", "lmList")" to a data.frame

我也试过

pairs(fm1Theo.lis, ~ ranef(., level = 2), id = 0.1)

但也有同样的错误。有什么想法吗？

Answer 1

在这种情况下，人们可能会这样想。误差

Error in as.data.frame.default(x) : 
  cannot coerce class ‘c("nlsList", "lmList")’ to a data.frame

表示c("nlsList", "lmList")类的某个对象正被胁迫到数据帧中。现在，由于fm1Theo.lis是使用nlsList的结果，因此错误中的对象似乎确实是nlsList。这意味着pairs不知道如何处理此类对象。要确认这一点，我们可以运行

pairs.default(fm1Theo.lis, id = 0.1)

当没有找到fm1Theo.lis的特定方法时，就会发生这种情况。事实上，错误是一样的。以某种方式确认nlsList和来自nlme的消息，很明显问题在于加载nlme包。加载它或重新启动会话几乎肯定会有所帮助。