如何獲取累積分佈函數

Question

我有一個隨機變量，我可以得到的

density(r)->pdf 

如何獲得從該PDF與R中的累積分佈函數（CDF）的密度ρ

我試圖

approxfun(pdf$x,pdf$y,yleft=0,yright=0)->f 
cdf<-integrate(f,-Inf,2) 

但是當我繪製CDF我只獲得

1.000396絕對錯誤< 7.9e-05

Answer 1

您可以使用R的logspline包，用於創建和可視化平滑的非參數累積分佈函數（和其他相關量）如下。

library(logspline) # load the package 
x <- rnorm(500) # generate some data 
x.lsp <- logspline(x) # create the logspline object 
plot(x.lsp, what="p") # plot the distribution function 
plot(x.lsp, what="d") # plot the density function