data(Iwasaki_Big_Five)
str(Iwasaki_Big_Five)
mod <- lm(cbind(N, E, O, A, C) ~ Group, data=Iwasaki_Big_Five)
car::Anova(mod)
# heplots
labs <- c("Neuroticism", "Extraversion", "Openness", "Agreeableness", "Consientiousness" )
heplot(mod,
fill = TRUE, fill.alpha = 0.2,
cex.lab = 1.5,
xlab = labs[1], ylab = labs[2])
heplot(mod, variables = c(2,5),
fill = TRUE, fill.alpha = 0.2,
cex.lab = 1.5,
xlab = labs[2], ylab = labs[5])
pairs(mod,
fill = TRUE, fill.alpha = 0.2, var.labels = labs)
Run the code above in your browser using DataLab