print_multinomial function

Prints descriptive statistics for multinomial variables