다변수 로지스틱 분석에서 범주형 설명 변수의 빈도표는 unitab 명령어를 써보자. OR과 95% CI는 logistic 명령어와 같으므로 비교해보자. (unitab 명령어는 사용자 작성 명령어므로 찾아서 설치해야 동작한다.)

. webuse lowbirth
 

. unitab low age lwt, c(ptd smoke ht ui)
------------------------------------------------------------------------------
     low    |    low=1(%)      Total(%)     OR   [95% Conf. Interval]  p-value
------------+-----------------------------------------------------------------
     age    |      56(50)       112(100)   1.000    0.917     1.090     1.000
------------+-----------------------------------------------------------------
     lwt    |      56(50)       112(100)   0.989    0.976     1.002     0.093
------------+-----------------------------------------------------------------
     ptd  0 |      38(44)        87(78)     .        .         .        0.013
          1 |      18(72)        25(22)    3.316    1.256     8.750
------------+-----------------------------------------------------------------
   smoke  0 |      26(39)        66(59)     .        .         .        0.007
          1 |      30(65)        46(41)    2.885    1.319     6.307
------------+-----------------------------------------------------------------
      ht  0 |      49(48)       102(91)     .        .         .        0.185
          1 |       7(70)        10( 9)    2.524    0.618    10.308
------------+-----------------------------------------------------------------
      ui  0 |      42(46)        92(82)     .        .         .        0.048
          1 |      14(70)        20(18)    2.778    0.981     7.864
------------------------------------------------------------------------------

unitab 명령어 바로 뒤는 결과 변수이고 연속형 변수는 결과 변수 뒤에, 범주형 변수는 옵션으로 c( )에 넣으면 끝. 연속형 변수의 low=1 (%)는 low 변수가 1로 코딩된 빈도와 백분율이다.

'통계 연습' 카테고리의 다른 글

[Stata] batplot 명령어  (0) 2011.12.26
[Stata] oddsrisk 명령어  (0) 2011.10.20
[Stata] 표준화비 신뢰구간 계산  (0) 2011.10.14
[Stata] 위험 차(risk difference) 계산  (2) 2011.10.14
[Stata] optmatch2 명령어  (0) 2011.10.14
Posted by cyberdoc
: