にサバイバル オブジェクトがありR
ます。
print(surv)
私にくれます
> print(surv)
Call: survfit(formula = Surv(TAGE, EVENT) ~ 1, data = data_LTC[data_LTC$TYPE ==
"Job", ])
records n.max n.start events median 0.95LCL 0.95UCL
299510 299510 299510 252884 177 173 180
ただし、quantile(surv)
機能せず、出力されます
> quantile(surv)
Error in is.na(y) : 'y' is missing
R
は中央値 (177) を計算できますが、他の四分位数は計算できないため、私にとっては奇妙なことです。
サバイバル オブジェクトの何が問題になっていますか?
[編集]
おそらくこれは問題の最小の例です: ドキュメント ( http://stat.ethz.ch/R-manual/R-patched/library/survival/html/quantile.survfit.html ) には、この例があります:
> fit <- survfit(Surv(time, status) ~ ph.ecog, data=lung)
> quantile(fit)
$quantile
25% 50% 75%
ph.ecog=0 285 394 655
ph.ecog=1 181 306 550
ph.ecog=2 105 199 351
ph.ecog=3 118 118 118
...
さて、この出力の最初の行だけを繰り返したい場合は、
> fit <- survfit(Surv(time, status) ~ 1, data = lung[lung$ph.ecog == 0,])
> quantile(fit)
Error in is.na(y) : 'y' is missing
@Edwinは、quantile(fit$time)
代わりに使用することを以下に提案しました
> quantile(fit$time)
0% 25% 50% 75% 100%
5.00 224.25 301.50 439.75 1010.00
ただし、明らかに異なる結果につながります。
[終了]以下の回答はパッケージ内で
使用されていないため無視してください。ただし、組み込みの関数を使用しています。quantile.survfit
survival
R
quantile
survival
この問題を解決するには、パッケージの最新バージョンに更新してください。
これを行う
update.packages()
そのためには root 権限が必要な場合があることに注意してください。