r语言 - GLM 输出故障

Question

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

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 输出中显示的日期之一？

Answer 1

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

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

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

如果您提供一个可重现的示例，我可以提供更具体的帮助。