Парные сравнения с emmeans

Question

У меня есть вопрос по поводу вывода emmeans.Я сделал MCMCglmm, а затем emmeans для парного сравнения.Выходные данные содержат «оценка», «нижний.HPD» и «верхний.HPD».Какова оценка здесь?Это журнал?

Спасибо!Джен

Answer 1

MCMCglmm использует байесовские методы, и, соответственно, emmeans по умолчанию выдает байесовскую сводку, показывающую интервалы задней средней и максимальной плотности для каждой EMM.Если вы хотите (возможно, более знакомое) резюме для часто используемых пользователей, используйте что-то вроде

emm = emmeans(...)
summary(emm, freq = TRUE)
summary(pairs(emm), freq = TRUE)

Если вы оснастили модель функцией связи, тогда сводки будут в масштабе функции связи (например, log) если только вы не добавите type = “response” к этим summary вызовам.

Надеюсь, это поможет.