使用lmer测试，以pmer表示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) 

谢谢 S