我正在尝试为纵向数据集中的不同组的原型个体创建固定效果的ggplot2图(标明种族和受教育年限对工资增长的影响)。我没有使用数据集本身,而是尝试根据混合效应模型的截距和斜率系数创建线形图。因此,我创建了一个空白框架,并向其添加了ablines。但是,我想在图例中添加说明这些线条的图例。这是我用来创建情节的代码。
df <- data.frame()
fixef.c <- list(intercept0 = 1.72, slope0 = 0.04, intercept1 = 0.038, slope1 = -0.016)
ggplot(df) + xlim(0, 12) + ylim(1.6, 2.4) +
geom_abline(intercept = fixef.c[[1]], slope = fixef.c[[2]]) +
geom_abline(intercept = fixef.c[[1]] + fixef.c[[3]]*3, slope = fixef.c[[2]], linetype = 2) +
geom_abline(intercept = fixef.c[[1]], slope = fixef.c[[2]] + fixef.c[[4]], linetype = 3) +
geom_abline(intercept = fixef.c[[1]] + fixef.c[[3]]*3, slope = fixef.c[[4]] + fixef.c[[2]], linetype = 4) +
xlab("Log Wages") +
ylab("Years Experience")
我需要知道如何在此图上添加自定义图例来解释线条。理想情况下,我希望将图例放置在相同的位置,并使用与使用基础图形包创建的下图相似的标签。(注意:我尝试添加show.legend = TRUE
到每个Abline函数中,但没有用。
exper.seq <- seq(0, 12)
x.w9 <- fixef.c[[1]] + fixef.c[[2]]*exper.seq
x.w12 <- fixef.c[[1]] + fixef.c[[2]]*exper.seq + fixef.c[[3]]*3
x.b9 <- fixef.c[[1]] + fixef.c[[2]]*exper.seq + fixef.c[[4]]*exper.seq
x.b12 <- fixef.c[[1]] + fixef.c[[2]]*exper.seq + fixef.c[[3]]*3 +
fixef.c[[4]]*exper.seq
plot(exper.seq, x.w9, ylim=c(1.6, 2.4), ylab="LNW.hat", xlab="EXPER", type="l", lwd=2)
lines(exper.seq, x.w12, lty=3)
lines(exper.seq, x.b9, lty=4, lwd=2)
lines(exper.seq, x.b12, lty=5)
legend(0, 2.4, c("9th grade, White/Latino", "9th grade, Black",
"12th grade, White/Latino", "12th grade, Black"), lty=c(1, 4, 3, 5))
如果您可以允许创建一个包含截距和坡度以及分组的对象,那么这是一项相当简单的任务。
xy <- as.data.frame(
rbind(line1 = c(intercept = fixef.c[[1]], slope = fixef.c[[2]]),
line2 = c(intercept = fixef.c[[1]] + fixef.c[[3]]*3, slope = fixef.c[[2]]),
line3 = c(fixef.c[[1]], slope = fixef.c[[2]] + fixef.c[[4]]),
line4 = c(fixef.c[[1]] + fixef.c[[3]]*3, slope = fixef.c[[4]] + fixef.c[[2]]))
)
xy$design <- rownames(xy)
ggplot() +
xlim(0, 12) + ylim(1.6, 2.4) +
geom_abline(data = xy, aes(intercept = intercept, slope = slope, linetype = design))
您还可以修改哪些线型属于哪个级别。
scale_linetype_manual(values = c(line1 = "solid", line2 = "dashed", line3 = "dotted", line4 = "dotdash")) +
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句