I'm trying to apply the Oaxaca-Blinder decomposition to the difference of two probabilities: Prob(y=1|A=1) for group A and Prob(y=1|B=1) for group B with y=(0,1) and Prob(y=1)=1/(1+e^-XBeta). My aim is separating the share of the probabilities difference due to cross-group differences in the explanatory variables from the part that remains unexplained by these differences (due to discrimination). I'd like to know if there is a package in R that perfoms that operation.
Thanks in advance for your help.