Значения p в lmer с использованием теста lmer. Запускать модель несколько раз, чтобы получить надежные значения p?

Question

Я использую lmerTest для получения значений p для моих моделей lmer в R. К сожалению, многократный запуск модели дает мне разные значения 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