#4-#6 on Homework 3

Here are the codes I used for #4-#6 on Homework 3 (you have to do #3 to do #4 so I included it). I think the graph in #6 could be better, but it is a start.

3.
. generate logdol = log(dollars)

. generate logincid = log(incid)
(4 missing values generated)

. generate logprev = log(preval)
(4 missing values generated)

. generate loghospdays = log(hospdays)

. generate logmort = log(mort)

. generate logyrslost = log(yrslost)

. generate logdisbil = log(disabil)

. sw regress logdol logincid logprev loghospdays logmort logyrslost logdisbil, forward pe(.1) pr(.2)


4. . sw regress logdol logincid logprev loghospdays logmort logyrs
> lost logdisbil, pe(.1) pr(.2)

5. regress logdol loghospdays


6. regress logdol loghospdays logmort logyrslost logdisbil


. predict yhat, xb

.predict h, leverage

. predict std_yhat
(option xb assumed; fitted values)

. predict std_f, stdf

. generate ciu_f = yhat + invttail(25-5-1, .025)*std_f

. generate cil_f = yhat - invttail(25-5-1, .025)*std_f

. generate cil_sbpf = exp(cil_f)

. generate ciu_sbpf = exp(ciu_f)

. display invttail(19, 0.025)
2.0930241

. predict t, rstudent

lowess t yhat, bwidth(0.8) mcolor(gs10) clwidth(thick) yline(-1.96 0 1.96) x
> label(4.7(.1)5.1)