# Which statitical modeling in more appropriate when there are several yi? (MANCOVA,Multivariate regression, canonocal correlation)

I have these variables in my data set:

1. a categorical variables with 5 levels (Country)
2. sex
4. different continuous variables as dependent such as y1,y2,y3,...,y25 (I have 25 columns and each column is one dependent variable)

Research question: I want to find the relationship between country and each one of yi by controlling the value of grade in men and women separately.

Questions:

1. Which statistical method I can use? Instead of just having separate analysis for each of the yi, can I analyze them in just one model?

2. I want to have the analysis for men and female separately but instead of split the data set into two subset, can I bring sex as the second independent variables into the model?

3. Can I use MANCOVA, multivariate regression, or canonical correlation? Which one is better?

• What are the 25 dependent variables, and why do you think they're related? All the methods I've ever seen assume a single dependent variable. If you think they're different ways of measuring the same underlying thing, then you might try principal components to construct a single dependent variable. Commented Jun 4, 2017 at 15:42
• Do you perhaps have 25 observations of $y_i$? I find the current specification to be a little odd and am not sure I understand you fully. Commented Jun 4, 2017 at 15:43

### Dimension Reduction?

What do these 25 yi represent? What is their real-world representation? Soil quality? Happiness? Return on investment? Are they all measures of completely different things? If so, then you could treat them all as different dependent variables. If they are related, you should try to do some type of dimension reduction. Depending on what the real-world content of your dependent variables yi are, you could try a principal components analysis or exploratory factor analysis to see if you could reduce your 25 items to a smaller number.

For example, consider the Symbolic Racism 2000 scale. It is 8 items, so in my dataset I have y1, y2, ..., y8. However, each of these items taps into the same theoretical construct: prejudice. The original authors of the scale show that all of those items load on one single factor in both exploratory and confirmatory factor analyses. So what do I do? I just get the mean of the 8 items and use that one variable as my dependent variable. I would be very, very surprised if the 25 items you have do not reduce to some more parsimonious structure than just 25 completely different things that you have measured.

### Univariate Analysis on Reduced Data

Let's say you run an EFA or PCA and see that the 25 items are really all measuring one concept. What most people do is simply average across all 25 items (reverse-scoring when necessary) and run a univariate analysis on that. So:

$\overline{y} = (y_1 + y_2 + ... + y_{25})/25$

Then you would specify the linear regression equation:

$\overline{y} = \beta_0 + \beta_1Country + \beta_2Sex + \beta_3Grade + \epsilon_i$

Except that Country would have a number of dummy-coded variables since it is categorical with more than two levels.

Your question about Sex is confusing. If you are interested in the effect of Sex, then you must keep them in the same analysis. The coefficient for Sex will tell you if there is an influence of Sex on the outcome. Are you thinking that the relationship between some of the other variables and yi will be different depending on Sex? If so, then you can specify and interaction between Sex and the other variables of interest.

### Structural Equation Modeling

Let's say that, after doing dimension reduction, you find that you can reduce all the yi to 3 factors. I would then suggest doing a structural equation model. With these models, you can specify various dependent variables at once. It also let's you specify the influence of each of the other predictors (e.g., Sex, Country, Grade on each of the three factors that are made up of different yi simultaneously). This is an example of what one of these models looks like:

You can see from this that there is one IV and three DVs. Each of these three DVs is made up of 3 measured items each. In theory, you could just specify all 25 variables as their own DV predicted from each of the other three IVs.

The problem with this is the model becomes so large and unwieldy. I would be shocked if you aren't able to reduce your data down to a couple constituent factors or componenets.

### Separate Analyses

Let's say that you really do have 25 completely different dependent variables you wish to analyze (a lot of the decisions will depend on the content of the measures). First, I would ask: Why did you choose to investigate 25 totally different things in the same study? If these dependent variables are different enough to not be considered to load on some smaller number of factors (i.e., y1 is test anxiety, y2 is political conservatism, y3 is neuroticism, etc.), then you could run a bunch of different models separately—I think that is fine, given that the dependent variables aren't related to one another at all.

• @ Mark White , great advice, thanks so much for your time. Some of the dependent variables are correlated and some are not. The dependent variables are related to proteins (around 15 variables ) and fertility factors(around 10 variables). Should I reduce the number of DV which are related to proteins? (an the same for fertility)? Also what do you think about applying MANCOVA? Commented Jun 5, 2017 at 12:50
• also regarding sex, if I want to see the relationship between these variable in men and female population separately, does it mean that sex has some effect? There is a possibility to have different relationship between factor and yi by controlling the covariate in men and female populations. Commented Jun 5, 2017 at 12:59
• I think it would be reasonable to create two DVs: those related to proteins, and those related to fertility. You should try a PCA and see what you get. As for sex, it does not sound like you want to run the models separately (I've found that this is rarely a good idea). Instead, just use sex as a covariate. Commented Jun 5, 2017 at 14:59
• Why should I consider sex as a covariate and not second independent variable? PCA is the great suggestion. But what do you think about MANCOVA or multivariate regression since we have different observations as dependent? Commented Jun 5, 2017 at 15:30
• Having it as a second independent variable is the same thing as considering it as a covariate. I would stick to using some type of path analysis, but that is just my taste, I suppose. I have no experience with MANCOVA. You could also fit a multilevel model where observation is at Level 1 and participant is at Level 2. Each participant would have two observations, which could be coded as Fertility or Protein in another column. Commented Jun 5, 2017 at 15:32

Usually in statistic, you use the notation $$y = f(x)$$ for dependent and independent variables, you seem to have this switched.

But if your $y$, your dependent variable, is a categorial model with more than two classes, a multinomial logistic regression should be your starting point.

• @ Shawn1987 , each yi is continious Commented Jun 4, 2017 at 14:47
• Just from your question, no one can tell you if it is a good advice to treat those $y$'s independently. If you suspect some dependence look into seemingly unrelated regression. Commented Jun 4, 2017 at 15:16
• It looks to me as though the OP has this the right way round and wants to do a multivariate regression using a 25-valued outcome with country, sex and grade as predictors. Commented Jun 4, 2017 at 15:19
• @mdewey, Hi, I have 25 columns and each column refers to one dependent variable. Instead of having 25 separate analysis I am looking for a way to analyze them all together. I mean we can have ancova for 25 times but I am looking for a better way. Commented Jun 4, 2017 at 15:33