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

import delimited "R:\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