Как извлечь значения наклона и перехвата для разных групп из взаимодействующих участков в jtools при построении линейных моделей с фиксированными эффектами - PullRequest
Новая
       11

Как извлечь значения наклона и перехвата для разных групп из взаимодействующих участков в jtools при построении линейных моделей с фиксированными эффектами

0 голосов
/ 05 февраля 2019

Я пытаюсь извлечь наклон и перехват для каждой из моих групп из моих линейных моделей смешанных эффектов.Модель была построена с использованием lmer в библиотеке lme4, и я могу просмотреть результаты для каждой группы, используя interactive_plot из библиотеки jtools.Как мне получить наклон и перехват для каждой из этих линий?

Я знаю, что могу использовать summary () или summ (), чтобы увидеть оценки для фиксированных эффектов и дисперсий случайных эффектов, но я не могу увидеть оценки случайных эффектов.Поэтому я не могу точно рассчитать наклон и точки пересечения моделей.

>library(lme4)
> cond_waterxsilver <- lmer(LnAg ~ LnVolume + (LnVolume | FilterID) + SilverType + WaterType + SilverType*WaterType + SilverType*LnVolume +  WaterType*LnVolume, data=capwater_removed.data)


> library(jtools)
> interact_plot(cond_ranin_waterxsilver, pred = LnVolume, modx = WaterType, mod2 = SilverType)

Я просто пытаюсь получить наклон и точки пересечения для шести линий модели (два разных WaterType и три разных SilverType).Есть ли в jtools инструмент или другой пакет, который может помочь мне извлечь наклон и перехваты из моей модели?

1 Ответ

0 голосов
/ 15 февраля 2019

Я разработчик этого пакета!

Краткое примечание: эта и другие функции, о которых я собираюсь упомянуть, были только что перемещены в новый пакет, который называется interactions, который находится в процессе добавления вКРАН.Предполагая, что вы не обновились до последней версии jtools (2.0.0; только что вышла несколько дней назад), эти функции все еще доступны в пакете jtools.Если вы обновитесь до jtools 2.0.0, вам нужно будет перейти по по этой ссылке для получения инструкций о том, как загрузить interactions, прежде чем он попадет в CRAN.

Там должно бытьпростой ответ на ваш вопрос.Функция sim_slopes (сокращение от «простые уклоны») должна дать вам то, что вы ищете. 

sim_slopes(cond_ranin_waterxsilver, pred = LnVolume, modx = WaterType, mod2 = SilverType, cond.int = TRUE)

Это распечатает условные уклоны и перехваты (перехваты печатаются только тогда, когда cond.int = TRUE.

Если вам нужно запрограммировать эти значения, вы можете сохранить объект sim_slopes. 

ss <- sim_slopes(cond_ranin_waterxsilver, pred = LnVolume, modx = WaterType, mod2 = SilverType, cond.int = TRUE)

ss$slopes # Matrix of slopes with test statistics, etc.
ss$ints # Matrix of intercepts with test statistics, etc.
...