從 lme 模型中提取置信區間

Question

我必須對使用nlme包中的lme()函數制作的多個大型模型的置信區間進行一些轉換。 我正在使用intervals()函數來獲取間隔，但是無法將其轉換為數據幀。 有什么辦法可以把它變成可訪問的號碼嗎？

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

Answer 1

您可以使用CI$fixed獲取值，然后存儲在數據框中。 例如

df <- data.frame(CI$fixed)