3-group comparison on more than 1 outcome variables

Hello,

Being an RStatsnewbie, I’ve spent the last 3 hours going through various ANOVA exercises. Unfortunately, everyone I’ve seen uses 1 outcome variable and 1 or 2 grouping/categorical variables.
For example, I have no issues with:

model1 = aov(bdi_total, ~ group.factor, data = dat1)
TukeyHSD(model1)

However, I have more than 1 outcome variable.
My data:
group.factor = 3 (HC, FM, CLBP)
outcome variables = 3 (bdi_total, pdi_total, psqi_total)

All I want to do right now is to test whether for a group difference on the 3 outcome variables at one time. I know the answer will no be overly complicated, but I honestly cannot find an appropriate example.

Jason

Do you want a multivariate analysis of variance? If so, this example shows how it works using one of the R in-built datasets:

data(mtcars)
head(mtcars)
help(manova)

model1 <- manova(cbind(mpg, qsec, wt) ~ cyl * gear * carb, data = mtcars)
summary(model1)
model1

HTH

2 Likes

:man_facepalming: YES! Thank you.
I was having a total brain fart and was fixated on ANOVA, completely forgetting MANOVA.
Color me embarrassed.

Jason

This topic was automatically closed 7 days after the last reply. New replies are no longer allowed.