表形式とグラフビルダーを使う利点の1つは、データフィルターに含まれる観測データの母集団に合わせて、分析結果を動的に更新できることです。母集団が変更されると要約結果がどのように変化するかを見るのは興味深いものです。特定の特性が結果にどのように影響するかについての洞察が得られます。
ただし、グラフの動的更新のメリットから除外されるエンドポイントのセットが1つあります。それは、イベント発生までの時間エンドポイントです。Survivalプラットフォームから分析を生成し、推定値をデータテーブルに保存し、グラフでカプラン・マイヤープロットを生成することは可能ですが、 サバイバルから取得した推定値の基となるデータテーブルを最初に更新しない限り、このプロットは元のデータテーブルの観測値のフィルタリングに反応しません。これによりループが発生します。このブログでは、患者について話しているものとします。
- Analyze patients using サバイバル and output estimates.
- Produce Kaplan-Meier plot from estimates using グラフビルダー.
- Filter patients in original table using データフィルター.
- Update analysis in サバイバル and output new estimates.
- Rebuild plot from updates estimates using グラフビルダー.
カプラン・マイヤー曲線の継続的な更新は、シームレスかつ容易に行えるべきです。基礎となる推定値に合わせて更新されるリスク数表を追加することも重要です。なお、リスク数表は、多くの医学雑誌(このページで「生存曲線」を検索してください)や米国食品医薬品局(FDA)において、カプラン・マイヤー曲線に必須の要件となっています。
サバイバルエクスプローラー アドインを入力してください。
サバイバルエクスプローラーアドインは、JMPサバイバルプラットフォームを使用してイベント発生までの時間分析を実行し、 GraphPadを使用してリスク数表付きのKaplan-Meier曲線を作成します。統計検定と信頼区間を使用して曲線を修正および注釈付けするオプションがあります。最も重要なのは、このアドインを使用すると、付属のデータフィルターに共変量を追加できるため、関心のある観測値にサブセット化することで分析がどのように動的に変化するかを把握できることです。
Survival Explorerのダイアログを図1に示す。
Figure 1. Dialog for Survival Explorer図1. サバイバルエクスプローラーのダイアログ
役割の説明:
- Y、イベント発生までの時間:各観測値について、イベント発生または打ち切りまでの時間を示す数値列。
- グループ化:オプションの文字列で、グループ化列の各レベルごとに個別のイベント発生時間曲線を作成するために使用されます。
- 検閲:各観測値について、イベントまたは検閲結果が発生したかどうかを示す数値列。
選択肢は以下のとおりです。
- 検閲コード:検閲列の値で、観測値が検閲されているかどうかを示します。その他の値は、観測値に何らかの事象が発生したことを示します。
- リスク人数間隔:リスク人数表は、この間隔に基づいて値が表示されます。例えば、間隔が25の場合、リスク人数表の値は25日ごとに表示されます(時間列がその単位で測定されていると仮定します)。間隔がすべてのグループレベルで最大時間を超える場合、最後のエントリは最大時間に表示されます。
- グループ別リスク表の色分け:このオプションを選択し、グループ分け列を指定すると、リスク表の行がカプラン・マイヤー曲線の線の色に合わせて色分けされます。選択されていない場合、またはグループ分け列が指定されていない場合は、リスク表の値は黒色になります。このオプションを使用するには、グループ分け列に対して値の色を表す列プロパティが定義されている必要があります。
- リスクのある人数の相対的な大きさ:この値は、カプラン・マイヤー曲線におけるリスクのある人数の相対的な大きさを調整します。グループ分けのレベル数が多い場合は、図の下部にレベル数を表示するのに十分なスペースを確保するために、この値を大きくする必要があります。
- リスク数値オフセット:この値は、リスク数値に表示されるデータの位置を調整し、より見やすくします。
本アドインの使用方法を、Collett (2023) のデータセット「肝移植レシピエントの生存率」を用いて説明します。このデータセットには、移植片不全までの日数、不全の原因、年齢、性別、および原疾患(ALD:アルコール性肝疾患、PSC:原発性硬化性胆管炎、PBC:原発性胆汁性胆管炎)に関するデータを持つ1,761人の患者が含まれています。
図1に示すダイアログは、我々の分析を表しています。「失敗」列は患者が治療失敗であるかどうかを示しており、そのため「打ち切りコード」は0に設定されています。最大期間は4,720日(約13年!)であるため、リスク期間の人数は365に設定されています。その他のオプションにはデフォルト値が使用されます。
Figure 2 presents the results of the analysis from the サバイバル platform after adding gender and age to the データフィルター. The Add Simultaneous Intervals option is switched on.
Figure 2. Analysis of liver transplant data図2.肝移植データの解析
出力画面のオプション領域には、最大4つのボタンが表示されます。ボタンの名前は選択すると変わります。
- 失敗/生存:このオプションを選択すると、プロットが生存設定から失敗設定に切り替わります。生存設定では、Y軸は1から始まり、時間とともに確率が減少します。Y軸は、生存した(つまり、失敗しなかった)観測値の割合として解釈できます。失敗設定では、Y軸は0から始まり、時間とともに確率が増加します。Y軸は、失敗を経験した(つまり、生存しなかった)観測値の割合として解釈できます。
- 区間の追加/区間の削除:このオプションは、各曲線の周囲に点ごとの95%信頼区間を追加または削除します。
- 同時区間の追加/同時区間の削除:このオプションは、各曲線の周囲に同時95%信頼区間を追加または削除します。
- テストの追加/テストの削除:このオプションは、生存曲線におけるグループ分けレベル間の差を検定するログランク検定およびウィルコクソン検定のp値を追加または削除します。なお、このボタンはグループ分け列が選択されている場合にのみ表示されます。
出力のグラフ領域には、移植片不全までの時間を要約した3行が表示されます。縦のパイプ記号は、少なくとも1人の患者が打ち切りになった時点を示します。リスクのある患者数の表は、ALD、PBC、PSCの患者数がそれぞれ933人、469人、359人で研究が開始されたことを示しています。値は365日ごとに表示されますが、最後の値のセットは4745日(13 × 365)がこの最大日数を超えているため、4720日に表示されます。ダイアログオプションにより、リスクのある患者数の表は、主要疾患グループ列の線と信頼区間の色に一致する色で印刷されます。
「テストを追加」をクリックすると、ログランク検定とウィルコクソン検定のp値(それぞれ0.0004と0.0014)が追加されます。
Subsetting the Data Filter to Gender = 2 instantly updates the グラフビルダー results to those presented in Figure 3.
Figure 3. Analysis of liver transplant data for gender = 2図3. 性別=2の肝移植データの分析
ここで、リスクのある患者数の表によると、ALD、PBC、PSCの患者数はそれぞれ200人、400人、94人です。このサンプルサイズの変更により、移植片不全までの時間は有意とはみなされません(ログランク検定:0.7094、ウィルコクソン検定:0.8926)。
The speed at which グラフビルダー updates is especially notable when subsetting using a continuous variable such as age. See the video below.
参考文献
Collett D. (2023).医学研究における生存データのモデリング、第4版。フロリダ州ボカラトン:CRC Press。