Identifying differentially expressed genes using the Polya urn scheme
Erlandson Ferreira Saraiva, Adriano Kamimura Suzuki, Luís Aparecido Milan · Communications for Statistical Applications and Methods · 2017
A common interest in gene expression data analysis is to identify genes that present significant changes in expression levels among biological experimental conditions.In this paper, we develop a Bayesian approach to make a gene-by-gene comparison in the case with a control and more than one treatment experimental condition.The proposed approach is within a Bayesian framework with a Dirichlet process prior.The comparison procedure is based on a model selection procedure developed using the discreteness of the Dirichlet process and its representation via Polya urn scheme.The posterior probabilities for models considered are calculated using a Gibbs sampling algorithm.A numerical simulation study is conducted to understand and compare the performance of the proposed method in relation to usual methods based on analysis of variance (ANOVA) followed by a Tukey test.The comparison among methods is made in terms of a true positive rate and false discovery rate.We find that proposed method outperforms the other methods based on ANOVA followed by a Tukey test.We also apply the methodologies to a publicly available data set on Plasmodium falciparum protein.