Multivariate binary responses - advice on regression strategy I would be grateful for advice on how to approach the following situation: 
I have a count variable X and four binary variables A, B, C, D.  The count variable is the independent variable (it refers to the number of adverse experiences in childhood) and the binaries are dependent variables (they refer to certain adverse outcomes in adulthood). A respondent in the dataset can have any combination of outcomes, e.g. A, AC, BCD etc.  I want to measure the strength of the association between the count variable X and the outcomes A, B, C, D conditional on the levels of the other outcomes.  
I’m not sure how best to approach this. Would it be justified to reverse the role of variables and treat the count variable X as the outcome and A-D as predictors? So this would be negative binomial regression (there is overdispersion). In this way the association between X and A (B, C…) would be estimated holding other binary variables constant. But it seems to me that logically it would be dodgy as we would be predicting something that happened earlier with something that happened later. 
Or should I use MANOVA instead (but I’ve read somewhere that the interpretation of results is not straightforward).
Or should I use a generalized linear mixed model (never tried it before) as suggested here https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2798811/ .    
 A: You are making a strong assumption that all the childhood events have equal weight in predicting adult outcomes.  But given that, there are several possible ways to proceed.  Here are three main approaches, one of which you've already mentioned.


*

*Turn the problem backwards to predict the number of childhood events given the outcome status of the 4 events.  Use a semiparametric model so as not to impose a distribution on the count, i.e., proportional odds ordinal logistic model.  The parameters of this backwards model will be hard to interpret but the overall test of association and overall measures of strength of association will be meaningful.  Backwards models, when there is only one original predictor (as in your case) are useful because the extent to which X predicts Y is the same as the extent to which Y predicts X in the purely statistical sense.

*Use a full multivariate model for the 4 binary outcomes.  There are several models from econometrics that will handle this situation.  See Greene's book Econometric Analysis.

*Create a hierarchical ordering of A,B,C,D and assign to each person the worst of the 4 events that happened to them.  Predict this ordinal outcome with a semiparametric ordinal response model.


You didn't mention your sample size but that could be an issue.  At least 96 observations are needed just to estimate a simple single proportion with no covariates.
A: The multivariate probit model might be considered, as described in the Greene book mentioned by Frank Harrell. See also (Lesaffre and Mohlenberghs, 1991 Stat. Med 10, 1391-1403). The idea is to think of a multivariate Normal (4 dimensions) distribution of propensity or tolerance towards each event. You model the multivariate normal mean vector as four functions of the independent variable(s). Estimate the probability of each event given the mean vector via probit link function.
Google the Greene book. You'll find some useful "links".
