How can I use gsem when I have dependent binary variables?

Using the example in the manual:
Code:
webuse womenwk
Code:
 
 generate selected = 0 if wage < . generate notselected = 0 if wage >= .
Code:
 
 gsem (wage <- educ age L)(selected <- married children educ age L@1,family(gaussian, udepvar(notselected))), var(L@1 e.wage@a e.selected@a)
I have generated a dummy variable (minimumwage) which I want to use in place of wage.

Code:
 
 generate minimumwage = (wage>20)
However, I get an error message when I estimate the equation using minimumwage instead of wage using probit

Code:
 
 gsem (wage <- educ age L, probit)(selected <- married children educ age L@1,family(binomial, udepv > ar(notselected))), var(L@1 e.wage@a e.selected@a) option family() invalid; suboptions not allowed with family bernoulli
What is the right syntax to consider?