Нанесите маржинальные эффекты в исходном масштабе, используя plot_model()

Я пытаюсь построить предельные эффекты из модели (сделано с помощью lme4), используя plot_model()

Модели (2 версии одной модели) выглядят как

Версия 1:

m1 <- glmer(Y~Xs*Zs+(1|random.factor), family=binomial(), data=dat)

Версия 2:

m2 <- glmer(Y~scale(X)*scale(Z)+(1|random.factor), family=binomial(), data=dat)

В обеих версиях непрерывные предикторные переменные X и Z масштабируются с использованием scale() функция, но в версии 1 это делается до создания модели m1 и хранится в виде отдельных столбцов (Xs и Zs) в dataframe=datтогда как в версии 2 это указано в модели m2

Мой вопрос сейчас: если я строю график m1, как я получаю

Участок без масштабированной оси абсцисс

Я хотел бы построить предельные эффекты для X, но вместо масштабированных значений для X на оси X я хотел бы иметь оригинальные, которые я пробовал

a <- plot_model(m1, type="pred", terms="Xs [myfun]")


b <- plot_model(m2, type="pred", terms="X [myfun]")

myfun была определена как обратная функция scale()

Участок

Участок б

Итак, это дает мне два разных сюжета. Поскольку график b выглядит аналогично графику m1 без масштабированной оси x, я бы сказал, что график b является правильным. Но я не понимаю, как "[myfun]" в термине ="..." действительно работает и что он делает. Спасибо за любую помощь!

0 ответов

Другие вопросы по тегам