Maximum likelihood estimation of mixed effect baseline category logit models for multinomial longitudinal data can be prohibitive due to the integral dimension of the random effects distribution. We propose to use multidimensional unfolding methodology to reduce the dimensionality of the problem. As a by-product, readily interpretable graphical displays representing change are obtained. The methodology can be applied to both nominal and ordinal response variables. Relationships to standard statistical models for multinomial data are presented. Several empirical examples are given to show the merits of the proposed modeling framework.