애니메이션 2

C. Huber의 Stata blog의 애니메이션 예제 파일 중 2번째 프로젝트

lowes 추정 과정을 보여줄 수 있는 애니메이션입니다.

sysuse auto, clear

local width = 1500

forvalues UB = 3500(500)15000 {

  local LB = `UB' - `width'

  twoway (scatter mpg price , color(%30) mlwidth(none))                          ///

    (lowess mpg price  if price < (`UB'-`width'/2), lcolor(pink%80) )    ///

    (lfit mpg price  if price > `LB' & price <`UB',      ///

         lwidth(medium) lcolor())                                      ///

    , xline(`UB' `LB', lwidth(medium) lcolor(red%50) lpattern(dash))  ///

    legend(on order(1 2 3) cols(3) pos(6)) scheme(538)

  graph export lowess_`UB'.png, as(png) width(1280) height(720) replace

}