JMP（SAS）與R在中位生存時間估計中存在奇怪的不一致嗨，大家好，
我在中位生存時間估計中發現了R和JMP（或SAS）之間的以下奇怪的不一致。
我知道CI可能因為不同的方式而有所不同，但我從未預料到中位生存點估計會有所不同。 有沒有人碰巧知道實際的技術原因？我想弄清楚，但沒有運氣。
對於以下簡單數據集（事件= 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 =事件，檢查代碼= 0）
R輸出:(代碼：庫（生存）; 摘要（倖存（Surv（時間，事件）~trt，data = data））$ table ）
| | ... | 中位數 | ... |
| TRT = 0 | ... | 116 | ... |
| TRT = 1 | ... | 150 | ... |
謝謝大家！
Re：JMP（SAS）與R在中位生存時間估計中的奇怪不一致非常感謝您的解釋！那講得通！
現在我必須看看我是否可以找到一個公式或者某個在SAS vs R中給出明確定義的東西，因為我需要讓R輸出SAS結果。
再次感謝！
最好，
嚴
Re：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個觀測值的值為中位生存時間。請注意，如果再更改一次t = 150觀察（第10行）進行右刪失，R會報告trt = 1的中位數缺失（NA）。
希望這可以幫助，
邁克爾
