GLM output故障

Question

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

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

Answer 1

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

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

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

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