0

如何对随机效应 beta 分布进行事后检验？

为了构建具有 beta 分布的混合模型，我使用了库 gamlss（我发现在 R 中没有其他方法可以做到这一点）。在这个模型中，有一个具有三个水平的因子。假设 level1、level2 和 level3。摘要显示了 level1 x level2 和 level1 x level3 之间的比较，因为 level1 是基线。

library(gamlss) mymodel <- gamlss(response~ph+dose, family = BE, random=~1|id, data = mydata) 摘要(mymodel)

它显示了 level1 x level2 和 level1 x level3 之间的比较

对于所有三个级别，摘要显示 level1 x level2 和 level1 x level3 之间的比较，因为 level1 是基线。

我想测试level2 x level3之间是否有区别。我怎样才能做到这一点？也许可能有另一个图书馆。谁能帮帮我吗？

4

1 回答 1

-1

我建议你看看包“emmeans”。

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

祝你好运

蒂博

于 2019-11-12T11:05:31.980 回答