GLM输出故障答案

【问题标题】：GLM output troubleGLM输出故障
【发布时间】：2020-05-30 03:19:45
【问题描述】：

所以，我想在 RStudio 中运行重复测量 GLM，我已经完成了......在大多数情况下......但是，并非所有日期都显示在我的输出中（缺少 2015 年 12 月 1 日）。这是输出的一部分以及我的模型代码，因此您可以了解我的意思：

CH4f1 <- glm(GC_CH4.flux~River*Site*Date*Hum.Hol, data = Rdata_w.o_OL_Date, family = gaussian)
 summary(CH4f1)

Call:
glm(formula = GC_CH4.flux ~ River * Site * Date * Hum.Hol, family = gaussian, 
    data = Rdata_w.o_OL_Date)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-37.307   -2.655   -0.341    0.314  163.247  

Coefficients:
                                 Estimate Std. Error t value Pr(>|t|)  
(Intercept)                       15.5643    75.6869   0.206   0.8372  
River                             -8.4098    55.1093  -0.153   0.8788  
Site                              -3.4776    25.9063  -0.134   0.8933  
Date12/1/2016                    112.5939    96.6623   1.165   0.2451  
Date4/1/2016                     -15.4780    96.9954  -0.160   0.8733  
Date4/1/2017                     -13.8752    94.5132  -0.147   0.8834  
Date6/1/2016                      12.5824    93.8721   0.134   0.8935  
Date6/1/2017                     -18.3304    94.5132  -0.194   0.8464  
Date9/1/2016                     170.2484    95.5697   1.781   0.0759 .
Date9/1/2017                     -38.4031    96.7184  -0.397   0.6916

如何让 12/1/2015 成为我的 GLM 输出中显示的日期之一？

【问题讨论】：

标签： r glm

【解决方案1】：

Date 列是一个因子，或者正在被glm 转换为一个因子。

默认情况下，glm 使用因子的处理对比，也称为虚拟编码。这意味着(Intercept) 是第一级的系数12/1/2015。

对于其他日期，这些系数是 12/1/2015 的变化。例如，6/1/2017 的截距实际上是(Intercept) + Date6/1/2017 = 15.5643 + ( -18.3304)。

如果您提供reproducible example，我可以提供更具体的帮助。

【讨论】：