Dear experts,
I want to do a mediation analysis using the type of data like below. exp is the exposure, med is mediator and out is the outcome. Can it be done using SEM? Does it make any sense? Or shall I just do logit with and without med and interpret the results. How to do mediation analysis when all variables are binary?
Thank you in advance.

Code:
* Example generated by -dataex-. For more info, type help dataex
clear
input float(exp med out)
0 0 0
0 1 1
0 0 1
0 0 0
1 0 0
0 1 0
0 1 0
1 1 0
1 1 0
1 0 0
1 0 1
0 1 0
1 0 0
1 1 0
1 0 0
0 1 0
0 0 1
1 0 0
0 0 1
1 1 0
1 1 1
0 1 0
1 1 1
1 0 0
1 1 0
1 0 1
1 0 1
0 1 0
0 0 0
0 0 1
0 0 1
1 0 0
1 0 1
1 1 0
0 0 1
0 0 1
1 1 0
1 0 0
0 1 0
1 1 0
0 0 0
0 0 1
0 1 0
0 0 0
1 1 0
0 0 0
0 0 0
1 1 0
1 0 1
1 0 0
0 0 0
0 1 0
1 1 0
1 0 0
0 1 0
0 0 1
1 1 0
1 0 0
1 1 0
0 1 0
1 0 1
1 1 0
1 0 0
0 0 0
0 1 1
1 0 1
0 1 0
1 0 0
0 1 1
0 0 1
0 0 1
1 0 0
0 1 0
1 0 0
1 0 0
0 0 0
0 0 1
0 0 0
0 0 0
1 0 0
1 0 0
1 0 1
1 1 0
1 0 0
1 1 0
1 1 1
1 1 0
0 1 0
1 0 0
1 0 0
1 0 0
0 1 1
0 1 0
1 1 1
0 0 0
1 1 0
1 1 0
1 1 0
0 0 1
0 0 1
end