モデルの仮定

前のコーディング演習でCox PHモデルを確認しましたね。この演習では、そのモデルがそもそも適切かどうかを確かめます。モデルは引き続きオブジェクトfitCPHに保存されています。

cox.zph()を使ってモデルfitCPHの比例ハザード仮定をチェックしてください。検定結果をtestCPHというオブジェクトに保存し、出力を表示します。
有意水準0.05で、1つの変数で仮定が破られているようです。どの変数でしょうか？その変数について、時間に対する係数betaのプロットを作成してください。
rmsパッケージのvalidate()関数で、クロスバリデーションを用いてモデルを検証してください。