Как да извърша post hoc тест за смесен модел в бета разпространение (използвам библиотеката gamlss на R)?

Как да извърша пост-хок тест за бета разпределение със случаен ефект?

За изграждане на смесения модел с бета разпространение използвах библиотеката gamlss (не намерих друг начин да направя това в R). В този модел има фактор с три нива. Да предположим ниво 1, ниво 2 и ниво 3. Обобщението показва сравнението между level1 x level2 и level1 x level3, тъй като level1 е базовата линия.

библиотека(gamlss) mymodel ‹- gamlss(отговор~ph+dose, семейство = BE, случайно=~1|id, данни = mydata) резюме (mymodel)

Показва сравнението между level1 x level2 и level1 x level3

И за трите нива резюмето показва сравнението между ниво 1 x ниво 2 и ниво 1 x ниво 3, тъй като ниво 1 е базовата линия.

Бих искал да тествам дали има разлика между level2 x level3. Как мога да направя това? Може би може да има друга библиотека за това. Може ли някой да ми помогне?


r mixed-models posthoc beta-distribution gamlss
person Omar    schedule 30.07.2019    source източник

Отговори (1)


arrow_upward
-1
arrow_downward

Предлагам ви да разгледате пакета "emmeans".

    library("emmeans")
    Post_Hoc <- emmeans (model, specs = pairwise ~ Var, type = "response")
    summary(Post_Hoc)

Късмет

Тибо

person Tibo    schedule 12.11.2019