Possibly using [statsmodels.formula.api.mnlogit](https://www.statsmodels.org/stable/generated/statsmodels.formula.api.mnlogit.html#statsmodels.formula.api.mnlogit) Currently categorical phenotypes are not supported. It will be important to have results in R that can be compared for accuracy.