R语言 ggplot2 画平滑图

最新推荐文章于 2024-10-02 15:13:13 发布

转载最新推荐文章于 2024-10-02 15:13:13 发布 · 1k 阅读

0 ·

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/emanlee/p/6308953.html

文章标签：

#r语言

本文通过R语言的ggplot2包演示了如何基于特定数据集绘制柱状图并添加平滑曲线。利用ns()函数进行样条平滑，并采用对数链接的广义线性模型(GLM)来拟合数据。

library(splines)

library(ggplot2)

dt1 <- structure(list(Age = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("o80", "u80"), class = "factor"), NoP = c(47L, 5L, 33L, 98L, 287L, 543L, 516L, 222L, 67L, 14L, 13L, 30L, 1L, 6L, 17L, 30L, 116L, 390L, 612L, 451L, 146L, 52L), pctOAC = c(0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100, 0, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100)), .Names = c("Age", "NoP", "pctOAC"), row.names = c(NA, -22L), class = "data.frame")

ggplot(dt1, aes(x=pctOAC,y=NoP, fill=Age)) +
    geom_bar(stat="identity",position=position_dodge()) +
    geom_smooth(aes(colour=Age), se=F,
                method="glm",
                formula=y~ns(x,8),
                family=gaussian(link="log"),
                show_guide = FALSE,lwd=0.7) +
    theme(legend.position=c(.2,0.8))

##Ref: http://stackoverflow.com/questions/13213451/smoothing-in-ggplot

转载于:https://www.cnblogs.com/emanlee/p/6308953.html