经验累积分布函数 标出四分位数，最大值，中值

大牛们好，

我想用R画出一组数的经验累积分布函数图，然后在该图上标出中值、最大值和 四分位数（75%），可以用程序自动实现吗？

多谢！！

原始数据比如：

pcdf=function(x){ plot(sort(x),1:length(x)/length(x),type="S",xlab="x",ylab="prob");
qq=quantile(x,c(0.25,0.5,0.75,1));abline(v=qq);
text(qq[1],0.05,paste(expression(q[25]), "=",format(qq[1], digits =3)),adj=c(1,1));
text(qq[2], 0.05,paste(expression(q[50]), "=",format(qq[2], digits =3)),adj=c(0.5,0));
text(qq[3], 0.05,paste(expression(q[75]), "=",format(qq[3], digits =3)), adj=c(0,1));
text(qq[4], 0.05,paste("max=",format(qq[4], digits =3)), adj=c(1,0))
}
看能不能满足你的要求，分位数和最大值都在图中标出。

楼上真棒啊！

学习了

1:length(x)/length(x)    这句话是不是 有问题啊， 不行啊

错误于xy.coords(x, y, xlabel, ylabel, log) :
  'x' and 'y' lengths differ
>

不知道怎么回事，不行啊，求大神

呵呵， 1楼是不是 不小心 把 （1:length(x)）/length(x)  写成了 1:length(x)/length(x)？

我下载不了你的数据，但是不是有quantile()这个函数吗，四个分为点都在，然后在经验函数图上标出来不就行了吗？经验函数图不是有plot(ecdf(x),verticals=T,do.p=F),还有中位数是median(x),都可以标出来啊

渲潇依旧 发表于 2012-11-25 11:15
我下载不了你的数据，但是不是有quantile()这个函数吗，四个分为点都在，然后在经验函数图上标出来不就行了 ...

不好意思，看成中位数了