我想使用plot()
和curve()
函数生成预测计数图。我运行了一个数据集相对较小的计数模型。
我可以使用哪些代码来生成预测计数?
我的因变量是政府部门收到的请愿数量。
我的自变量是一个从0到3的连续变量。
我的回归结果为all.glm
我已使用以下代码尝试此操作:
# Extract Coefficients and Alpha
all.alpha <- coef(all.glm)[1]
all.coeff <- coef(all.glm)[2]
# Generate synthetic values
varx <- seq(0, 3, 0.25)
# Generate Plot
plot(-10, -10,)
# Add curves
for (i in 1:length(corrx)) {
curve(pnorm(all.alpha*1 + all.coeff*varx[i,]), add = TRUE, lwd = 2)
}
正如您所看到的,我创建了一个假的自变量值向量,我试图迭代它们以生成各种预测结果。
不幸的是,我收到了这个错误:
Error in curve((all.alpha * 1 + all.coef * 1.5), add = TUE, lwd = 2) :
'expr' must be a function, or a call or an expression containing 'x'
我确定这是一个愚蠢的错误,但任何帮助都将不胜感激!