たとえば、次のような CDF を持つ離散関数があります。
cdf <- c(0.00, 0.35, 0.71, 0.92, 1.00, 1.00, 1.00, 1.00)
次の行で一種の分位関数を作成できます。. .
result <- which(cdf == min(cdf[cdf > x]))
. . . ここで、x は累積確率です。たとえば、qfunction(0.9) = 4 および qfunction(0.99) = 5 です。
このソリューションは、ベクトルを処理するまでは(エレガントではありませんが)問題ないように見えます。したがって、x = c(0.9, 0.99) の場合、関数は失敗します。これは、人々がRで多くのことを行うように思えますが、解決策はまだ見つかりません. R は私の第一言語ではありません。
どんな助けでも大歓迎です。