added robust multinomial calculation used in Jie's model (but not in a c...

...lass, darn)