如何用经验CDF绘制估计的CDF

Question

我正在拟合给定数据的分布。然后我估计了分布的参数。我还在R中使用ecdf（）命令绘制了经验cdf。现在我必须绘制估计分布的cdf和经验cdf 。我怎么能这样做？ ecdf（）命令可以帮助我吗？

Answer 1

是的，ecdf可以帮助您。它有一个绘图方法。下面是正态分布的可能代码。

x <- rnorm(100)
plot(ecdf(x))
lines(seq(-3, 3, by=.1), pnorm(seq(-3, 3, by=.1)), col=2)

编辑：你可以使用log-logistic分配函数来做同样的事情。它例如在包actuar中实现。

# load package
require(actuar)
# check out the parametrization! 
?dllogis
# estimate shape and scale from your data
shape <- 20
scale <- 1
# Don't do this. Use your own data instead. 
x <- rllogis(100, shape=shape, scale=scale)
# Plotting the empirical distribution function
plot(ecdf(x))
# x-values for plotting distribution function
xvals <- seq(min(x), max(x), length=100)
# plot estimated distribution function with your estimated shape and scale
lines(xvals, pllogis(xvals, shape=shape, scale=scale), col=2)