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库在glm模型上运行bootcov。问题是，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)

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

1 个答案: