在計算生存函數相等性的對數秩檢驗時，R 的生存 package 如何處理發生在 t0 的死亡？

Question

我在 Stata 中讀到，存在一個問題，即觀察可以在 t0 結束，因此無意中從分析中省略了。 解決方法是在 t0 時為每次死亡添加一些小的 e > 0 以防止這種情況發生。 有關 Stata 中問題的描述，請參閱此鏈接。 https://www.stata.com/support/faqs/statistics/time-and-cox-model/

https://www.statalist.org/forums/forum/general-stata-discussion/general/1609582-cox-regression-observations-end-on-or-before-enter

我的問題是 R 在運行長排名比較時如何使用生存 package 來處理這個問題？ 我找不到任何相關文檔，但通過運行一些測試，似乎 R 會自動添加一點點 epsilon >0 來解決這個問題。 任何人都可以確認或否認這一點，如果可以的話，提供一種方法來覆蓋默認行為？

我嘗試在 R 和 Stata 中運行並行分析，但得到了不同的結果。 Stat 給出了信息“

 45,024  total observations
          1  event time missing (GTIME>=.)                      PROBABLE ERROR
        341  observations end on or before enter()

“

而R則沒有給出任何這樣的信息。 卡方統計計算的最終結果不同，我認為這是因為這種差異。 任何人都可以確認嗎？

Answer 1

沒關系，找到了解決方案。 R 保留在 t = 0 時消失的觀察結果，而 Stata 則忽略它們。