* do-file for exam 2013.1 (same analyses as in listing)
version 13
set more off
cd "h:\vhm\vhm802"

insheet using data_csv\exam13_1.csv /* these data are not available */
encode tx, g(Tx)
gen lnhr=ln(hr)
anova lnhr Tx dog / Tx#dog Tx##time
gen event=Tx*100+dog
xtmixed lnhr Tx##time || dog: || event:, var reml
* recommended form for version 13+
* mixed lnhr Tx##time || dog: || event:, reml