Dear SAS users,
I am doing a mediation analysis using cross-sectional data (n=7,173):
X= family tension (binary)
M= sleep disturbance (continuous)
Y= Depression (binary and rare outcome in my population)
There is no interaction between X (family tension) and M (sleep disturbance)
The equations are the following:
Logit {P(depression=1|a,m,c)}=θ0+ θ1family_tension+ θ2sleep_disturbance+θ3family_tension*sleep_distrubance+θ’4covariates
E[sleep_disturbance|family_tension,covariates]=β0+ β1dispute + β’2covariates
When I run my mediation analysis using @PROCCAUSALMED, I obtain the following results:
Total effect (OR= 4.05, 95%CI, 3.06-5.05)
Natural direct effect (OR=2.39, 95%CI, 1.88-2.90)
Natural indirect effect (OR=1.70, 95%CI, 1.48-1.91).
If I run a logistic regression between X and Y (without the mediator and including covariates), I would expect to get the OR=4.05 as in the CAUSALMED result of the total effect, however I get a OR= 3.06, 95%CI, 2.49-3.76).
Does anyone know why I am not getting the same OR? I would expect to also get OR=4.05 from a logistic regression (x->y without the mediator), because I had understood the total effect in the causal mediation would be equivalent to the reduced model using traditional mediation (x->y without the mediator).
Perhaps I am mixing concepts, your help is very wellcome. Thank you so much in advance,
Ana
