Skip to menu

cph plot

mint 2016.12.14 12:59 Views : 606

바쁘신데 또 질문을 올리게 되었습니다.


제 data 에서 spss cox 로 돌리면 HPA_selection_IgA 라는 변수가 Acute rejection ever 라는 종속 변수에 대하여 통계적으로 유의한 값이 나오는데


다른 변수들을 같이 넣고 cph plot 을 그릴 수 있는지요


require(survival)

require(rms)

Rejection<-with(IgAN, Surv(TPL_to_everAR,Rejection_ever)

fit=coxph(Rejection~rcs(HPA_selection_IgA,4)+Age+Sex+Donor_age+Donor_sex+DSA+Desesntization,data=IgAN) 

이렇게 해보았는데 결과가 나오지 않습니다.


Cox proportional hazard ratio 그림에서 Acute rejection ever 를 outcome 으로 TPL_to_everAR 을 time 으로 해서

HPA_selection_IgA 를 주요 독립 변수로 하고

다른 변수들, Age, Sex, Donor age, Donor sex, DSA, Desentization 등의 다변수 요인으로 보정해서

그래프를 그리려면 어떠한 함수를 사용해야 할 지 고견을 부탁드립니다.


다른 방법으로 

surv<-Surv(IgAN$TPL_to_everAR,IgAN$Rejection_ever)

 

fit<-coxph(surv~pspline(IgAN$HPA_selection_IgA)+IgAN$Age+IgAN$Sex+IgAN$Donor_age+IgAN$Donor_sex+IgAN$TPL_ever+IgAN$DSA+IgAN$Desensitization,data=IgAN)

 

predicted<-predict(fit,type="terms",se.fit=TRUE,terms=1)

 

plot(0,xlab="HPA_selection_IgA",ylab="Hazard Ratio",main="Rejection_ever",type="n",xlim=c(0,100),ylim=c(0,4))


상기처럼 하고는 옆에 표가 x 축 HPA_selection_IgA 로 0 부터 100 까지 뜨고, y 도 Hazard ratio 표시되어 0 부터 4 까지 표시는 되는데

그링믄 안 그려지고


lines(sm.spline(IgAN$APRIL_selection,exp(predicted$fit)),col="black",lwd=0.8) 라고 했는데


Error in tapply(seq(along = y), match(x, ux), function(x, y, w) c(mean(y[x]),  : 

  arguments must have same length 라는 Error 가 뜹니다.

 

어떠한 함수를 사용해야 제가 원하는 것을 그릴 수 있을 지요.

No. Subject Author Date
1213 교수님 지난 번 문의 드렸던 PSM 재문의 드립니다. [3] secret 긍정적인마음가짐 2023.06.27
1212 안녕하세요 ^^ odds ratio 와 관련된 문의드리고자 합니다 ^^ [2] 성열 2023.06.26
1211 mytable() 함수에서 na처리와 pvalue 계산 문의드립니다. [1] hjk34 2023.06.23
1210 교수님 안녕하세요 ^^ web-R table 작성관련 문의드립니다. [2] file mkw7654 2023.06.21
1209 Methods for continuous variables과 P value corrected for multiple comparisons 방법 궁금합니다. [1] 호호셀 2023.06.19
1208 ordinal logistic regression analysis 문의 [1] 버블 2023.06.18
1207 e-mail 변경 문의 [1] secret kimdo 2023.06.16
1206 안녕하세요 number at risk 관련 문의 드립니다 [1] secret cebukim 2023.06.13
1205 Propensity matching score 에 대해 문의드립니다. [1] saga 2023.06.08
1204 콕스비례위험모델의 적용에 관련해서 질문 드립니다. [2] secret kng1032 2023.06.07
1203 기초통계 계산 [2] secret leesav 2023.06.07
1202 정회원 후원하였습니다. [2] secret toast 2023.06.04
1201 logistic regression analysis 가 에러가 뜹니다. [1] secret Ikarus86 2023.06.02
1200 PMS caliper 에 대해 여쭤봅니다. [2] secret pjhaaa 2023.05.30
1199 ggGam에 대하여 [2] leesav 2023.05.22
1198 R 에서 PFS 의 confidence interval 을 구하고 싶을 때 어떻게 하면 될까요? [1] ddeeny 2023.05.21
1197 메타분석 관련 문의드립니다 secret 화이팅 2023.05.12
1196 교수님, 새로운 앱 개발 문의드립니다! [4] swpapa 2023.05.08
1195 교수님, 안녕하세요? 전용 앱 관련 문의드립니다! [1] swpapa 2023.05.07
1194 PSM 시 범주형 변수가 총 10가지 일 때 어떻게 돌려야 하는지요 [1] 긍정적인마음가짐 2023.05.06