DF = read.csv("../../Data/swedish_hospitals.csv")
DF$Hospital = as.factor( DF$Hospital )

fit = aov( Bill ~ Hospital, data=DF )
sfit = summary(fit)
print( sprintf("Classical ANOVA (Hospital ME): prob. type I error: %10.6f", sfit[[1]][["Pr(>F)"]][1]) )