Separation plot for models with binary response variables

data(binary)
ggs_separation(ggs(s.binary, family="mu"), outcome=y.binary)