Hello,

I have a sample of 1239 participants. My exposure variable is binary (1 = yes, 0 = no), and this is experience of psychosis (168 said yes). I have a number of outcomes, and I plan to use separate logistic regression models within a generalised estimating equations (GEEs) to get odds ratio for each outcome (example of such outcomes are smoking, alcohol, BMI etc). Most outcomes are binary.
I have two covariates that I do not think are appropriate to simply adjust for (binary - depression score and anxiety score). I am wondering how best to know how to proceed. Do I do a stratified analyses by looking at whether my exposure is associated with any of my outcomes within the depressed group and/ anxious group? How do I find out whether I have power for this? Or do I do a test of co-linearity?
Thanks