R、 python还是octave：经验分位数（逆cdf）与置信区间？

Question

我在寻找一个内建函数，它返回样本分位数和估计的置信区间（MATLAB的ecdf会这样做）。在

我猜R有这个内置的，但我还没有找到它。在

如果您有任何独立的代码来执行此操作，您也可以在这里指向它，尽管我希望找到一些作为更大的开放代码库的一部分包含的内容。在

-试图摆脱MATLAB。在

Answer 1

利用survfit函数可以得到具有置信区间的生存函数。因为它只是1-ecdf，分位数之间有直接的关系。要使用此功能，您必须创建一个变量，该变量表示您的每个观测值都是完整的（不是经过审查的）：

library(survival)
x <- rexp(10)
ev <- rep(1, length(x))
sf <- survfit(Surv(x,ev)~1)

带输出：

^{pr2}$

实际上，survfit确实计算了中值及其置信区间，但不计算其他分位数：

>sf
Call: survfit(formula = Surv(x, ev) ~ 1)

records   n.max n.start  events  median 0.95LCL 0.95UCL 
 10.000  10.000  10.000  10.000  -0.205  -0.940      NA 

计算中值置信区间的实际工作很好地隐藏在survival:::survmean函数中，您可以使用它来推广到其他分位数。在