如何在glm模型上使用bootcov（R）

Question

我正在使用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库运行bootcov的GLM模型。 问题是，bootcov似乎不想使用glm fit。 它只想使用lrms等。 我尝试对数据进行lrm模型建模，但是当我在那运行bootcov时，RStudio崩溃了。 多次。

有谁知道如何让bootcov与glms一起使用？ 我还有其他应该做的事情吗？ 我是R（和一般的统计学）的新手，实际上只是来自语言学背景。

谢谢！！

Answer 1

这是一段较长的评论时间，但并不是完整答案，我以前也没有使用过rms软件包。

通过查看bootcov的代码，它还会寻找类Glm 。 这是一个调用glm.fit的包装函数，因此应该可以。 您还需要在Glm调用中将x = TRUE和y = TRUE设置为bootcov精确说明这些。

小例子

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