I am reading P.R. Rosenbaum, Model-based direct adjustment, Journal of the American Statistical Association, 1987 (82), pp. 387-394. However, there is one equation, which I can't understand.
$$ d = \frac{1}{N}\{\sum_{s=1}^{S}\sum_{i=1}^{N_s}\frac{z_{si}r_{si}}{\hat{e}_s} - \sum_{s=1}^{S}\sum_{i=1}^{N_s}\frac{(1-z_{si})r_{si}}{1-\hat{e}_s}\}, $$
where $d$ represents the difference, $\hat{e}_s$ is the probability, $z_{si}$ is the indicator function and $r_{si}$ is a binary response such as yes/no. My question is that How can a binary categorical variable $r_{si}$ be divided by a continuous variable $\hat{e}_s$? For example, how can a categorical variable (sex) male = 1 & female = 0 be divided by 0.23? It does not make any sense to me at all. many thanks in advance.