所以,我想在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输出中的一个日期?
栏目 Date
是一个因数,或正被以下因素转换为一个因数 glm
.
默认情况下。glm
使用治疗对比的因素,也称为Dummy Encoding。这意味着 (Intercept)
是第一级的系数。12/1/2015
.
就其他日期而言,这些系数是指与下列日期相比的变化。12/1/2015
. 例如,对 6/1/2017
其实 (Intercept) + Date6/1/2017 = 15.5643 + ( -18.3304)
.
如果你能提供一个更具体的帮助。可复制的例子.