Проблема с функцией пар () в NLME - PullRequest
0 голосов
/ 03 января 2019

У меня проблемы с настройкой функции 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)

Но получаю ту же ошибку.Есть идеи?

1 Ответ

0 голосов
/ 03 января 2019

Вот как можно думать в этом случае. Ошибка

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. Загрузка или перезапуск сеанса почти наверняка помогут.

...