密度関数から確率を計算する

Question

密度関数を構築しましたが、新しいデータポイントが選択された間隔（たとえばa = 3、b = 7）に落ちる確率を計算したいと思います。だから、私は探しています：

P(a<x<=b) 

いくつかのサンプルデータ：

df<- data.frame(x=c(sample(6:9, 50, replace=TRUE), sample(18:23, 25, replace=TRUE))) 

dens<- density(df$x) 

は、私はすべてのソリューションを聞いて幸せになるでしょうが、好ましくはベースR に

はありがとう

Answer 1

濃度を関数として取得し（approxfunを使用）、必要な制限を超えて関数を積分する必要があります。

integrate(approxfun(dens), lower=3, upper=7) 
0.258064 with absolute error < 3.7e-05 

## Consistency check 
integrate(approxfun(dens), lower=0, upper=30) 
0.9996092 with absolute error < 1.8e-05