*! version 1.0.0 STB-34 gr20
program define slrplot
	version 5.0

	local y="`1'"
	local x="`2'"
	tempvar yhat

	gph open, saving(slrplt1.gph,replace)

	quietly {
		graph `y' `x', xlab ylab
		gphsave

		regress `y' `x'
		predict `yhat'

		gphdt vpoly `y' `x' `yhat' `x'

		sum `x'
		local x1=_result(5)
		local x2=_result(6)
		local y1=_b[_cons]+_b[`x']*`x1'
		local y2=_b[_cons]+_b[`x']*`x2'

		gphdt line `y1' `x1' `y2' `x2'
	}

	gph close
end