computep function

Computes individual and pairwise treatment assignment probabilities.