r - 使用 lmer 测试的 lmer 中的 p 值。多次运行模型以获得稳健的 p 值？

Question

我使用 lmerTest 在 R 中为我的 lmer 模型获取 p 值。不幸的是，多次运行模型每次都会给我不同的 p 值。有没有办法多次运行模型来计算一个稳健的 p 值？任何其他获得更稳健 p 值的方法也非常受欢迎。

我的代码（这是根据 Aicc 具有最佳随机结构的模型，已删除与高变化膨胀因子的交互）：

model <- lmer(cmic~h2of+pH+TempGlob+exp.age+tree.density+log.div+
              h2of:pH+
              h2of:TempGlob+
              h2of:log.div+
              pH:log.div+
              TempGlob:log.div+
              exp.age:log.div+
              (1|experiment/site/block),data=micc2,REML=TRUE) 

谢谢小号