1Department of Biostatistics, School of Public Health2Jiangsu Key Laboratory of Preventive and Translational Medicine for Geriatric Diseases3Center for Genetic Epidemiology and Genomics, Medical College of Soochow University, Suzhou 215123, China4Department of Biostatistics, University of Alabama at Birmingham, Birmingham, AL 35294, USA5Department of Bioinformatics and Genomics, University of North Carolina at Charlotte, Charlotte, NC 28223, USA6Department of Bioinformatics7Department of Cell Biology, School of Biology & Basic Medical Science, Soochow University, Suzhou 215123, China
Checking for direct PDF access through Ovid
Motivation:Large-scale molecular data have been increasingly used as an important resource for prognostic prediction of diseases and detection of associated genes. However, standard approaches for omics data analysis ignore the group structure among genes encoded in functional relationships or pathway information.Results:We propose new Bayesian hierarchical generalized linear models, called group spike-and-slab lasso GLMs, for predicting disease outcomes and detecting associated genes by incorporating large-scale molecular data and group structures. The proposed model employs a mixture double-exponential prior for coefficients that induces self-adaptive shrinkage amount on different coefficients. The group information is incorporated into the model by setting group-specific parameters. We have developed a fast and stable deterministic algorithm to fit the proposed hierarchal GLMs, which can perform variable selection within groups. We assess the performance of the proposed method on several simulated scenarios, by varying the overlap among groups, group size, number of non-null groups, and the correlation within group. Compared with existing methods, the proposed method provides not only more accurate estimates of the parameters but also better prediction. We further demonstrate the application of the proposed procedure on three cancer datasets by utilizing pathway structures of genes. Our results show that the proposed method generates powerful models for predicting disease outcomes and detecting associated genes.Availability and implementation:The methods have been implemented in a freely available R package BhGLM (http://www.ssg.uab.edu/bhglm/).Contact:firstname.lastname@example.orgSupplementary information:Supplementary data are available at Bioinformatics online.