为逻辑回归添加预测（分组变量）

Question

我正在尝试通过分组进行逻辑回归 我尝试了已列出的所有方法，但仍然出现以下错误

UseMethod（“ predict”）中的错误：'predict'没有适用的方法 应用于类“ c（'tbl_df'，'tbl'，'data.frame'）”的对象

CharField

Answer 1

您可以将broom软件包中的augment()与nest / map / unnest结合使用。

library(purrr)
library(broom)

models <- df %>%
  group_by(organ) %>%
  nest() %>%
  mutate(model = map(data, ~ glm(IHC ~ Dose, data = ., family = binomial(logit)),
         augmented = map2(model, data, augment))

models %>%
  unnest(augmented)

预测值将在列.fitted中。

请注意，它们将成对数级：如果希望它们成为概率，则可能需要将type.predict = "response"添加为扩增函数的参数。