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
1093 Box plot 그리프 설정에 관해 문의합니다. papahong 2022.06.20
1092 cluster data 의 risk difference 에 대해 여쭙고자 합니다. lagrima 2022.06.19
1091 Pair plot을 할 경우 [1] PipeDragoN 2022.06.13
1090 PSM에서 matching 버튼 위치 문의 통계어려워 2022.06.13
1089 PSM (4그룹) 관련 문의드립니다. emksh77 2022.06.08
1088 정회원후원하였습니다. [1] secret rosemary 2022.06.08
1087 PSM 문의 [2] secret sfumatoi 2022.06.07
1086 multipleROC 패키지에서 multivariable의 cutoff 구하기 file 찬쓰 2022.06.07
1085 survival analysis관련 문의 드립니다 [1] cebukim 2022.06.04
1084 autoReg 설치 관련 오류 [5] 닥터강 2022.05.24
1083 survival tree (RPA) 분석 관련 문의 드립니다. [2] 곰곰11 2022.05.23
1082 통계 질문입니다. [1] file 봄봄봄 2022.05.23
1081 데이터 불러오기 문제 [3] strat 2022.04.25
1080 r프로그래밍 문자는 factor로 읽을 수 없나요? [1] file 채채채추ㅐㅊ채야채 2022.04.21
1079 R 코드와 관련하여 문의 드려도 될까요..? [5] file 꾸이뭉 2022.04.19
1078 R 업데이트 후 mytable 함수 에러가 있습니다. [1] 카이토르 2022.04.18
1077 Survival curve 에서 cumulative incidence % 로 표현 가능할까요 [1] nyknicks 2022.04.16
1076 multiple logistic regression test에서 설명변수 선택에 대해 여쭙습니다. [2] file kms9057 2022.04.16
1075 PSM 문의 드립니다 [1] caesar 2022.04.07
1074 Survival regression에서 Survival probability 그래프에서 회색음영이 무엇인지요? [2] file 케로스 2022.03.24