我的R代码,我有以下内容:
input.formula = 'some_valid_formula'
glm.model <- glm(inp.formula,data=data, family='quasipoisson')
summary <- capture.output(summary(glm.model))
print(summary)
问题在于摘要打印:
调用:[INFO] - 模型系数:glm(公式= input.formula,family = \&#34; quasipoisson \&#34) ...
请注意,公式将打印为公式变量名称,而不是公式本身。我在这里缺少什么?
答案 0 :(得分:1)
一种选择是使用do.call
来构建和执行glm
调用。这将返回已显示公式的已评估glm
调用。例如:
data <- data.frame(group = gl(3,3), type = gl(3,1,9), counts = rpois(9, 2))
input.formula = "counts ~ group + type"
glm.model <- do.call("glm", list(input.formula, 'quasipoisson', data))
summary(glm.model)