r - Rで適合モデルの欠落を見つける

Question

coxphたとえば、100 個のデータ ポイントでモデルを適合させることを検討してください。95 のみが分析に含まれますが、5 は存在するNA(つまり、欠落している) ために除外されます。当てはめたデータの残差を抽出するので、95 個の観測値を持つ残差ベクトルが得られます。残差を元のデータ フレームに戻したいのですが、長さが異なるため実行できません。

元のデータ フレームのどの観測値がモデルに含まれていないかを特定するにはどうすればよいですか? そのため、それらを除外/削除して 2 つの長さを同じにすることができます。

(元のデータの方がはるかに大きいため、データが欠落している場所を特定するのは困難です...)

Answer 1

na.action引数をに設定して、モデルを再調整しますna.exclude。これにより、近似オブジェクトの一部である残差と近似値がNAs でパディングされます。元のモデルが次の場合zn50:

zn50_na <- update(zn50, na.action=na.exclude)

これにより、適切な長さが得られresiduals(zn50_na)ますfitted(zn50_na)。詳細については、を参照?na.omitしてください。