# How do we obtain p values from a robust mixed regression model in R?

I have yet to find an answer to this problem, so here goes.

I am fitting a robust mixed regression model using the robustlmm package on R. Unlike the lme4 package, from which we can obtain p values using the lmerTest package, the robustlmm package provides no method of obtaining p values from my model.

Does anyone have any experience with this issue? How can we calculate p values from a robust mixed model?

• Don't. P values are the devil's work and are not outputted by the software for good reasons (the same reasons lme4 does not output them). Resist any request from reviewers or editors for their inclusion and interpret the effect sizes instead. Aug 31 '20 at 13:04

If you really, really want to obtain p values, then one approach you can take is to fit the model with robustlmm::rlmer and an identical model with lmerTest::lmer and use the degrees of freedom reported by the latter to compute the p value for the t statistic reported by the former:
 pt(t.from.robustlmm, df.from.lmerTest, lower.tail = FALSE) * 2