当我在glm模型中使用交互加上family = binomial参数时以及当我将其遗漏时,我有一个关于重要性和重要性差异的问题 . 我对逻辑回归很陌生,过去只进行了更简单的线性回归 .

我有一个 observations of tree growth rings 的数据集,其中包含 two categorical explanatory variablesTreatmentOrigin ) . 治疗变量是一个实验性干旱治疗,有四个级别(控制,第一次干旱,第二次干旱和两次干旱) . Origin变量有三个级别,指的是树的原点(给定代码颜色表示不同的来源为红色,黄色和蓝色) . 我的观察结果是否存在增长环( 1 = growth ring present, 0 = no growth ring ) .

就我而言,我对治疗效果,Origin的效果以及治疗和起源的可能相互作用感兴趣 .

有人建议二项Logistic回归是分析该数据集的好方法 . (希望这是合适的?也许有更好的方法?)

我有n = 5(每个原始治疗组合有5个观察结果 . 因此,例如,对照处理蓝色原生树的5个生长环观察结果,对照处理黄色原生树的5个观察结果等)所以总共数据集中有60个生长环观测值 .

在R中,我使用的代码是glm()函数 . 我把它设置如下:growthring_model < - glm(growthringobs~Relation Origin Treatment:Origin,data = growthringdata,family = binomial(link =“logit”))

我已经考虑了我的解释变量,因此控制处理和蓝色原点树是我的参考 .

我注意到的是,当我从代码中留下“family = binomial”参数时,它给出了我在给定数据结果时合理预期的p值 . 但是,当我添加“family = binomial”参数时,p值为1或非常接近1(例如,1,0.98,0.99) . 这看起来很奇怪 . 我可以看到它的重要性很低,但是根据我的实际数据,这些值接近于1会让我产生怀疑 . 如果我在不使用“family = binomial”参数的情况下运行模型,我会得到似乎更有意义的p值(即使它们仍然相对较高/无意义) .

Can someone help me to understand how the binomial argument is shifting my results so much? (我知道它指的是分布,即我的观察结果是1或0)它在模型中究竟发生了什么变化?这是样本量小的结果吗?我的代码中有什么东西吗?也许那些非常高的值是正确的(或不是?)?

以下是我的模型摘要的读出,其中包含二项式参数:调用:glm(公式=生长期〜治疗原点治疗:起源,家庭=二项式(link =“logit”),data = growthringdata)

Deviance Residuals: 
     Min        1Q    Median        3Q       Max  
-1.79412  -0.00005  -0.00005  -0.00005   1.79412  

Coefficients:
                                       Estimate Std. Error z value Pr(>|z|)
(Intercept)                          -2.057e+01  7.929e+03  -0.003    0.998
TreatmentFirst Drought               -9.931e-11  1.121e+04   0.000    1.000
TreatmentSecond Drought               1.918e+01  7.929e+03   0.002    0.998
TreatmentTwo Droughts                -1.085e-10  1.121e+04   0.000    1.000
OriginYellow                          1.918e+01  7.929e+03   0.002    0.998
OriginRed                            -1.045e-10  1.121e+04   0.000    1.000
TreatmentFirst Drought:OriginYellow  -1.918e+01  1.373e+04  -0.001    0.999
TreatmentSecond Drought:OriginYellow -1.739e+01  7.929e+03  -0.002    0.998
TreatmentTwo Droughts:OriginYellow   -1.918e+01  1.373e+04  -0.001    0.999
TreatmentFirst Drought:OriginRed      1.038e-10  1.586e+04   0.000    1.000
TreatmentSecond Drought:OriginRed     2.773e+00  1.121e+04   0.000    1.000
TreatmentTwo Droughts:OriginRed       2.016e+01  1.373e+04   0.001    0.999

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 57.169  on 59  degrees of freedom
Residual deviance: 28.472  on 48  degrees of freedom
AIC: 52.472

Number of Fisher Scoring iterations: 19

这里是没有二项式参数的模型摘要的读出:调用:glm(公式=增长〜治疗原点治疗:原点,数据= growthringdata)

Deviance Residuals: 
Min      1Q  Median      3Q     Max  
-0.8     0.0     0.0     0.0     0.8  

Coefficients:
                                   Estimate Std. Error t value Pr(>|t|)  
(Intercept)                          -4.278e-17  1.414e-01   0.000           1.0000  
TreatmentFirst Drought                3.145e-16  2.000e-01   0.000   1.0000  
TreatmentSecond Drought               2.000e-01  2.000e-01   1.000   0.3223  
TreatmentTwo Droughts                 1.152e-16  2.000e-01   0.000   1.0000  
OriginYellow                          2.000e-01  2.000e-01   1.000   0.3223  
OriginRed                             6.879e-17  2.000e-01   0.000   1.0000  
TreatmentFirst Drought:OriginYellow  -2.000e-01  2.828e-01  -0.707   0.4829  
TreatmentSecond Drought:OriginYellow  2.000e-01  2.828e-01   0.707   0.4829  
TreatmentTwo Droughts:OriginYellow   -2.000e-01  2.828e-01  -0.707   0.4829  
TreatmentFirst Drought:OriginRed     -3.243e-16  2.828e-01   0.000   1.0000  
TreatmentSecond Drought:OriginRed     6.000e-01  2.828e-01   2.121   0.0391 *
TreatmentTwo Droughts:OriginRed       4.000e-01  2.828e-01   1.414   0.1638  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for gaussian family taken to be 0.1)

    Null deviance: 8.9833  on 59  degrees of freedom
Residual deviance: 4.8000  on 48  degrees of freedom
AIC: 44.729

Number of Fisher Scoring iterations: 2

(我提前为我的问题可能的简单性道歉 . 我已经尝试阅读逻辑回归并试图遵循一些例子 . 但我一直在努力寻找解决我的特定情况的答案)

非常感谢 .