0

我对 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)

有没有办法“手动”绘制一条线？

提前致谢！！

4

0 回答 0