私は、3 種類のサービス (酒場、レストラン、急行) の 10 年間の調査における死亡率を表すこの生存データを持っています。
データには、サービスの種類 (1 = サロン、2 = レストラン、3 = 特急)、年 (1 から 11 までの整数、11 は 10 年以上を意味する)、および検閲の 3 つの変数が含まれています。
2 つの質問があります。
1) Cox 比例ハザード モデルを当てはめましたが、比例ハザードの仮定を確認する方法は何ですか。つまり、各個人のハザード比とベースライン ハザードは時間に依存しないと仮定します。
2) 時間依存の Cox モデルを R に適合させる方法は?
ここに私のコードがあります:
#Cox Proportional Hazards
cox <- coxph(Surv(Years, Censor) ~ data$`Service Type` )
summary(cox)