Stata command review

use http://www.stata-press.com/data/mlmus/hsb, clear

keep if schoolid <=1317

egen meanses=mean(ses) ,by(schoolid)

tab schoolid ,gen(s)

desc

sum mathach, detail

tab female

tabstat mathach ses, by(schoolid) statistics(mean sd)

histogram mathach

twoway (scatter mathach ses) (lfit mathach ses), by(schoolid) legend (order (1 “observed” 2 ”fitted”)) xtitle(ses) ytitle(math achievement)

reg mathach ses

predict infrisk_fit

predict infrisk_stdres, rstandard

version 7

graph infrisk_stdres infrisk_fit, yline(0)

hettest

graph infrisk_stdres, box

qnorm infrisk_stdres

swilk infrisk_stdres

kdensity infrisk_stdres, normal

anova mathach ses schoolid, continuous(ses)

xtmixed mathach sector meanses devses sec_ses mn_ses || schoolid: devses, cov(unstr) mle