바쁘신데 또 질문을 올리게 되었습니다.
제 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 가 뜹니다.
어떠한 함수를 사용해야 제가 원하는 것을 그릴 수 있을 지요.
Comment 0
No. | Subject | Author | Date |
---|---|---|---|
1094 | 로지스틱 회귀분석 관련 질문 [2] | alse | 2016.06.18 |
1093 | 나무 분석 결과 문제 [3] | 매일이5월 | 2016.06.19 |
1092 | linear regression에서 회귀모형 적합성 여부 [3] | 매일이5월 | 2016.06.19 |
1091 | 학생강의project2 | cardiomoon | 2016.06.21 |
1090 | ROC 분석중 문제가 생겨서 질문 드립니다 [2] | candid00 | 2016.06.26 |
1089 | moonBook2 install중 error message [2] | hawk1227 | 2016.06.27 |
1088 | 표 만들기에서 에러가 생깁니다 [2] | jaljary | 2016.06.27 |
1087 | 선그래프 질문드립니다. [2] | 도치맘 | 2016.07.02 |
1086 | plot.or그래프 그리는 방법 [4] | 두치님 | 2016.07.04 |
1085 | 교수님 결측치 관련하여 문의드립니다. [1] | tunneled_vision | 2016.07.07 |
1084 | 질문드립니다. [1] | 죵죵죵 | 2016.07.10 |
1083 | ROC에서 여러설명변수 갖는 model간 비교 [5] | alse | 2016.07.11 |
1082 | excel 변환 부탁드립니다. [3] | 죵죵죵 | 2016.07.12 |
1081 | 다중 선형 회귀 분석에서 "relative importance of predictor variables"에 관해 질문드립니다. [2] | Ikarus86 | 2016.07.13 |
1080 | 시계열 회귀와 회귀 [1] | 반디12 | 2016.07.19 |
1079 | 데이타 선택 [4] | 김현정 | 2016.07.19 |
1078 | 로그변환 후 회귀분석 시 결과물 [1] | 반디12 | 2016.07.19 |
1077 | 데이터 분석 후 예상과 결과가 다를 때, [1] | 반디12 | 2016.07.19 |
1076 | 샘플사이즈 문제 [3] | 반디12 | 2016.07.20 |
1075 | 표만들기에서 이런 에러가 생기는데 뭐가 문제인지요? [2] | 지훈이아빠 | 2016.07.22 |