Comment 1 for bug 422216

I'm afraid the description is imprecise. Instead of 2d data it should be 2d with shape[1] = 2. This is for binomial data that has a different number of trials for each observation and so the response data is (success, failure) and the argument data_weights then needs to be used in fit with data_weights set to success+failure. There should be an example of this in the examples directory.

Unfortunately, there is not support for multinomial data yet. There was some discussion about extending this, and someone even offered some code, but I just haven't had the time to extend it. I think it would be pretty straightforward to implement though.

Skipper