use hernes,clear
des
summ
hernes age rate,method("hernes")
list
return list

use hernes,clear
hernes age rate,method("nl")
list

gr rate rate_hat age, c(.l) s(Oi) xlab(15(1)38) ylab(0(.1)1) gap(6)saving (hernes,replace)


use hernes,clear
hernes age rate if age<=24