#r
#r
Вопрос:
Я построил модель со смешанным эффектом, используя lmer (), с целью сравнения роста показателей чтения для четырех разных групп детей по мере их старения. Я хотел бы построить график 4 разных наклонов с доверительными интервалами в R, чтобы визуализировать эту взаимосвязь, но я продолжаю застревать.
Я пытался использовать функцию plot и некоторые версии ggplot, как я делал для предыдущих моделей lm (), но пока это не работает. Вот моя попытка построения модели, которая, я надеюсь, показывает, как изменение показателей чтения с течением времени (возраст) взаимодействует с группой SESDLD ребенка (это указывает, есть ли у ребенка языковые проблемы и имеют ли они высокий или низкий доход).
AgeSES.model <- lmer(ReadingMeasure ~ Age.c*SESDLD1 (1|childid), data = reshapedomit, REML = FALSE)
Значение ReadingMeasure представляет собой непрерывную оценку, возраст.c — возраст по центру, измеряемый в месяцах. SESDLD1 — это категориальная мера, которая имеет 4 уровня. Я бы ожидал четырех положительных наклонов роста ReadingMeasure с разными перехватами и, вероятно, разными наклонами.
Я был бы очень признателен за любые указания о том, как это сделать!
Большое вам спасибо!!
Тип графика, которого я хотел бы достичь — это было сделано в Stata
Комментарии:
1. Добро пожаловать в Stack Overflow! Вы упомянули, что уже пробовали несколько вещей, не могли бы вы поделиться кодом, чтобы мы увидели, что именно вы пробовали? Не могли бы вы также поделиться (образцом) данных, которые вы используете?
2. Чаще всего я подхожу к такого рода задаче, получая прогнозы из модели, а затем строя график. Я написал об этом пост , в который я включаю пример из смешанной модели, соответствующей
lme(). Работа с lme4 была бы довольно похожей; основное отличие заключается вpredict.merMod()vspredict.lme().