In many biomedical applications, the response consists of multiple categorical random variables and the predictor is high-dimensional. In this talk, we discuss new models and methods for modeling such ...