*
* GRFUNC.RPF
* Example from User's Guide section 3.10
*
open data prsmall.xls
cal(q) 1947:1
data(format=xls,org=obs) 1947:1 1988:1
*
set ydiff = gnp-gnp{1}
set rsum  = rate+rate{1}
set mdiff = m-m{1}
*
set rhos   1 101 = .5+.005*(t-1)
set rss    1 101 = 0.0
do i=1,101
   ar1(noprint,rho=rhos(i)) invest 1950:1 1985:4
   # constant ydiff{1} gnp rate{4}
   compute rss(i)=%rss
end do i
scatter(style=lines,vlabel="Residual Sum of Squares",hlabel="Value of rho",header="Multiple Mode in AR1") 1
# rhos rss