Hi all,

I have an issue trying to estimate multiple outcomes simultaneously using hierarchical models in Stata (I am using version 15.1).

The research question is looking at different health outcomes (eg. anxiety and depression levels) and whether some patient characteristics (eg. awareness of their disease and prognosis) affects their health outcomes. The patients of the study are from multiple countries, as such, we want to use a hierarchical model to nest the individuals within the country they are from.

The issue now arises as we are trying to test multiple outcomes simultaneously (using meglm or mixed). I am unable to find a method for doing so in Stata, and would like to ask if anyone is familiar with such estimation methods in Stata.

Thank you.