추정결과 그리기

추정계수 그리기

기존의 관행에 따르면 추정결과를 표로 보고하는 것이 보통이다.

하지만 최근에는 한 눈에 알아볼 수 있게 그림을 이용하여 보고하는 경우가 늘어나고 있다. 모형간 비교를 해야할 때에는 그래프로 추정결과를 보여주면 가독성이 개선되는 것도 사실이다. 예컨대, 분위회귀에서는 다양한 분위에 대한 추정결과를 표로 정리하는 것보다는 그래프로 정리하면 훨씬 보기가 좋다.

Stata에서는 margins를 이용해도 되지만 최근에 ssc 더해진 coefplot 을 이용하면 쉽게 그릴 수 있다.

필요한 경우 다음과 같이 coefplot을 설치하자.

ssc install coefplot, replace

다음과 같이 자료를 읽어들이자. 그리고 분석의 편의를 위해 표본을 rep78이 3이상인 경우로 한정하자.

. sysuse auto, clear
(1978 Automobile Data)

. keep if rep78>=3
(10 observations deleted)

아래와 같이 회귀분석을 하자. mpg를 여러가지 변수에 대해 회귀분석을 하자.

. reg mpg headroom i.rep##i.foreign

      Source |       SS           df       MS      Number of obs   =        59
-------------+----------------------------------   F(6, 52)        =      5.33
       Model |  829.833037         6  138.305506   Prob > F        =    0.0002
    Residual |  1350.40425        52  25.9693125   R-squared       =    0.3806
-------------+----------------------------------   Adj R-squared   =    0.3091
       Total |  2180.23729        58  37.5902981   Root MSE        =     5.096

-------------------------------------------------------------------------------
          mpg |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
--------------+----------------------------------------------------------------
     headroom |  -1.120753   .9886418    -1.13   0.262    -3.104608    .8631023
              |
        rep78 |
           4  |  -.3064994   1.973721    -0.16   0.877    -4.267059    3.654061
           5  |   11.91038   3.856209     3.09   0.003      4.17233    19.64843
              |
      foreign |
     Foreign  |   3.710693   3.149593     1.18   0.244    -2.609427    10.03081
              |
rep78#foreign |
   4#Foreign  |   1.675263   3.941705     0.43   0.673    -6.234348    9.584873
   5#Foreign  |  -8.972643   5.129298    -1.75   0.086    -19.26533    1.320046
              |
        _cons |   22.61131    3.33317     6.78   0.000     15.92282    29.29981
-------------------------------------------------------------------------------

이상의 추정결과를 아래와 같이 coefplot을 이용하여 그림을 그릴 수 있다.

. coefplot, xline(0) drop(_cons) omitted baselevels cismooth grid(none) ///
> headings(  ///
>         3.rep78 ="{bf:repair record}" 0. ///
>         0.foreign = "{bf:car type}" ///
>         3.rep78#0.foreign = "{bf:interaction}" ///
>         )

아래와 같은 그림을 얻을 수 있다. coefpot의 옵션은 비교적 다양해서 옵션을 조정하면 다양한 그림을 그릴 수 있다.