R: Поверните график плотности [KDE] в cdf?

Question

данных: 34,46,47,48,52,53,55,56,56,56,57,58,59,59,68

Плотность участка

ECDF

То, что я хотел бы сделать, это производный график плотности и превратить его в кумулятивную частоту распределения, чтобы получить% 's из. И наоборот. Моя надежда заключается в том, чтобы использовать оценку плотности ядра специально для получения сглаженной функции кумулятивного распределения. Я не хочу полагаться на исходные точки данных, чтобы сделать ECDF, но использовать KDE для создания CDF.

Edit:

Я вижу, есть 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") 

мне еще нужно, чтобы выяснить, как получить данные Y х, но это было серьезным подспорьем