R：[KDE]密度プロットをcdfに変換しますか？

Question

データ： 34,46,47,48,52,53,55,56,56,56,57,58,59,59,68

密度プロット

関数ecdf私がやりたいのは何

は％さんからを導出する導出密度プロットを取り、累積分布頻度にそれを回すあります。およびその逆。 私は、平滑化された累積分布関数を導出するために、特にカーネル密度推定を使用することを望んでいます。私は、ECDFを行うために生データポイントに頼るのではなく、KDEを使ってCDFを行うことを望んでいます。

編集：私はKernelSmoothing.CDFがある

参照、これは解決策になるかもしれませんか？そうだとすれば、これまでの実装方法はわかりません。

Mathworksには、ECDFからKECDFに変換して「計算し、評価されたcdfを特定の値のセットでプロットする」という私のやりたいことの例があります。

http://www.mathworks.com/help/stats/examples/nonparametric-estimates-of-cumulative-distribution-functions-and-their-inverses.html?requestedDomain=www.mathworks.com

私は実装がかなりずさんだと思うが。多項式回帰直線を考慮することは、よりよく適合する。

Answer 1

library("DiagTest3Grp", lib.loc='~/R/win-library/3.2") 

data <- c(34,46,47,48,52,53,55,56,56,56,57,58,59,59,68) 
bw <- BW.ref(data) 
x0 <- seq(0, 100, .1) 
KS.cdfvec <- Vectorize(KernelSmoothing.cdf, vectorize.args = "c0") 
x0.cdf <- KS.cdfvec(xx = data, c0 = x0, bw = bw) 
plot(x0, x0.cdf, type = "l") 

あなたはパッケージの生存 "にカプラン・マイヤー生存関数を見てすることができます私はまだ指定されたxのyを導出する方法を把握する必要がありますが、これは大きな助け