Проблема с выводом GLM - PullRequest
Новая
       72

Проблема с выводом GLM

0 голосов
/ 29 мая 2020

Итак, я хочу выполнить повторные измерения GLM в RStudio, что я сделал ... по большей части .... Однако не все даты отображаются в моем выводе (отсутствует 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 ?

1 Ответ

2 голосов
/ 29 мая 2020

Столбец Date является коэффициентом или преобразуется в коэффициент glm.

По умолчанию, glm использует контраст обработки для факторов, также называемый фиктивным кодированием. Это означает, что (Intercept) - это коэффициент для первого уровня 12/1/2015.

Для других дат эти коэффициенты отличаются от 12/1/2015. Например, перехват для 6/1/2017 на самом деле (Intercept) + Date6/1/2017 = 15.5643 + ( -18.3304).

Я мог бы дать более подробную c помощь, если бы вы предоставили воспроизводимый пример .

...