Я относительно новичок в R и, очевидно, не очень опытный. Однако я использовал многоуровневое моделирование для определения влияния голоса на параметры сна. Например, TST в этом фрагменте - это общее время сна, интенсивность - это интенсивность голоса (в данном случае как среднее).
Мне удается получить диаграмму рассеяния, в зависимости от номера участника, как я хочу. Однако теперь я хочу включить линии регрессии для моей модели, которые показывают пересечение и наклон для модели против нулевой модели (исключая мою независимую переменную).
Тем не менее, независимо от того, что я пытаюсь, я не могу отобразить линии регрессии, основанные на перехвате и наклоне, даже после ввода их значений вручную!
Вот мой код для расчет и мой расчет и для участка.
Модель:
library(lme4)
TST_RE_Intensity.model = lmer(Intensity_mean ~ TST_re + Day + (1+ TST_re|Participant_ID) + ( 1|Filename), data=my.df, REML = FALSE)
TST_RE_Intensity.null = lmer(Intensity_mean ~ Day + (1+ TST_re|Participant_ID) + (1|Filename), data=my.df, REML = FALSE)
Участок:
library(ggplot2)
p <- ggplot(my.df, aes(x=my.df$TST_re, y=my.df$Intensity_mean, colour=my.df$Participant_ID))+
theme(legend.position = "none")+
geom_point(shape=20) +
geom_abline(aes(intercept=64, slope = - 0.0167, size= 1.5)+
geom_abline(aes(intercept=61, slope = - 0.0162, size=1.5)+
scale_size_manual(values = c(0.3, 0.3))+
scale_y_log10(name="Log10(TST)", limits=c(40,80)) +
scale_x_log10(name="Log10(Intensity)")
Я не получаю никаких сообщений об ошибках, но также никогда не вижу любые строки. Я пытался следовать инструкциям руководства ggplot2, как это, но в итоге получился короткий
p + geom_abline(intercept = 37, slope = -5)
Есть ли способ просто нарисовать линию "вручную"?
Заранее спасибо !!