compute_probs function

Estimating marginal and joint probabilities in imputed or synthetic datasets