我在R中估计了以下泊松回归:
poissmod2 <- glm(aerobics$y ~ factor(aerobics$x1) + aerobics$x3 + aerobics$x4,
family = poisson))
到目前为止没问题。但是,下一步是确定aerobics$x1
的95%置信区间。
我尝试过以下两个代码:
confint.default(poissmod2,"aerobics$x1")
和
confint.default(poissmod2,"factor(aerobics$x1)")
对于这两个代码,我得到一个置信区间'Na-Na'...... 但是,这段代码似乎适用于非因子变量。
有没有人知道我做错了什么?