r - 如何用经验 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)

编辑：您可以使用对数逻辑分布函数来做同样的事情。例如，它在 package 中实现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)