このウィッシュリストのリクエストのきっかけは何ですか?
複数の治療群(層)を対象に、日常的に生存分析を実施しています。カプラン・マイヤー分析では、ログランク検定とウィルコクソン検定が提供され、全群間の差異を評価できます。モデルが有意差を示した場合、適切なα補正を適用した全群の対比較による事後分析を実施します。これは、分散分析(ANOVA)に続いてTukeyのHSD検定(Bonferroni補正)を実行するのと類似しています。SASでは、Sidak補正を適用した層間でログランク検定を実行できます。
どのような改善を望みますか?
ペア比較を行うことで、どの治療グループが他の治療グループと異なるかをより細かく判定できるため、生存分析がより強力になります。
SAS の場合: STRATA 'グループ化変数'/ test=logrank adjust=sidak;
出力: (xxxx と yyyy はさまざまなグループ (層) を表します)
| 階層間の平等性のテスト |
| テスト |
カイ二乗 |
DF |
Pr >カイ二乗 |
| ログランク |
259.7179 |
4 |
<.0001 |
| ログランク検定における多重比較の調整 |
| 地層の比較 |
|
p値 |
| グループ |
グループ |
カイ二乗 |
生 |
突然の検査 |
| xxxxx |
ええと |
234.6 |
<.0001 |
<.0001 |
| xxxxx |
ええと |
17.1415 |
<.0001 |
0.0003 |
| xxxxx |
ええと |
3.5276 |
0.0604 |
0.4634 |
| xxxxx |
ええと |
1.2451 |
0.2645 |
0.9537 |
| xxxxx |
ええと |
82.0587 |
<.0001 |
<.0001 |
| xxxxx |
ええと |
155.3 |
<.0001 |
<.0001 |
| xxxxx |
ええと |
169.1 |
<.0001 |
<.0001 |
| xxxxx |
ええと |
4.8642 |
0.0274 |
0.2427 |
| xxxxx |
ええと |
8.5731 |
0.0034 |
0.0336 |
| xxxxx |
ええと |
0.5678 |
0.4511 |
0.9975 |
SAS の個別の分析では、各グループに文字を割り当てて、重要度レベルを示すことができます。
このアイデアはなぜ重要なのでしょうか?
複数の治療群と対照群を用いた実験では、最終的には群間の差異の分析が必要になります。少数の群であれば手作業で行うことができますが、多数の群になると非常に面倒な作業となります。
@DensityOfSpaces様、ご要望ありがとうございます。新しいサバイバルプラットフォームにこの機能を追加する予定がある場合は、改めてお知らせいたします。