生存関数グラフの重ね合わせ

[質問]

カテゴリカル変数を含む、比例ハザードモデルをPHREGプロシジャで推定しています。カテゴリカル変数の各水準における生存関数のODS統計グラフを一つにまとめて表示できますか。

[回答]

カテゴリカル変数を含め、説明変数の任意の値における生存関数を求めるには、BASELINEステートメントを用います。この際、説明変数の値を含むデータセットを予め作成します。

例)

                      /* 説明変数の値の組合せデータセットを作成 */
DATA cov;
  INPUT id$ drug$ x;
  DATALINES;
NAME1 A  10
NAME2 B  10
;
RUN;

この例では、説明変数として変数DRUG、Xがあり、変数DRUGの水準として、A, B があると仮定しています。変数IDは、組み合わせに対する名称に該当します。

ODS統計グラフを作成するためには、PROC PHREGステートメントにてPLOTS=オプションを指定します。生存関数のグラフの場合、オプションの値として SURVIVAL を用い、一つのグラフとして表示するため、OVERLAYを追記します。

例)

                                            /* PLOTS=オプションを追加 */
PROC PHREG DATA=test PLOTS(OVERLAY)=(SURVIVAL);
  CLASS drug;
  MODEL time*censor(1)=drug x;
                                       /* BASELINEステートメントを追加 */
  BASELINE COVARIATES=cov / ROWID=id;
RUN;