GLM输出故障

问题描述 投票:0回答:1

所以,我想在RStudio中运行重复测量GLM,我已经做了......大部分......。然而,并不是所有的日期都显示在我的输出中(1212015不见了)。这里是部分输出和我的模型代码,你可以看到我的意思。

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  

我如何让1212015成为我的GLM输出中的一个日期?

r rstudio glm
1个回答
2
投票

栏目 Date 是一个因数,或正被以下因素转换为一个因数 glm.

默认情况下。glm 使用治疗对比的因素,也称为Dummy Encoding。这意味着 (Intercept) 是第一级的系数。12/1/2015.

就其他日期而言,这些系数是指与下列日期相比的变化。12/1/2015. 例如,对 6/1/2017 其实 (Intercept) + Date6/1/2017 = 15.5643 + ( -18.3304).

如果你能提供一个更具体的帮助。可复制的例子.

© www.soinside.com 2019 - 2024. All rights reserved.