如何在 beta 发行版中对混合模型执行事后测试(我正在使用 R 的 gamlss 库)？

Question

我如何执行随机效果beta分布的事后测试？

为了构建与beta分发的混合模型，我使用了库gamlss(我没有其他方法可以在R中执行此操作(。在此模型中，有一个具有三个级别的因素。假设Level1，Level2和Level3。摘要显示了Level1 X Level2和Level1 X Level3之间的比较，因为Level1是基线。

库(gamlss(mymodel＆lt; - gamlss(响应〜ph 剂量，family = be，andan = 〜1 | id，data = mydata(摘要(mymodel(

它显示了Level1 X Level2和Level1 X Level3

之间的比较

对于所有三个级别，摘要都显示了Level1 X Level2和Level1 X Level3之间的比较，因为Level1是基线。

我想测试Level2 x级别3之间是否存在差异。我怎样才能做到这一点？也许可能还有另一个库。任何人都可以帮我吗？

Answer 1

我建议您看看" emmeans"软件包。

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

祝你好运

tibo