안녕하십니까?
"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 | 생존분석 [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 |