我的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 = \“quasipoisson \”)...
请注意,公式将打印为公式变量名称,而不是公式本身。我在这里错过了什么?
一种选择是使用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)