Не удается заставить вычисление p-значения kr в tab_model (sjPlot) работать

Я пытаюсь обобщить модель lmer с помощью sjPlot::tab_model. Это работает, пока я не попытаюсь получить p-значения «kr», чтобы также получить степени свободы. У меня установлен pbkrtest, и я могу использовать его отдельно, но он безуспешно вызывается tab_model. Кто-нибудь знает, что мне не хватает?

Я запустил линейную модель смешанных эффектов. Я могу создать таблицу, используя

tab_model(BWRTmod1b,show.stat = TRUE)

Когда я изменяю вызов, чтобы получить p-значения на основе аппроксимированных степеней свободы Кенворда-Роджера, например:

tab_model(BWRTmod1b,p.val= "kr" ,show.stat = TRUE) 

Я получаю следующее сообщение об ошибке:

Ошибка в раунде (attr (pv, "df.kr", точно = TRUE)) : нечисловой аргумент математической функции

Я погуглил проблему, и единственное, что я нашел, это отсутствующий пакет pbkrtest, но здесь это не так. Я также погуглил само сообщение об ошибке и обнаружил, что это может произойти, когда фрейм данных передается в функцию, которая требует числового ввода. Я также обновил R и переустановил/обновил все пакеты, но это не помогло.

Я не знаю, могу ли я это исправить и как, и вполне возможно, что я упускаю что-то еще. Любые подсказки будут оценены.

Ссылка на выходные данные модели

# working
tab_model(BWRTmod1b,show.stat = TRUE)

# not working
tab_model(BWRTmod1b,p.val= "kr" ,show.stat = TRUE) 

# long term goal
tab_model(BWRTmod1b,p.val= "kr", show.df=TRUE ,show.stat = TRUE) 

Я ожидал создать таблицу результатов моей линейной модели смешанных эффектов, которая включает степени свободы. Я могу сгенерировать сводную таблицу на основе приближения Вальда, но не могу получить версию Кенворда-Роджера со степенями свободы.