안녕하십니까?
"Survival" 패키기로 log-rank test를 시행하는데 아래와 같은 결과가 나왔습니다.
===================================
> #-- Progression Free Survival by Stage
> survdiff(Surv(T_PRG, O_PRG==1)~Stage, data=dset1)
Call:
survdiff(formula = Surv(T_PRG, O_PRG == 1) ~ Stage, data = dset1)
N Observed Expected (O-E)^2/E (O-E)^2/V
Stage=3 50 19 24.4 1.188 1.74
Stage=4 103 72 66.6 0.435 1.74
Chisq= 1.7 on 1 degrees of freedom, p= 0.2
==================================
위에서 P value가 0.2 로 소수점 한자리까지만 표시가 되는데, 혹시 소수점 둘째자리나 셋째자리까지 표시할 수 있는 방법이 있는지 궁금합니다.
감사합니다!
Comment 3
-
마술샘
2019.02.18 00:05
안녕하세요 마술샘입니다다른 부분은 쉽게 자료에 접근이 가능한데 p값은 접근이 어렵네요먼저 x2값에 접근하는 방법입니다 이 값도 소수점 한자리죠...survdiff(Surv(futime, fustat) ~ rx,data=ovarian)[["chisq"]]survdiff(Surv(time, status) ~ pat.karno + strata(inst), data=lung)[["chisq"]]그 다음이 p값의 접근인데... 이거는 직접 구해줘야 합니다명령어는 다음과 같습니다pchisq(여기에 x2입력, df=여기에 df입력, lower.tail=FALSE)df는 통계결과에 나옵니다...이를 조합하면pchisq(survdiff(Surv(futime, fustat) ~ rx,data=ovarian)[["chisq"]], df=1, lower.tail=FALSE)pchisq(survdiff(Surv(time, status) ~ pat.karno + strata(inst), data=lung)[["chisq"]], df=7, lower.tail=FALSE)그리고 소수점 3자리까지만 표기하고 싶으시면round(pchisq(survdiff(Surv(futime, fustat) ~ rx,data=ovarian)[["chisq"]], df=1, lower.tail=FALSE),3)round(pchisq(survdiff(Surv(time, status) ~ pat.karno + strata(inst), data=lung)[["chisq"]], df=7, lower.tail=FALSE),3)해당 예시는 기본 생존분석 예시입니다도움이 되셨으면 합니다... -
begabung
2019.02.21 21:09
복잡하지만... 상세하게 알려주셔서 성공했습니다! 감사합니다.
-
Moncher
2019.02.25 10:45
저도 도움 많이 받았습니다!
혹시 cox분석에서는 어떻게 하는게 좋을까요??
코드를 붙여서 해보니 잘 안 되는 것 같아서요 ㅎㅎ