* do-file for Additional multivariate exercise 12, VHM 802, Winter 2021
version 16 /* works also with versions 14-15 */
set more off
cd "r:\"

* Manly Example 4.3
import delimited r:\skull.csv, clear
encode period, gen(Period)
foreach var of varlist maximumbreadth-nasalheight {
  oneway `var' Period, tab bon  
  }
manova maximumbreadth-nasalheight=Period
matrix errorsscp=e(E)
matrix list errorsscp
matrix errorcorr=corr(errorsscp)
matrix list errorcorr
mvtest covariances maximumbreadth-nasalheight, by(Period)
mvreg maximumbreadth-nasalheight=i.Period
predict res, eq(#1) residual
margins Period
margins Period, predict(equation(#2))
pwcompare Period, equation(#2) pv mcomp(bon) /* equation 1 is the default */
lincom [#2]1.Period - [#2]2.Period /* equation 1 is the default */

* comparison with univariate analysis
anova maximumbreadth Period
regress
predict res1, res
margins Period