如何在 glm 模型上使用 bootcov (R)答案

【问题标题】：How to use bootcov on a glm model (R)如何在 glm 模型上使用 bootcov (R)
【发布时间】：2014-04-04 20:43:21
【问题描述】：

我正在使用 R 来分析语言数据。我有一个看起来像这样的数据框 -

phoneme  onset  voice  ident
b        TRUE   TRUE   TRUE
b        TRUE   TRUE   FALSE
b        TRUE   TRUE   FALSE
b        TRUE   TRUE   FALSE
b        TRUE   TRUE   FALSE
... etc

我在上面运行了glm，现在我需要使用rms 库在glm 模型上运行bootcov。问题是，bootcov 似乎不想使用 glm fit。它只想使用 lrms 之类的。我尝试对我的数据进行 lrm 模型，但是当我在上面运行 bootcov 时，RStudio 崩溃了。多次。

有谁知道如何让 bootcov 与 glms 一起工作？还有什么我应该做的吗？我是 R 的新手（以及一般的统计数据），实际上只是从语言学背景开始。

谢谢！！

【问题讨论】：

标签： r statistics linguistics

【解决方案1】：

评论有点长，但不是完整的答案，我之前也没有使用过rms 包。

通过查看bootcov 的代码，它还查找Glm 类。这是一个调用 glm.fit 的包装函数，所以应该没问题。您还需要在 Glm 调用中将 x=TRUE 和 y=TRUE 设置为 bootcovexpacts 这些。

小例子

# data
set.seed(1)
mydf <- data.frame(replicate(3,rbinom(100,1,0.3)))
names(mydf) <- letters[1:3]

# ------------------------------------------------------------
library(rms)

(mod1 <- Glm(a ~ b + c  , data = mydf , family="binomial", x=TRUE , y=TRUE))

# Check with glm
(summary(mod2 <- glm(a ~ b + c  , data = mydf , family="binomial", x=TRUE , y=TRUE)))

# Run bootstrap
set.seed(1)
b <- bootcov(mod1, B=200)

【讨论】：