GEE solver for correlated nominal or ordinal multinomial responses
