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
равны. Вероятно, это относится к вашим данным.