生存期間中央値推定におけるJMP（SAS）とRの間の奇妙な矛盾皆さんこんにちは、
生存期間中央値の推定において、RとJMP（またはSAS）の間に次のような奇妙な矛盾がありました。
CIはそれを実行する方法が異なるために異なる可能性があることを理解していますが、生存期間の中央値の推定値が異なる可能性があるとは思っていませんでした。 実際に技術的な理由を知っている人はいますか？私はそれを理解しようとしましたが、運がありませんでした。
次の単純なデータセット（event = 1は死亡、0は検閲を意味します）、
時間 | イベント | trt |
116 | 1 | 0 |
116 | 0 | 0 |
110 | 1 | 0 |
116 | 1 | 0 |
116 | 1 | 0 |
110 | 1 | 0 |
116 | 1 | 0 |
116 | 1 | 0 |
150 | 0 | 1 |
150 | 1 | 1 |
138 | 1 | 1 |
138 | 1 | 1 |
150 | 0 | 1 |
150 | 0 | 1 |
150 | 0 | 1 |
138 | 1 | 1 |
Rはそうであるが、JMPはtrt = 1群の生存期間中央値の推定値を提供することができない。
JMP出力：（JMP 12.1、分析 - >信頼性と生存 - >生存、Y =時間、グループ化= trt、cencor = event、censor code = 0を使用）
| グループ | 時間の中央値 | ... |
| 0 | 116 | ... |
| 1 | 。 | ... |
Rの出力：（コード：ライブラリ（生存）。概要（survfit（Surv（時間、イベント）〜trt、データ=データ））$ table ）
| | ... | 中央値 | ... |
| trt = 0 | ... | 116 | ... |
| trt = 1 | ... | 150 | ... |
よろしくお願いします。
生存期間中央値推定におけるJMP（SAS）とRの間の奇妙な矛盾説明をどうもありがとう！それは理にかなっている！
RにSASの結果を出力させる必要があるため、SAS vs Rで明確な定義を与える式または何かを見つけることができるかどうかを確認する必要があります。
再度、感謝します！
ベスト、
ヤン
生存期間中央値推定におけるJMP（SAS）とRの間の奇妙な矛盾現在起こっているのは、中央値を定義するには右打ち切り観測が多すぎるということです。trt = 1グループの8つの観測に対して、3つはt = 138での失敗、1つはt = 150での失敗、そしてそれらの半分（4）はt = 150での右打ち切りです。8つの値の中央値は4番目と5番目の値の平均なので、150と150よりも大きい（未知の）値を平均する必要があります。その結果、サバイバルプラットフォームは中央値の欠損値を報告しています。
なぜRが値を与えるかについては、私は推測することしかできません、しかし私の推測は4番目の順序付き値が打ち切られず、8のうち4番目の観測値が0.50分位数になるので、Rは4番目の観測値を生存期間中央値。もう1つt = 150の観測値（行10）を右打ち切りに変更すると、Rはtrt = 1の中央値が欠落している（NA）と報告します。
お役に立てれば、
マイケル
