R中取百分位比用quantile()
1. 求某个百分位比
> x<-rnorm(200)
> quantile(x,0.9)
90%
1.377796
> quantile(x,c(0.1,0.9))
10% 90%
-1.317574 1.377796
2. 产生一序列百分位比值
> quantile(x,seq(0.1,1,0.1))
10% 20% 30% 40% 50% 60% 70% 80% 90%
-1.31757448 -0.94861519 -0.67421389 -0.33618527 -0.03578113 0.24425065 0.50656752 0.99848083 1.37779613
100%
3.83505973
2. 只取数值:unname()
> y<-quantile(x,0.9)
> y
90%
1.377796
> y-1
90%
0.3777961
会带着百分值,只取
举例:
upbound.p[i]<-unname(quantile(x$spread_positive[(i-n):(i-1)],upbound,na.rm=T))
> y<-unname(quantile(x,0.9))
> y
[1] 1.377796