r - Cox モデルの欠落の理由を理解する

Question

次の Cox モデルを実行したところ、1526679 個の削除された観測が得られました。これは私のデータの大部分です。

Call: coxph(formula = Surv(time1sec, time2sec, event) ~ gain + 
Buy + Lev + TP + frailty(ID), data)

n= 73322, number of events= 73322 (1526679 observations deleted due to missingness)

これらの観測が削除された理由はわかりません。これらの値が存在し、空ではないことは確かです。これは、脆弱性用語として ID を追加したときに発生し始めました。

ここで何が起こっているのでしょうか？

Answer 1

標準的な回帰 (および cox 回帰の場合の回帰風) タイプの問題では、欠損値を処理するためのデフォルトの方法は単純にそれらを無視することです。これは、予測しようとしている係数と値の両方で発生します。

必要なものがすべてそこにあることを明示的に確認することから始めます。これは、次のコードで実行できます。

apply(data, 2, function(x) length(which(is.na(x))))

また、次のコードを使用して、Surv オブジェクトに NA がないことを確認します。

length(which(is.na(Surv(time1sec, time2sec, event))))