cox regression을 통해서 다변량 분석을 했을 때, 결과는 유의하게 나왔습니다. (p=0.047)
1분위수 (ref) 대비 2분위수와 (HR=4.145), 1분위수 (ref)대비 3분위수는(HR=4.392) 각각 유의하지만, 큰 차이는 나지 않습니다.
그러나 같은변수로 kaplan-meier로 분석했을때는 유의하게 나오지 않았습니다. (p=0.126)
제가 알기론 보통 결과가 일치한다는데, 혹시 이렇게 나온게 2분위수와 3분위수의 결과차이가 크지 않아서 그런건가요?
그렇다면 이렇게 나와도 의미있는지 여쭤보고싶습니다.
Comment 2
-
cardiomoon
2018.05.01 23:20
-
아빠가들려주는통계김지형
2018.05.31 03:05
"kaplan-meier로 분석했을때" 이라는 말은 옳지 않습니다.
아주 쉽게 말하면, kaplan-meier 는 생존 곡선을 그리는 방법이라고 이해할 수도 있습니다.
즉, Cox regression을 시행하고 나오는 생존 곡선도 KM plot입니다.
흔히들 연구자들이. kaplan-meier로 분석했을때 라고 말하는 것은 아마도, 두 집단, 혹은 세 집단을 비교한
log rank test를 말하는 것이라고 생각됩니다. (log-rank 말고도 몇 가지 다른 것이 있지만...대표적으로)
이것은 단변수 분석에 해당하죠.
그래서 Cox했을 때 단변수 분석을 하면 이것은 log- rank test와 비슷한 결과가 나옵니다.
비유하자면, 카이제곱 검정과 단변수 로지스틱 회귀분석과 비슷한 결과가 나오는 것과 유사합니다.
그런데, 보통 cox 하게 되면 다변수 분석을 하게 되지요. 그래서, 단변수 분석한 것과는 차이가 나는 것은 오히려 당연합니다.
위의 경우를 요약하면 아마도 단변수 분석하니 유의하지 않고, 다변수 분석을 하니 0.047 정도 나왔다. 이렇게 요약될 것같습니다.
흔히 있는 경우라고 생각됩니다.
No. | Subject | Author | Date |
---|---|---|---|
1254 | 생존분석 [1] | toby | 2015.04.29 |
1253 | RStudio v0.99 Preview: Code Diagnostics [1] | fermat39 | 2015.04.30 |
1252 | Rstudio 사용법` [3] | 방랑자 | 2015.05.02 |
1251 | need finite 'xlim' values [3] | anthleemd | 2015.05.02 |
1250 | 샤이니 앱에서 SQL 명령어 사용 방법 [2] | 김브링이 | 2015.05.06 |
1249 | 표만들기 관련 [2] | 현승협 | 2015.05.07 |
1248 | MS word문서 저장관련 [1] | 아자 | 2015.05.09 |
1247 | R서버 설치 문의드립니다. [4] | 화모기 | 2015.05.13 |
1246 | 웹에서 하는 R통계 - R코드 [2] | 김영진 | 2015.06.07 |
1245 | Plot inset 만들기 | cardiomoon | 2015.06.09 |
1244 | 오류 질문 - 표작성시 [1] | 미소천사 | 2015.06.15 |
1243 | 오류 질문 - PDF 변환 안되는 file입니다. [1] | 미소천사 | 2015.06.15 |
1242 | print ztable [1] | 김영진 | 2015.06.17 |
1241 | 표 작성시 질문 [1] | simsim | 2015.06.17 |
1240 | 나무 분석 질문입니다. | 미소천사 | 2015.06.19 |
1239 | Interactive ggplot2 [2] | cardiomoon | 2015.06.22 |
1238 | Rstudio server 서버 오류 관련해서 문의드립니다. [2] | 노상래 | 2015.07.22 |
1237 | Rstudio server 무한 로딩 관련해서 문의 드립니다. [6] | 노상래 | 2015.07.24 |
1236 | http://www.r-statistics.com/ 번역소식 [1] | fermat39 | 2015.07.25 |
1235 | 생존분석 변수선택 [3] | nash | 2015.08.11 |
질문 내용이 정확이 이해되지 않습니다. 데이터와 R코드를 올려주시면 한번 살펴보겠습니다.