glht() дает тот же std. ошибка для каждого лечения

Я пытаюсь провести тест Даннета на линейной смешанной модели с использованием lme4 и glht. Я настроил и запустил модель, как показано ниже

Untransformed.lmer <- lmer(Sum ~ Treatment + (1|Block), data = EggCounts_poolSUM)
anova(Untransformed.lmer)
summary(glht(Untransformed.lmer, linfct = mcp(Treatment = 'Dunnett'), alternative = 'less'))

И когда я запускаю это, я получаю следующий результат

     Simultaneous Tests for General Linear Hypotheses

Multiple Comparisons of Means: Dunnett Contrasts


Fit: lmer(formula = Sum ~ Treatment + (1 | Block), data = EggCounts_poolSUM)

Linear Hypotheses:
             Estimate Std. Error z value Pr(<z)  
75 - 0 >= 0    -914.2      911.6  -1.003  0.372  
150 - 0 >= 0  -1207.4      911.6  -1.325  0.243  
300 - 0 >= 0  -2162.2      911.6  -2.372  0.030 *
600 - 0 >= 0  -1446.3      911.6  -1.587  0.160  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Adjusted p values reported -- single-step method)

Может кто-нибудь объяснить, как все методы лечения могут приводить к одному и тому же Std. ошибка? Что-то я делаю не так?

1 ответ

Стандартная ошибка Даннета для лечения i является sqrt(s2) * sqrt(1/ni + 1/n0) где s2 - объединенная оценка дисперсии, ni количество наблюдений для лечения i а также n0- количество наблюдений для контрольной группы. Так что стандартные ошибки одинаковы в случае, когдаniравны. Вероятно, это относится к вашим данным.

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