* 교수님, 아래와 같이 model을 만들어 보았을때, 207개의 AchRAbtiter의 자료 중에 7개의 자료가 missing value가 있어서 variable이 0과1로 두개로 나누어지는 거 같습니다. 그래서 web-R의 결과와 똑같이 R console에서 구현이 안되는데요...
coxph는 na.rm=T와 같은 추가할 수 있는 항목이 없는지요?? Web-R의 결과를 R에서 구현해보고자 합니다.
* web-R에서는 다름과 같이 표현되었습니다.
For stepwise backward elimination, data containing NA value were removed. Total nrow: 207-> After na.omit : nrow: 200 Total 7 row(s) were removed.
> cox2=coxph(Surv(MMfreeSurvival,MM==1)~AchRAbTiter_6.7_or_more,method="breslow")
> summary(cox2)
Call: coxph(formula = Surv(MMfreeSurvival, MM == 1) ~ AchRAbTiter_6.7_or_more, method = "breslow")
n= 207, number of events= 134
coef exp(coef) se(coef) z Pr(>|z|)
AchRAbTiter_6.7_or_more0 0.66269 1.94001 0.52583 1.260 0.208
AchRAbTiter_6.7_or_more1 0.09358 1.09810 0.51150 0.183 0.855
exp(coef) exp(-coef) lower .95 upper .95
AchRAbTiter_6.7_or_more0 1.940 0.5155 0.6922 5.437
AchRAbTiter_6.7_or_more1 1.098 0.9107 0.4029 2.993
Concordance= 0.555 (se = 0.02 )
Rsquare= 0.039 (max possible= 0.998 )
Likelihood ratio test= 8.28 on 2 df, p=0.01596
Wald test = 9.01 on 2 df, p=0.01105
Score (logrank) test = 9.26 on 2 df, p=0.009772
Comment 2
-
PipeDragoN
2017.01.18 16:51
-
cardiomoon
2017.01.18 23:06
혼자 해결하셨다니 다행입니다. 축하드립니다.^^
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 |
아..해결하였습니다. read.xlsx로 읽으니 공란이 NA처리가 안되어서 공란을 NA로 바꾸어준후 read.csv로 읽으니 NA처리가 되어 R console에서도 web-R과 같이 결과가 구현이 되었습니다.