r - lmer Test を使用した lmer の p 値。堅牢な p 値を取得するためにモデルを複数回実行しますか?

Question

私は lmerTest を使用して、R で lmer モデルの p 値を取得しています。残念ながら、モデルを複数回実行すると、毎回異なる p 値が得られます。モデルを複数回実行して、1 つの堅牢な 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