2

MATLAB 以外のサンプル分位点と推定信頼区間を返す組み込み関数を探しています(MATLABecdfはこれを行います)。

Rにはこれが組み込まれていると思いますが、まだ見つけていません。

これを行うためのスタンドアロン コードがある場合は、ここでそれを参照することもできますが、より大きなオープン コード ベースの一部として含まれているものを見つけたいと考えています。

-MATLAB から離れようとしています。

4

1 に答える 1

4

このsurvfit関数を使用して、信頼区間で生存関数を取得できます。これは単なる 1-ecdf であるため、分位数の間には直接的な関係があります。これを使用するには、各観測が完全である (打ち切られていない) ことを示す変数を作成する必要があります。

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

出力あり:

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

     time n.risk n.event survival std.err lower 95% CI upper 95% CI
-1.4143     10       1      0.9  0.0949       0.7320        1.000
-1.1229      9       1      0.8  0.1265       0.5868        1.000
-0.9396      8       1      0.7  0.1449       0.4665        1.000
-0.4413      7       1      0.6  0.1549       0.3617        0.995
-0.2408      6       1      0.5  0.1581       0.2690        0.929
-0.1698      5       1      0.4  0.1549       0.1872        0.855
 0.0613      4       1      0.3  0.1449       0.1164        0.773
 0.1983      3       1      0.2  0.1265       0.0579        0.691
 0.5199      2       1      0.1  0.0949       0.0156        0.642
 0.8067      1       1      0.0     NaN           NA           NA

実際、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関数にうまく隠されているため、おそらく他の分位数に一般化するために使用できます。

于 2010-08-09T20:30:07.540 に答える