我无法让pairs()函数在nlme中工作。从Pinhiero和Bates的S和S-Plus中的混合效果模型中获取此示例。
模型本身运行良好
fm1Theo.lis <- nlsList(conc ~ SSfol(Dose, Time, lKe, lKa, lCl), data = Theoph)
但是两人情节...
pairs(fm1Theo.lis, id = 0.1)
...返回此错误
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)
但是得到同样的错误。有什么想法吗？
答案 0 :(得分：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软件包。几乎肯定可以加载它或重新启动会话。