有没有办法计算R中GLM的估计值？

Question

我正在尝试获得R中的GLM模型的剩余系数估算。

这是我简单的循环版本实现：

data("mtcars")
res <- matrix(nrow = nrow(mtcars), ncol = 3L)
for (i in seq_len(nrow(mtcars))) {
  loo.data <- mtcars[-i, ]
  model <- glm(vs ~ mpg + disp, family = binomial, data = loo.data)
  res[i, ] <- model$coefficients
}

循环后，我得到了一个系数的矩阵，每行都包含一对一的估算。 res中的行数与观察总数相同。 res不必是矩阵格式。

是否有更清洁，也许更快地获得相同结果的方法？它将在其他C ++代码中使用。我在启动软件包中找到了cv.glm，但在我的情况下似乎不适用。使用RCPP 函数类可能是一个选项吗？在C ++中编写Fisher评分算法似乎只是浪费时间和精力。谢谢。