R lsmeans指调整多重比较

Question

我使用lme4在R中运行混合效果的逻辑回归(通过调用glmer)，现在我正在尝试进行临时比较。因为他们是对的，图基应该是好的，但我想手动调整多少测试应该作出修正-现在它是为12个测试，但我只是在6个比较的兴趣。

到目前为止，我的代码看起来像这样

    for (i in seq_along(logmixed_ranks)) {
    print(lsmeans(logmixed_ranks[[i]], pairwise~rating_ranks*indicator_var, adjust="tukey"))
    }

不知何故，我可能需要使用以下内容，但我不确定如何使用。

      p.adjust(p, method = p.adjust.methods, n = length(p))

有人能帮忙吗？谢谢!劳拉

Answer 1

一定有一个原因，你只想调整6比较，我猜是因为你想打破比较，你所做的有条件的因素之一。这很容易使用lsmeans

lsmeans(logmixed_ranks[[i]], 
    pairwise ~ rating_ranks | indicator_var, adjust = "tukey")

或

lsmeans(logmixed_ranks[[i]], 
    pairwise ~ indicator_var | rating_ranks, adjust = "tukey")

顺便说一句，如果您使用adjust = "mvt"，您将获得与glht用于其单步过程所使用的完全相同的调整。因此，我认为唯一不受lsmeans支持的lsmeans特性是多步测试。

我对您为什么有一个glmer对象列表感到困惑，但这似乎与我的回答无关。