Web-R과 R에서 logistic regression을 한 값이 다르게 나옵니다.
R과 SPSS에서 돌린 거의 같습니다.
1. R로 돌린 값입니다.
> result <-glm(AE_90D~age+sex+IHD+HF+BA2+FEV1+Home_O2+CRP, family=binomial, data=AE)
> reduced.model2=step(result)
Start: AIC=85.32
AE_90D ~ age + sex + IHD + HF + BA2 + FEV1 + Home_O2 + CRP
Df Deviance AIC
- IHD 1 67.378 83.378
- CRP 1 67.678 83.678
- sex 1 68.143 84.143
- age 1 68.429 84.429
<none> 67.316 85.316
- Home_O2 1 69.693 85.693
- FEV1 1 70.795 86.795
- BA2 1 70.820 86.820
- HF 1 70.858 86.858
Step: AIC=83.38
AE_90D ~ age + sex + HF + BA2 + FEV1 + Home_O2 + CRP
Df Deviance AIC
- CRP 1 67.707 81.707
- sex 1 68.238 82.238
- age 1 68.541 82.541
<none> 67.378 83.378
- Home_O2 1 69.731 83.731
- BA2 1 70.833 84.833
- FEV1 1 70.849 84.849
- HF 1 71.261 85.261
Step: AIC=81.71
AE_90D ~ age + sex + HF + BA2 + FEV1 + Home_O2
Df Deviance AIC
- sex 1 68.661 80.661
- age 1 68.869 80.869
<none> 67.707 81.707
- Home_O2 1 70.545 82.545
- FEV1 1 70.897 82.897
- BA2 1 71.099 83.099
- HF 1 71.877 83.877
Step: AIC=80.66
AE_90D ~ age + HF + BA2 + FEV1 + Home_O2
Df Deviance AIC
- age 1 69.326 79.326
<none> 68.661 80.661
- FEV1 1 71.177 81.177
- BA2 1 71.417 81.417
- Home_O2 1 72.511 82.511
- HF 1 73.026 83.026
Step: AIC=79.33
AE_90D ~ HF + BA2 + FEV1 + Home_O2
Df Deviance AIC
<none> 69.326 79.326
- FEV1 1 72.081 80.081
- BA2 1 72.483 80.483
- Home_O2 1 73.297 81.297
- HF 1 73.298 81.298
> extractOR(reduced.model2)
OR lcl ucl p
(Intercept) 0.48 0.06 3.66 0.4750
HFYes 7.04 1.03 49.27 0.0415
BA2Yes 3.30 0.89 14.08 0.0837
FEV1 0.29 0.05 1.23 0.1195
Home_O2Yes 5.69 1.03 35.40 0.0497
그런데 Web-R로 돌린 값은 아래와 같습니다.
(Intercept) | 0.1140 | 0.8930 | 0.13 | 0.8984 | 1.12 | 0.20 | 6.87 |
p-value threshold entering multiple regression가 달라서 그런가 해서 0.1로 바꿔봤는데, 변화가 없습니다.
논문 몇개를 Web-R로 logistic regression을 했는데 맨붕이 되었습니다.
도와주셨으면 합니다. ㅜㅜ
(파일 첨부했습니다.)
감사합니다.
Comment 2
-
cardiomoon
2016.04.27 20:23
-
namuhanayeyo
2016.04.27 22:22
감사합니다.
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 |
단변량분석에서 p 값을 보면 FEV1과 Home_O2Yes를 제외한 값들이 모두 0.2를 훨씬 넘습니다. 0.2 이하인 값만 multiple logistic regression을 해서 나온 결과 입니다. 만일 모든 변수를 다중로지스틱회귀분석을 하시려면 p vlaue threshold를 1로 하시면 같은 결과를 얻으실 수 있습니다