안녕하십니까?
"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분석에서는 어떻게 하는게 좋을까요??
코드를 붙여서 해보니 잘 안 되는 것 같아서요 ㅎㅎ
No. | Subject | Author | Date |
---|---|---|---|
1254 | moonbook package mycsv 관련 오류 문의드립니다. [2] | qwetyui | 2023.11.30 |
1253 | logistic regression 명령어 문의드립니다. [2] | 로렌쵸메디치 | 2023.11.22 |
1252 | spss 파일 불러오기. [1] | 통계린 | 2023.11.22 |
1251 | PSM 매칭변수 숫자 관련 문의입니다. [1] | 오하나 | 2023.11.22 |
1250 | Multiple imputation 관련 질문 및, age adjusted incidence rate 관련 질문을 드립니다. | YSKBH | 2023.11.15 |
1249 | semMediation 설치 재문의드립니다 [3] | 마술샘 | 2023.11.11 |
1248 | 안녕하세요. 정회원 문의 [1] | lionking | 2023.11.06 |
1247 | PSM 관련 질문드립니다. | dr. | 2023.11.05 |
1246 | ggplot [1] | 통계린 | 2023.10.30 |
1245 | 정회원 [1] | 통계린 | 2023.10.29 |
1244 | PSM 관련 질문 드립니다. [2] | 로렌쵸메디치 | 2023.10.25 |
1243 | forest plot 그림이 전체가 다 나오지 않습니다. [1] | joyguni | 2023.10.08 |
1242 | PSM caliper default 값 관련 문의드립니다. [1] | Ooppps | 2023.10.05 |
1241 | reliability test tool,,, [1] | FAdavid | 2023.10.05 |
1240 | 정회원 가입 확인 요청드립니다. [1] | osman | 2023.10.03 |
1239 | 생존분석, cox 회귀 관련 기초적인 질문 입니다. | 만니톨 | 2023.09.27 |
1238 | Twang-mnps에서 매칭된 데이터를 추출하는 방법 좀 알려주세요 | Dorocy | 2023.09.26 |
1237 | 영수증 발급을 부탁드립니다. [1] | 규빈어멈 | 2023.09.26 |
1236 | 정회원 확인 [2] | lionking | 2023.09.24 |
1235 | autoReg 에러 문의드립니다 [1] | 소소 | 2023.09.23 |