Categorical data analysis becomes challenging when high-dimensional sparse covariates are involved, which is often the case for omics data. We introduce a statistical procedure based on multinomial logistic regression such scenarios, including variable screening, model selection, order selection response categories, and selection. perform our gene expression with 801 patients, 2426 genes, five ...