我有一些计数数据分为两类,并且已经进行了泊松GLM以寻找两个类别的差异:
mod <- glm(count ~ category, family = "poisson")
对于海报,我想要一个简单的结果图 - 可能是带有误差条的条形图。有没有一种简单的方法可以从mod
获取情节数据?误差线是否对应于模型Std. Error
?对于它的价值,我打算用ggplot2
来制作数字。
另外,我认为,对于Poisson GLM,模型单位是“实际”单位的对数 - 所以如果Estimate
的{{1}}是0.4,那么这个对应于B类中的观察结果,平均而言,exp(0.4)= 1.5比A类更高?