jointprobs function

Estimate joint event probabilities for multinomial data