我有要进行一些计算的数据(egd 和 di)。
数据:
d= 1:20
xi= 2:21
我想使用这个等式(对 d 和 xi 分别使用第 95 个和第 5 个百分位数):
res= d(95th)-d(5th)/xi(95th)+xi(5th)
res2= res + xi(5th)
我们使用这个函数得到第 95 和第 5 个百分位数(由 stackoverflow 的用户提供):
fun <- function(x){
quantiles <- quantile( x, c(.05, .95 ) )
x[ x < quantiles[1] ] <- quantiles[1]
x[ x > quantiles[2] ] <- quantiles[2]
x
}
但这不适用于我的情况,因为我想在上面的等式中分别使用第 95 个和第 5 个百分位数。请有任何想法。