이 위시리스트 요청은 무엇에서 영감을 받았나요?
우리는 여러 치료군(층)에 대한 생존 분석을 일상적으로 수행합니다. Kaplan-Meier 분석은 모든 그룹 간의 차이를 평가할 수 있는 로그 순위 검정 및 Wilcoxon 검정 통계량을 제공합니다. 모델에서 유의미한 결과가 나타나면, 적절한 알파 보정을 적용하여 모든 그룹 간의 쌍별 비교에 대한 사후 분석을 수행하고자 합니다. 이는 Bonferroni 보정을 적용한 Tukey의 HSD 검정을 수행하는 ANOVA 분석과 유사합니다. SAS에서는 Sidak 보정을 적용하여 층별 로그 순위 검정을 실행할 수 있습니다.
어떤 개선점을 보고 싶으신가요?
쌍별 비교를 수행하면 생존 분석의 정확도가 향상되고, 어떤 치료 그룹이 다른 치료 그룹과 다른지 더욱 세밀하게 파악할 수 있습니다.
SAS에서: STRATA '그룹화 변수'/ test=logrank adjust=sidak;
출력: (xxxx와 yyyy는 다양한 그룹(계층)을 나타냅니다.)
| 계층에 따른 평등성 검증 |
| 시험 |
카이제곱 |
DF |
Pr >카이제곱 |
| 로그랭크 |
259.7179 |
4 |
<.0001 |
| 로그랭크 검정에 대한 다중 비교 조정 |
| 계층 비교 |
|
p값 |
| 그룹 |
그룹 |
카이제곱 |
날것의 |
갑작스러운 검사 |
| xxxxx |
yyyyyyy |
234.6 |
<.0001 |
<.0001 |
| xxxxx |
yyyyyyy |
17.1415 |
<.0001 |
0.0003 |
| xxxxx |
yyyyyyy |
3.5276 |
0.0604 |
0.4634 |
| xxxxx |
yyyyyyy |
1.2451 |
0.2645 |
0.9537 |
| xxxxx |
yyyyyyy |
82.0587 |
<.0001 |
<.0001 |
| xxxxx |
yyyyyyy |
155.3 |
<.0001 |
<.0001 |
| xxxxx |
yyyyyyy |
169.1 |
<.0001 |
<.0001 |
| xxxxx |
yyyyyyy |
4.8642 |
0.0274 |
0.2427 |
| xxxxx |
yyyyyyy |
8.5731 |
0.0034 |
0.0336 |
| xxxxx |
yyyyyyy |
0.5678 |
0.4511 |
0.9975 |
SAS에서 별도의 분석을 통해 각 그룹에 문자를 할당하여 유의 수준을 나타낼 수 있습니다.
이 아이디어가 왜 중요한가요?
다수의 처리군과 대조군이 있는 모든 실험은 궁극적으로 군 간의 차이를 분석해야 합니다. 군의 수가 적을 때는 수작업으로 분석할 수 있지만, 군이 많을 때는 매우 번거로운 작업이 됩니다.
@공간밀도 님, 요청해 주셔서 감사합니다. 새로운 서바이벌 플랫폼에 해당 기능을 추가하는 것을 고려하게 되면 알려드리겠습니다.
