Значения p в lmer с использованием теста lmer. Запускать модель несколько раз, чтобы получить надежные значения p?
Я использую 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