Извлечение доверительных интервалов из модели IME - PullRequest
0 голосов
/ 21 ноября 2018

Я должен сделать некоторые преобразования на доверительных интервалах нескольких больших моделей, выполненных с помощью lme ().Я использую функцию interval (), чтобы получить интервалы, однако невозможно превратить его в информационный фрейм. Есть ли какой-нибудь способ, которым я мог бы превратить это в доступные номера?

library(nlme)
data(mtcars)
model = lme(mpg ~ disp + cyl, random=~1|gear,  data=mtcars)
CI<-intervals(model, which = "fixed")
Approximate 95% confidence intervals

 Fixed effects:
              lower        est.         upper
(Intercept) 29.43497446 34.66099475 39.8870150353
disp        -0.04163023 -0.02058363  0.0004629675
cyl         -3.04786061 -1.58727681 -0.1266930105
attr(,"label")
[1] "Fixed effects:"

data.frame(CI)
Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = 
stringsAsFactors) : 
  cannot coerce class ‘"intervals.lme"’ to a data.frame

1 Ответ

0 голосов
/ 21 ноября 2018

Вы можете использовать CI$fixed, чтобы получить значения, а затем сохранить их во фрейме данных.Например,

df <- data.frame(CI$fixed)
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...