Hi all,

I am trying to build logits to measure gender employment gaps among lone parents in Canada in Feb 2020, in the Labour Force Survey data. I've created a dummy variable for employment (1=employed, or absent; 0=unemployed or not in labour force), and want to measure the difference (gap) in mean employment rates between subgroups.

What syntax can I use to capture the difference in proportions of employment between male and female lone parents? How can I compare this gender gap among lone parents, e.g., within two sets of (ordinal) sub-groups: 1) older or younger child (<6 or 6-12) and 2) education?

Please below an excerpt from dataex:

Code:
* Example generated by -dataex-. For more info, type help dataex
clear
input float lfs byte sex float(loneyg edu)
0 1 . 1
1 1 . 1
1 1 . 0
1 1 . 2
1 1 . 1
0 1 . 2
1 1 . 0
1 1 . 1
1 1 . 1
1 1 . 0
1 1 . 1
1 1 . 2
1 1 . 1
1 1 . 2
1 1 . 2
0 1 . 0
1 1 . 0
1 1 . 0
1 1 . 1
1 1 . 0
1 1 . 1
1 1 . 2
1 1 . 2
0 1 . 1
1 1 . 1
0 1 . 0
1 1 . 1
1 1 . 0
1 1 . 1
1 1 . 2
1 1 . 1
0 1 . 1
1 1 . 1
1 1 . 1
1 1 . 0
1 1 . 0
1 1 . 1
1 1 . 1
1 1 . 1
0 1 . 0
1 1 . 0
1 1 . 2
1 1 . 2
1 1 . 1
1 1 . 0
1 1 . 0
0 1 . 0
1 1 . 1
1 1 . 1
0 1 . 0
1 1 . 1
1 1 . 0
1 1 . 2
1 1 . 0
1 1 . 1
1 1 . 0
1 1 . 1
1 1 . 1
1 1 . 0
1 1 . 1
1 1 1 1
1 1 . 2
1 1 . 0
0 1 . 0
1 1 . 1
0 1 . 0
1 1 . 2
1 1 . 1
1 1 . 1
1 1 . 1
1 1 . 1
1 1 . 2
1 1 . 2
1 1 . 0
1 1 . 2
1 1 . 2
1 1 . 1
1 1 . 1
1 1 . 1
0 1 . 0
1 1 . 2
1 1 . 2
1 1 . 0
1 1 . 1
0 1 . 1
1 1 . 1
1 1 . 1
1 1 . 0
1 1 . 2
1 1 . 1
1 1 . 1
1 1 . 1
1 1 . 2
1 1 . 0
0 1 . 2
1 1 . 0
1 1 . 2
1 1 . 2
1 1 . 1
1 1 . 0
end
label values lfs lfs
label def lfs 0 "not", modify
label def lfs 1 "Employed", modify
label values sex SEX
label def SEX 1 "Male", modify
label values loneyg loneyg
label def loneyg 1 "Lone parents, yg child", modify
label values edu edu
label def edu 0 "(<)HS", modify
label def edu 1 "some uni/college deg/trades", modify
label def edu 2 "BA degree+", modify