Accounting for nine out of ten kidney cancers, kidney renal cell carcinoma (KIRC) is by far the most common type of kidney cancer. In view of limited and ineffective available therapies, understanding the genetic basis of disease becomes important for better diagnosis and treatment. The present studies are based on a single type of genomic data. These studies do not consider interactions between genomic data types and their underlying biological relationships in the disease. However, the current availability of multiple genomic data and the possibility of combining it have facilitated a better understanding of the cancer’s characterization. But high dimensionality and the existence of complex interactions (within and between genomic data types) are the two main challenges of integrative methods to analyze cancer effectively. In this paper, we propose a method to build an integrative model based on Bayesian model averaging procedure for improved prediction of clinical outcome in cancer survival. The proposed method initially uses dimensionality reduction techniques to generate low-dimensional latent features for the predictive models and then incorporates interactions between them. It defines the latent features using principal components and their sparse version. It compares the predictive performance of models based on these two latent features on real data. These models also validate several ccRCC-specific cancer biomarkers previously reported in the literature. Applied on kidney renal cell carcinoma (KIRC) dataset of The Cancer Genome Atlas (TCGA), the method achieves better prediction with sparse principal components model by including latent feature interactions as compared to without including them.