We propose a sparse Bayesian hierarchical model for the analysis of data including radiomic features for characterization of head and neck squamous cell carcinoma. The proposed model facilitates radiomic feature selection, handling of missing values in key predictors as well as prediction in a unified framework. The fully Bayesian approach enables adequate incorporation of uncertainty arising from various aspects of the inference and prediction procedure. The prediction performance of the model is assessed via cross validation and compared with two frequentist methods. [ABSTRACT FROM AUTHOR]