메뉴 건너뛰기

웹에서 하는 R 통계

안녕하세요 교수님 

항상 많은 도움 주셔서 감사하고

 

특히 최근에 나온 autoReg package를 통해 손쉽게 publication ready 상태로 출력하게 해주셔서 정말 

정말 많은 도움 받았습니다. 또한  어려워서 R로 구현하지 못하고 있었던 Grey test도 가능하게 해주셔서 논문 작성에 많은 도움이 될꺼같습니다. 

 

 

 

제가 최근에 관심 있는건 예후 예측 점수 만들기 인데요. 사실 점수 만드는 거 자체는 cox를 이용하여 하였는데요, 

그 이후 이점수로 survival을 예측할때 필요한 값을 어떻게 추출해야하는지 모르겠습니다. 

 

공부해보니... 

The probability of death at time ‘‘t’’ can be estimated by the equation:

(여기서 CLIF-C ACLF 장소에 다른 예측점수를 넣으려고 하는데요) 

 

인데요, 여기에 적용하기 위해서는 

 

28day, 90day, 180day, 365day각각의 cumulated baseline hazard (CI(t) and the score coefficient(베타(t)가 필요하는데요, 

 

이 CI(t)와 베타 (t)를 어떻게 구해야할까요?

 

제가 구글링해서 떠듬떠듬 만든 코드는 다음과 같습니다. 

 

fit1<-coxph(Surv(SUD,Final_Pt)~new_score, data=abc)

fit1

 

#cumulated baseline hazard estimate

sv0<-survfit(fit1,newdata=data.frame(new_score=0))  <- 이게 맞을까요...?;;

summary(sv0)

 

#28-days

s28<-min(summary(sv0)$surv[summary(sv0)$time<=28]) <- 이게 맞을까요...?;;

s28

#baseline hazard maybe for 28-days is 

CL28<-1-s28

CL28

 
## beta for 28-days
fit28<-coxph(Surv(SUD28,Final_Pt28)~new_score, data=abc)
B28<-fit28$coefficients
 
 
CL28
B28
...
 
개인적 의견으로는 베타를 구하는 방식은 맞는거같은데요,,, cumulated baseline hazard estimate을 구하는 방법이 정말 저게 맞는지 궁금합니다..;
basehaz()이러한 기능도 흔히 쓴다고 하는데 그럼 각각 28day, 90day, 180day, 365day날짜에 맞는 cumulated baseline hazard estimate는 그결과에서 무엇을 봐야하는걸까요...?
 
항상 많은 도움주셔서 감사하고 
이번에도 미리 감사의 말씀 올립니다. 
 
 
 
 
번호 제목 글쓴이 날짜 조회 수
1083 통계 질문입니다. [1] file 봄봄봄 2022.05.23 52
1082 데이터 불러오기 문제 [3] strat 2022.04.25 93
1081 r프로그래밍 문자는 factor로 읽을 수 없나요? [1] file 채채채추ㅐㅊ채야채 2022.04.21 79
1080 R 코드와 관련하여 문의 드려도 될까요..? [5] file 꾸이뭉 2022.04.19 105
1079 R 업데이트 후 mytable 함수 에러가 있습니다. [1] 카이토르 2022.04.18 191
1078 Survival curve 에서 cumulative incidence % 로 표현 가능할까요 [1] nyknicks 2022.04.16 107
1077 multiple logistic regression test에서 설명변수 선택에 대해 여쭙습니다. [2] file kms9057 2022.04.16 86
1076 PSM 문의 드립니다 [1] caesar 2022.04.07 122
1075 Survival regression에서 Survival probability 그래프에서 회색음영이 무엇인지요? [2] file 케로스 2022.03.24 87
1074 Bonferroni's adjustment 관련 문의드립니다. [1] kjptrn 2022.03.23 100
1073 파일 업로드 문의드립니다. [1] ulysses11 2022.03.07 153
1072 x 축 변수가 어려개인 그래프를 그릴 수 있을까요? [1] 나현둥 2022.02.20 100
1071 독립변수가 이상한 (?) 경우의 로지스틱 회귀분석 [4] papahong 2022.02.20 107
1070 survival tree 분석을 하였는데 node가 한 개만 나옵니다 [1] 뒤통수조심 2022.02.10 186
» 예후예측 점수 만든후 그것으로 survival 예측하는 데에 문제 겪고있습니다. [1] 권혜미 2022.02.08 100
1068 기초적인 통계에서 평균, 중간값, 최빈값에 대한 질문입니다. 우엉우왕 2022.02.04 197
1067 메타분석 plot 출력 문제와 관련해 답변해 주신 내용에 따라 다시 글을 올립니다. [2] secret endostat 2022.01.27 5
1066 메타분석 후 plot을 다운로드했을 때 일부 plot이 나오지 않습니다. [1] endostat 2022.01.27 106
1065 에러가 났는데 어떻게 해야 하나요? [1] 세렌 2022.01.21 1477
1064 multipleROC package 질문 [2] 구현아빠 2022.01.13 157