data(schoolsites)
str(schoolsites)
schools <- cbind(schooldata, schoolsites)
schools.mod <- lm(cbind(reading, mathematics, selfesteem) ~
education + occupation + visit + counseling + teacher +
type + region, data = schools)
car::Anova(schools.mod)
Run the code above in your browser using DataLab