메뉴 건너뛰기

웹에서 하는 R 통계



 안녕하세요. 


 여기에는 처음 인사드립니다. 


 올해 들어 R을 접하게 되어 R로 데이터 분석을 하고 있는 신참 연구자입니다. 


 논문 작성을 위해 데이터 분석을 진행 중에 다음과 같은 결과를 얻었습니다. 


glm(formula = Ep8 ~ factor(IFG) + factor(sex) + factor(HTN) + 

    age + LVMI + bmi + heartrate, family = binomial, data = result)

 

Deviance Residuals: 

    Min       1Q   Median       3Q      Max  

-2.7671  -0.6920  -0.3915   0.6164   3.3010  

 

Coefficients:

               Estimate Std. Error z value Pr(>|z|)    

(Intercept)  -15.677580   0.236574 -66.269  < 2e-16 ***

factor(IFG)1   0.128478   0.032657   3.934 8.35e-05 ***

factor(sex)1   0.450754   0.040347  11.172  < 2e-16 ***

factor(HTN)1   0.562605   0.041035  13.710  < 2e-16 ***

age            0.143539   0.002363  60.736  < 2e-16 ***

LVMI           0.017034   0.001117  15.255  < 2e-16 ***

bmi            0.164362   0.005370  30.609  < 2e-16 ***

heartrate      0.046502   0.001797  25.880  < 2e-16 ***

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

 

(Dispersion parameter for binomial family taken to be 1)

 

    Null deviance: 36957  on 30942  degrees of freedom

Residual deviance: 27425  on 30935  degrees of freedom

  (87 observations deleted due to missingness)

AIC: 27441

 

Number of Fisher Scoring iterations: 5

 

> extractOR(out6)

               OR  lcl  ucl     p

(Intercept)  0.00 0.00 0.00 0e+00

factor(IFG)1 1.14 1.07 1.21 1e-04

factor(sex)1 1.57 1.45 1.70 0e+00

factor(HTN)1 1.76 1.62 1.90 0e+00

age          1.15 1.15 1.16 0e+00

LVMI         1.02 1.01 1.02 0e+00

bmi          1.18 1.17 1.19 0e+00

heartrate    1.05 1.04 1.05 0e+00



 여기서 0e+00 라는 표현이 해석이 잘 되지 않아서 여쭤봅니다. 동일한 내용을 ORplot 으로 분석하면 *** 로 모두 유의한 P 값으로 나오게 됩니다. 분석 내용으로 봤을 때는 P값이 0.001 이하일 것 같은데요.... 


 제가 잘 몰라서 여쭤보니 혹시 너무 기초적인 내용을 모르더라도 양해 부탁드립니다. 


 

그리고 moonBook 패키지는 아주 잘쓰고 있습니다. R을 처음 공부할 때 문건웅 선생님 패키지를 알게 되어 많은 어려움이 해결되었습니다. 


 정말 감사드립니다.