confidence ellipse 그리기

confidence ellipse 그리기

1변수인 경우에는 신뢰구간이 되지만 2변수인 경우에는 신뢰타원이 된다.

신뢰타원을 쉽게 그리기는 좀 어려운데 Stata에서는 ellip 을 이용하여 그릴 수 있다.

필요한 경우 다음과 같이 ellip를 설치하자.

ssc install ellip, replace

다음과 같이 자료를 읽어들이자.

. sysuse auto, clear
(1978 Automobile Data)

먼저 ellip 을 이용하여 confidence ellipse를 그려보자.

. ellip mpg weight, ///
> means color(green%50) lp(solid) ///
> by(foreign , note("") ) ///
> plot(scatter mpg weight, color(pink%30) ms(O) mlwidth(none)) ///
> legend(row(1))

그러면 아래와 같은 confidence ellipse를 그릴 수 있다.

한편, 추정치의 경우에도 유사하게 그릴 수 있다.

. reg price mpg weight

      Source |       SS           df       MS      Number of obs   =  
>       74
-------------+----------------------------------   F(2, 71)        =  
>    14.74
       Model |   186321280         2  93160639.9   Prob > F        =  
>   0.0000
    Residual |   448744116        71  6320339.67   R-squared       =  
>   0.2934
-------------+----------------------------------   Adj R-squared   =  
>   0.2735
       Total |   635065396        73  8699525.97   Root MSE        =  
>     2514

----------------------------------------------------------------------
> --------
       price |      Coef.   Std. Err.      t    P>|t|     [95% Conf. I
> nterval]
-------------+--------------------------------------------------------
> --------
         mpg |  -49.51222   86.15604    -0.57   0.567    -221.3025    
>  122.278
      weight |   1.746559   .6413538     2.72   0.008      .467736    
> 3.025382
       _cons |   1946.069    3597.05     0.54   0.590    -5226.245    
> 9118.382
----------------------------------------------------------------------
> --------

. ellip mpg weight , coefs color(blue%50) lp(dash)