안녕하십니까?
"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 | R서버 설치 문의드립니다. [4] | 화모기 | 2015.05.13 |
1253 | 초초 초보자 문의 드립니다^^ [2] | 매직하마 | 2017.12.22 |
1252 | mytable() 함수에 대한 질문입니다. [3] | 브리디온 | 2017.11.19 |
1251 | 묻고답하기 게시판을 운영합니다 [1] | cardiomoon | 2015.04.20 |
1250 | ANOVA, Kruskal-Wallis 검정 및 다중비교 [3] | 소오름 | 2016.03.31 |
1249 | ios | bnmy6581 | 2018.12.28 |
1248 | 예제 데이터 열 제목 편집방법 [2] | 신돌석 | 2016.04.25 |
1247 | rstudio 사이트에 연결할 수 없음 문제 [1] | 지원 | 2019.09.10 |
1246 | 로지스틱 회귀분석에 대해 질문 올립니다. [12] | SOL | 2017.09.15 |
1245 | Cox proportional analysis에서 subgroup분석에 대해 여쭙니다 | 김규 | 2017.09.30 |
1244 | rJava 에러 문제 해결에 대해 문의드립니다 [1] | mdhjcho | 2017.06.29 |
1243 | 3D 형태 그래프에 대해 질문드립니다. | 백지백지 | 2019.04.26 |
1242 | 메타분석의 변수지정 [1] | 야스오 | 2021.05.24 |
1241 | Logistic regression [3] | namuhanayeyo | 2017.03.02 |
1240 | 생존 분석 시 Cox-regression vs. Log-rank [2] | begabung | 2017.05.22 |
1239 | 표본설계시 subset 명령어 문의드립니다. [1] | 네잎클로바 | 2020.05.09 |
1238 | Rstudio 사용법` [3] | 방랑자 | 2015.05.02 |
» | Survdiff로 log-rank test 시 P-value 소수점 [3] | begabung | 2019.02.15 |
1236 | length 오류질문 [1] | 촙춉 | 2020.05.26 |
1235 | 메타분석-출판편향 [1] | 날아라 | 2019.08.26 |