計算したいデータ (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 パーセンタイルを別々に使用したいので、これは私の場合には当てはまりません。