如何对 Beta 分布中的混合模型进行事后测试（我正在使用 R 的 gamlss 库）？

Question

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

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

图书馆（gamlss） 我的模型<- gamlss(response~ph+dose, family = BE, random=~1|id, data = mydata) summary(mymodel)

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

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

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

Answer 1

emmeans 尚不适用于 gamlss 对象

您可以尝试 mgcv::anova.gam(model) 它将为您打印摘要，但输出列表很难使用，并且无法轻松转换为数据框。

Answer 2

我建议你看看“emmeans”包。

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

祝你好运

蒂博