* do-file for Exercise 26.37 of PSLS 3e
version 13 /* works also with version 14 */

import delimited "H:\VHM\VHM801\Datasets\csv\Chapter 26\ex26_037.csv"
table rate hybrid, contents( mean percent sd percent)

encode hybrid, gen(Hybrid)
egen group_rate=group(rate)
anova percent group_rate Hybrid group_rate#Hybrid
predict pred, xb
sort Hybrid rate
overlay pred rate, by(Hybrid) c(l)
* note: overlay is an add-on command (findit overlay)

predict stdres, rstandard
scatter stdres pred
qnorm stdres
tabulate stdres