基于DNA微阵列数据分析的分级Bayes模型

如何分离出少量区别不同组织类型的特异性基因是DNA微阵列数据分析中的主要问题,特别是构建恰当的统计模型来刻画这些不同组织类型的DNA表达形式尤为重要.为此,基于基因DNA微阵列数据的特点,我们假定对数变换后的微阵列数据服从混合正态分布.我们采用分级Bayesian先验刻画不同基因的相关性,利用分级Bayesian方法构建模型,给出了刻画不同组织基因表达的差异的一个标准,用MCMC迭代计算该标准.模拟计算表明我们的模型具有较好的识别能力.

A Bayesian Hierarchical Model for DNA Microarray Data

In DNA microarray analysis, there is often interest in isolating a few genes that best discriminate between tissue types. In particular, it is critical to develop suitable models to explain the patterns of DNA expression for these different types of tissues. Toward this goal, we propose a methodology for the analysis of high-density oligonucleotide arrays.The log-transformed data are assumed to follow a mixture normal distribution based on the characteristic of gene itself. The variation in the data can reasonably be thought to arise from the effects of genes, tissue types, and their interactions. We introduce a hierarchical Bayesian priors for the parameters and propose a model selection criterion for identifying subsets of genes that show different expression levels between normal and tumor types. In addition, we develop Markov chain Monte Carlo algorithms for sampling from the posterior distribution of parameters and for computing criterion. The proposed methodology is evaluated via simulations studies.