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

我用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。它只想使用LRM等等。我试着对我的数据做一个lrm模型，但当我在上面运行bootcov时，RStudio崩溃了。多次

有人知道如何让bootcov与glms合作吗？还有什么我应该做的吗？我对R（以及一般的统计学）是个新手，我只是从语言学的角度来研究这个问题的

---

谢谢

这篇评论有点长，但不是完整的回答，我以前也没有使用过

rms

软件包

通过查看

bootcov

的代码，它还可以查找类

Glm

。这是一个调用glm.fit的包装函数，所以应该是ok。您还需要在Glm调用中将x=TRUE和y=TRUE设置为

bootcov

expacts

小例子

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

#数据
种子（1）
mydf