Two-way MANCOVA

R CODE

# two-way mancova
annualincome <- c(50, 100, 120, 70, 250)
studentsloan <- c(200, 300, 350, 120, 100)
educationlevel <- c(2, 3, 2, 1, 2)
gender <- c(1, 2, 2, 1, 1)
parentincome <- c(80, 120, 100, 80, 90, 200, 250, 300, 350, 400)
data <- data.frame(group=c(annualincome, studentsloan), educationlevel=educationlevel, gender=gender, parentincome=parentincome)
manova <- manova(cbind(educationlevel, gender) ~ group + parentincome, data=data)
summary(manova)
##              Df   Pillai approx F num Df den Df Pr(>F)
## group         1 0.197487  0.73826      2      6 0.5168
## parentincome  1 0.043389  0.13607      2      6 0.8754
## Residuals     7