penalized-spline regression fitting

Penalized-spline line 추정

Rupper et al.(2003)에 기반하여 y를 x에 penalized-spline regression을 하고 

이를 이용하여 함수를 플롯팅을 하자.

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

ssc install pspline, replace

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

. use http://fmwww.bc.edu/repec/bocode/l/lidar.dta, clear
(Light detection and ranging (LIDAR) data)

그리고 다음과 같이 pspline 를 이용하여

. pspline logratio range, g(fitted) nograph  
(pilot goodness-of-fit chi2(35) = 411.98; p = 0.0000)
(using penalized model ...)

. tw (sc logratio range, color(green%30) mlwidth(none) ) (line fitted 
> range, color(red%50)), legend(off)

그러면 아래와 같은 그래프를 얻을 수 있다. 경우에 따라 통제변수를 추가할 수도 있다.