Skip to menu

R과 SPSS에서 시행한 logistic regression값은 거의 일치하는데, Web-R에서 돌린 것은 일치하지가 않습니다.

namuhanayeyo 2016.04.27 17:23 Views : 681

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.11400.89300.130.89841.120.206.87
FEV1-1.31180.7604-1.730.08450.270.051.07
Home_O2Yes1.15930.78891.470.14173.190.6715.77


p-value threshold entering multiple regression가 달라서 그런가 해서 0.1로 바꿔봤는데, 변화가 없습니다.

논문 몇개를 Web-R로 logistic regression을 했는데 맨붕이 되었습니다.

도와주셨으면 합니다. ㅜㅜ

(파일 첨부했습니다.)


감사합니다.


I want to
List
No. Subject Author Date
1174 Web R studio 초기화 후 문의드립니다. [2] secret 화모기 2016.02.23
1173 [정보] 구조방정식 모형 분석 관련 참고자료 소개 [2] 보석 2016.02.24
1172 그래프가 이상하게 나와서 어떻게 하면 좋을까요? [3] file osr33260216 2016.02.25
1171 데이터 업로드가 되지 않습니다.. [3] jenni 2016.02.28
1170 어제 강의 너무 잘 들었습니다 [5] 마술샘 2016.02.29
1169 lavaan의 공식 구글 토론실입니다 마술샘 2016.02.29
1168 [오류]web-r 구조방정식 [2] 마술샘 2016.02.29
1167 저 이 홈피 체험 1주일 지난것같은데... [3] osr33260216 2016.03.03
1166 이럴 땐 어떤 통계 기법을 쓰나요? (기초) [1] file 송이버섯 2016.03.03
1165 남북한별 인구를 나타내는 막대 그래프를 만들려는데 왜 바가 않나타나죠? [1] file osr33260216 2016.03.04
1164 제가 금시세 표를 만들었는데 왜 업로드 파일은 이상하게 나왔나요? [1] file osr33260216 2016.03.04
1163 Standard error 질문 드립니다. [7] file R까기 2016.03.05
1162 object 'res' not found [3] 민자 2016.03.06
1161 제가 선생님 말씀대로 금시세표를 csv쉼표분리로 변경했는데 이상하게되었습니다. [1] file osr33260216 2016.03.07
1160 인공지능, 인지 컴퓨팅, 머신 러닝, 딥 러닝 관련 IBM 자료 [1] file Ben 2016.03.09
1159 [RStudio Server] 사용 연장 문의 [2] KB 2016.03.09
1158 descriptive 표만들기 질문입니다. [4] 난둥 2016.03.10
1157 그래프 작성 문의드립니다. [4] file 도치맘 2016.03.10
1156 표만들기 관련 질문입니다. [1] tympa 2016.03.10
1155 머신러닝 레볼루션 행사안내 [1] mike 2016.03.10
Write Tag