R 泊松GLM:相互作用的公式表示
我有一个关于泊松GLM和公式表示的问题: 考虑到数据集:R 泊松GLM:相互作用的公式表示,r,glm,poisson,R,Glm,Poisson,我有一个关于泊松GLM和公式表示的问题: 考虑到数据集: p <- read.csv("https://raw.githubusercontent.com/Leprechault/PEN-533/master/bradysia-greenhouse.csv") 最后一个公式是:bradysia=exp(??)和任何帮助,请?y~x1*x2在R中的意思与y~x1+x2+x1:x2相同,或者在数学上是y=a+b1*x1+b2*x2+b3*x1*x2。谢谢@DanY,但是关
p <- read.csv("https://raw.githubusercontent.com/Leprechault/PEN-533/master/bradysia-greenhouse.csv")
最后一个公式是:
bradysia=exp(??)
和任何帮助,请?y~x1*x2
在R中的意思与y~x1+x2+x1:x2
相同,或者在数学上是y=a+b1*x1+b2*x2+b3*x1*x2。谢谢@DanY,但是关于不显著系数的问题呢?因为在一个GLM中很难找到所有系数都具有显著性的分类和定量因子。
m1 <- glm(bradysia ~ area + mes, family="quasipoisson", data=p)
summary(m1)
#(Intercept) 4.36395 0.12925 33.765 < 2e-16 ***
#areaCV -0.19696 0.12425 -1.585 0.113
#areaMJC -0.71543 0.08553 -8.364 3.11e-16 ***
#mes -0.08872 0.01970 -4.503 7.82e-06 ***
m2 <- glm(bradysia ~ area*mes, family="quasipoisson", data=p)
summary(m2)
#(Intercept) 4.05682 0.15468 26.227 < 2e-16 ***
#areaCV 0.15671 0.35219 0.445 0.6565
#areaMJC 0.54132 0.31215 1.734 0.0833 .
#mes -0.03943 0.02346 -1.680 0.0933 .
#areaCV:mes -0.05724 0.05579 -1.026 0.3052
#areaMJC:mes -0.22609 0.05576 -4.055 5.57e-05 **