Как интерпретировать коэффициент лог-преобразованной модели из смешанной модели в исходном масштабе в R?

Полная модель для роста растений выглядит следующим образом:

lmer(log(growth) ~ nutrition + fertilizer + season + (1|block)

где питание (азот / фосфор), удобрения (нет / добавлено), сезон (сухой / влажный)

Краткое изложение модели выглядит следующим образом:

REML criterion at convergence: 71.9

Scaled residuals: 
     Min       1Q   Median       3Q      Max 
-1.82579 -0.59620  0.04897  0.62629  1.54639 

Random effects:
 Groups   Name        Variance Std.Dev.
 block     (Intercept) 0.06008  0.2451  
 Residual             0.48633  0.6974  
Number of obs: 32, groups:  tank, 16

Fixed effects:
                        Estimate Std. Error      df t value Pr(>|t|)    
(Intercept)               3.5522     0.2684 19.6610  13.233 3.02e-11 ***
nutritionP                0.2871     0.2753 13.0000   1.043  0.31601    
fertlizeradded           -0.3513     0.2753 13.0000  -1.276  0.22436    
seasonwet                 1.0026     0.2466 15.0000   4.066  0.00101 ** 
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Рост растений здесь зависит только от сезона, и увеличение роста составляет 1,0026 в логарифмическом масштабе. Как мне интерпретировать это в масштабе исходных данных, если я хочу узнать, каким было увеличение фактической высоты растения? Это всего лишь e(1.0026) ~ 3 см или есть какой-то другой способ интерпретировать это?

1 ответ

Решение

exp(1.0026) действительно около 3 (2,72), но это значение представляет собой пропорциональное изменение. Во влажный период рост в три раза выше, чем в сухой, при прочих равных условиях.

Другие вопросы по тегам