从 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)

从 lme 模型中提取置信区间

Extracting confidence intervals from lme model

r

confidence-interval

nlme