标题说明了一切。我想做的是绘制中位数,但也显示每个时间点的两个四分位数和第90个百分点。
答案 0 :(得分:1)
首先,这是一个小样本数据集,其中包含3个时间点y
的观察time
(每个时间点有100个观测值)。
dd< - data.frame(y = c(rnorm(100,10,2),rnorm(100,15,2),rnorm(100,20,2)), time = rep(c(1,2,3),每个= 100))
接下来,根据您的需要总结您的数据(这里只有中位数和第一个分位数,但原理是相同的,以获得其他分位数)
dd1< - ddply(dd," time",总结, MED =中间值(Y), firstquart = quantile(y,probs = .01))
最后,用ggplot
绘制线条 - 第一个分位数为红色。
ggplot(dd1,aes(x = time,y = med))+ geom_line()+ geom_line(aes(x = time,y = firstquart),color =" red")