A full Bayesian hierarchical mixture model for the variance of gene differential expression
 Samuel OM Manda^{1}Email author,
 Rebecca E Walls^{1, 2} and
 Mark S Gilthorpe^{1}
DOI: 10.1186/147121058124
© Manda et al; licensee BioMed Central Ltd. 2007
Received: 11 October 2006
Accepted: 17 April 2007
Published: 17 April 2007
Abstract
Background
In many laboratorybased high throughput microarray experiments, there are very few replicates of gene expression levels. Thus, estimates of gene variances are inaccurate. Visual inspection of graphical summaries of these data usually reveals that heteroscedasticity is present, and the standard approach to address this is to take a log_{2} transformation. In such circumstances, it is then common to assume that gene variability is constant when an analysis of these data is undertaken. However, this is perhaps too stringent an assumption. More careful inspection reveals that the simple log_{2} transformation does not remove the problem of heteroscedasticity. An alternative strategy is to assume independent genespecific variances; although again this is problematic as variance estimates based on few replications are highly unstable. More meaningful and reliable comparisons of gene expression might be achieved, for different conditions or different tissue samples, where the test statistics are based on accurate estimates of gene variability; a crucial step in the identification of differentially expressed genes.
Results
We propose a Bayesian mixture model, which classifies genes according to similarity in their variance. The result is that genes in the same latent class share the similar variance, estimated from a larger number of replicates than purely those per gene, i.e. the total of all replicates of all genes in the same latent class. An example dataset, consisting of 9216 genes with four replicates per condition, resulted in four latent classes based on their similarity of the variance.
Conclusion
The mixture variance model provides a realistic and flexible estimate for the variance of gene expression data under limited replicates. We believe that in using the latent class variances, estimated from a larger number of genes in each derived latent group, the pvalues obtained are more robust than either using a constant gene or genespecific variance estimate.
Background
The recent advancement of deoxyribonucleic acid (DNA) microarray technology allows the measurement of expression levels of tens of thousands of genes simultaneously [1, 2]. A DNA microarray experiment measures the abundance of messenger ribonucleic acid (mRNA) present in a set of cells, and a high concentration of mRNA for a given gene indicates a high expression level for that gene [3, 4]. The solution of mRNA is either radiolabelled or fluorescently labelled and then allowed to hybridize to spots on the array. Further information about early DNA microarray experiments can be found in [5]. There is a wide variety of arrays, but the two main kinds are short and long oligonucleotide arrays. In oligonucleotide arrays, there are approximately 20 probes pairs of a perfect match (PM) and a mismatch (MM) for each gene. The PM probe contains a match to a small subsequence of a gene's polynucleotide, about 25 bases long, while the MM probe acts as a control, being a copy of the PM probe with the central position flipped to its complement. The amount of mRNA present in the target gene in a sample is estimated from a combination (usually, the average) of the PMMM intensity differences over the 20 probe pairs [6].
In a spotted complimentary DNA (cDNA) or long oligonucleotide arrays (the kind featured in this article) thousands of spots of cDNA from the genes are printed onto a glass slide or some other form of substrate. Two different mRNA samples are separately reversetranscribed into cDNA and labelled with different fluorescent dyes, green (cyanine 3) or red (cyanine 5). The mixture of labelled cDNA is cohybridized on the same microarray, and the labelled cDNA molecule will bind to the complementary fragments of cDNA sequence on the slide. A laser scanner is then used to measure both fluorescent signals emitted at each spot on the chip. The general idea assumed behind the technology is that if a particular gene is highly expressed in the sample, it produces many molecules of mRNA [4]. These in turn will hybridize to the probes on the microarray and generate a very bright fluorescent area. Genes that are less expressed produce less mRNA, which results in dimmer fluorescent spots. If there is no fluorescence, no messenger molecules have hybridized to the probes, indicating that the gene is inactive. By comparing the intensity levels of the emitted fluorescent lights between the samples, it is hoped that one might be able to identify any differences in the gene expression profiles of the various samples.
In a typical cDNA microarray experiment, we are looking to ascertain whether gene i displays differential expression between two samples T and C, labelled with differential colours red and green. For instance, these experiments include comparing tumor and normal tissue cells, treated and untreated cells, or cells from knockout and wildtype organisms. The samples can be compared on the same slide (i.e. same hybridization), resulting in a direct estimate of differences in expression levels since the measurements come from the same hybridization. An alternative is when expression levels T_{ i }and C_{ i }are estimated in two different hybridizations, with T together with reference sample R and C with another reference sample R'. This is an indirect estimate of the gene's differential expression since the T and C expression levels are from different hybridizations [7]. Sometimes, common reference samples are hybridized with each mRNA sample of interest (T or C), resulting in what are known as common reference designs. The common reference sample could be tissues from wildtype organisms or control tissue or a pool of all the samples of interest.
Standard statistical analyses
The simplest way, used by many in the field, to ascertain a gene's differential expression, is on the basis of a foldchange criterion, defined by the logratio δ_{ i }= log_{2}(T_{ i }/C_{ i }), under direct comparison or δ_{ i }= log_{2}(T_{ i }/R_{ i })  log_{2}(C_{ i }/${{R}^{\prime}}_{i}$) for indirect comparison. It is expected that the majority of genes will have a δ_{ i }value close to 0 [4]. Those genes with a large positive δ_{ i }value (δ_{ i }> 1) are generally concluded to be overexpressed or upregulated in the T sample, and those whose δ_{ i }is negatively large (δ_{ i }< 1) are concluded to be underexpressed or downregulated in the T sample. However, the use of fold change is of limited use, as the intensities are associated with some biological, experimental and measurement error, and unless these error distributions can be derived, it is difficult to assess whether a ratio of 1.9, say, is worth noting or whether it has occurred by chance. Furthermore, the boundaries accepted for thresholding these foldchanges seem to be rather arbitrary and very little documentation can be found to support these criteria.
In recent times, the identification of genes that are differentially expressed between two conditions has been based on a standardised foldchange, which is the foldchange divided by an estimate of its standard deviation. A ttest, with a correction for multiple testing, is then used to test for significance of the standardised foldchange. This allows for the assessment of significance of the observed differences in the presence of all the sources of variation, which are not necessarily equal from gene to gene. Our contribution to the problem of identifying genuinely differentially expressed genes is on the estimation of reliable standard deviations of gene expression levels. Modelling of gene expression variability ranges between two extreme cases: a constant variance, which is too stringent an assumption, to independent genespecific variances. The latter option has low power as it is based on very few replications as a result of the relatively large cost of commercial microarray chips. This makes estimation of the sample standard deviation very unreliable and unstable. An adhoc solution to the problem includes discarding genes with a small foldchange and very small standard deviations [4]. A better method, called Significance Analysis of Microarrays (SAM), developed in [8], adds a constant a_{0} to each genespecific standard deviation, thus preventing the denominator of the tstatistic from getting too small. This was expounded in [3], using an empirical Bayesian procedure, by taking a_{0} equal to the 90th percentile of the standard deviations of all the genes. Shrinking the genespecific standard deviations in this way helps to minimize the false discovery rates (i.e. a large t statistic).
Mixture models
An intuitive approach to modelling gene expression data is to assume two groups of genes, one group with genes that are differentially expressed and the other with genes that are not differentially expressed. This approach has been used in analyses involving mixture models for gene expression levels. Formally, the basic assumption is that the distribution of the difference δ_{ i }can be flexibly modelled as a mixture with two components: a subgroup of genes with δ_{ i }around 0 and a subgroup of genes with non zero δ_{ i }[3, 4, 9, 10]. Using this approach, Lonnstedt and Speed [4], derived an empirical Bayesian statistic B, which is the log posterior odds of differential expression. Efron et al. [3] and Efron and Tibshirani [10] also consider a two component mixture model to model differential gene expression, the later using a rankbased nonparametric twosample test statistic. A similar approach was followedup in [11] based on a fully Bayesian hierarchical model, but with an unknown number of mixture components. The number of components was treated as a random variable and estimated with the other parameters based on the pioneering work of Richardson and Green [12]. Other approaches using mixture models for gene expression data from microarray experiments can be found in [13, 14].
Purpose of the paper
A number of methods, particularly based on full Bayesian hierarchical models, have been used to provide better estimates of variance for gene differential data. These methods provide estimates of genespecific variance, which are the weighted average of the empirical variance and a prior variance estimate [11, 15, 16]. In particular, Lewin et al. [16] provide a fully Bayesian approach combining estimation of gene differential expression, biological variability and array effects with a hierarchical prior distribution on genespecific variances. Other than modelling the genespecific variances with an exchangeable hierarchical prior, Delmar et al. [17] use a mixture model on the distribution of genespecific variances. Genes are grouped into latent classes based on homogeneity of their variances. A gene is assigned a variance based on its latent class membership and this variance is estimated with high power because of the large number of genes (hence a larger number of replicates) in that latent class. All these methods produce what are called regularized ttests.
Our work is closely related to that proposed by Delmar et al. [17], who used the EMalgorithm to fit a variance mixture model for gene expression data. We believe our Bayesian approach has certain advantages and adds value in comparison to the EMalgorithm approach. In using a Bayesian hierarchical model, we are able to model various sources of variability in a common model, thus propagating uncertainty. Within a Bayesian hierarchical model framework, it is far easier to borrow and share data across all genes in order to obtain more reliable estimates of their variance and at the same time allowing for some variability. In this approach, variances are stabilized and shrunk towards the average variance within each component of the mixture, in particular some small and large variance estimates that are incompatible with the overall distribution are increased and decreased, respectively.
Furthermore, in complex biological data exhibiting a lot of noise, traditional statistical methods, such as the EMalgorithm, can struggle to cope with complex nonlinear models when used to explore such data. In the Bayesian paradigm, on the other hand, all the unknown quantities are treated together in a consistent manner, to give fully probabilistic information on all unobserved variables, even their functions. Our method is based on a Chisquared (χ^{2})mixture model for the genespecific variances, with the number of components ranging from 1 to 5. In bioinformatics, as in many fields, mixture models have been fitted through the expectationmaximisation (EM) algorithm with different values of the number of mixture components [9, 17, 18]. In this paper, we present a Bayesian analysis of the variance mixture model, which we implement in the Bayesian software package, WinBUGS [19]. The methodology is applied to a dataset on diffuse large Bcell lymphomas. The data contains expression levels of chronic lymphatic leukaemia (CLL) and diffuse large cell (DLCL) malignancies for the 9216 genes under study. The results of applying our model to the dataset are presented in the Results from the mixture model section. A discussion of the model and the results are in the Conclusions section.
It is hoped that the work presented in this paper will contribute to a large volume of current research work aimed at minimizing the risk of false positives in microarray experiments. The two extreme assumptions on the variance of gene expression data are presented in the Standard models for gene variance section and a description of the variance mixture model is contained in A mixture model for gene variance section.
Results and discussion
Description of the data set
We analyse the data described in [20], which investigates the classification of diffuse large Bcell lymphomas into distinct groups by gene expression profiling. Diffuse large Bcell lymphoma is an aggressive malignancy of mature B lymphocytes. It is estimated that, with an annual incidence of over 25,000 cases, it accounts for approximately 40% of all cases of nonHodgkin's lymphoma. Currently, a combination of clinical parameters is used to provide an assessment of a patient's risk profile and to determine the most suitable clinical course of treatment. Whilst most patients initially respond well to chemotherapy, patients receiving the same diagnosis can have very different final outcomes in terms of remission achieved and their overall survival. It is suspected that the prognostic variables used are in fact proxies for the underlying cellular and molecular variation within diffuse large Bcell lymphomas. Their work considers whether gene expression profiling could subdivide this 'clinically heterogenous diagnostic category into molecularly distinct diseases that would possess more homogeneous clinical behaviours' [20]. The microarrays used in this experiment were specially designed complementary DNA microarrays, called 'Lymphochips', which included those genes with a known or suspected role in processes that are important in immunology or cancer, and those genes known to be preferentially expressed in lymphoid cells. The profiling of gene expression included the three most prevalent adult lymphoid malignancies, in addition to profiling the gene expression in purified normal lymphocyte subpopulations under a range of activation conditions, in normal human tonsil and lymph node, and in a variety of lymphoma and leukaemia cell lines. From each experimental mRNA sample, a cDNA sample was prepared and labelled with red Cy5 dye. Furthermore, a corresponding reference cDNA sample, labelled with green Cy3 dye, was prepared from a pool of mRNAs isolated from nine different lymphoma cell lines. The labelled samples were combined and hybridized to the microarray.
Approximately 1.8million measurements of gene expression were taken in all, across 96 normal and malignant lymphocyte samples, using 128 of the Lymphochip microarrays. To demonstrate the method proposed in this paper, we use a small subset, containing only eight slides from two conditions. Four slides quantify gene expression relating to chronic lymphatic leukaemia (CLL) malignancies and the other four to diffuse large cell lymphoma (DLCL) malignancies. Each slide contains measurements for 9216 genes. The redtogreen intensity ratio can be quantified for each gene and this reflects the relative abundance of mRNA in the experimental sample compared with the reference mRNA pool. By using a common reference sample, these fluorescent ratios can be considered a comparable measurement of the relative expression level of each gene across all of the samples. We want to compare different gene expression levels between the CLL and DLCL malignancy conditions.
Results from the mixture model
where ${y}_{ic{r}_{c}}$ is gene i's logratio of observed intensity for condition c to the reference mRNA pool in replicate r_{ c }(i = 1, ..., 9216; c = 1 for CLL, 2 for DLCL; r_{1} = 1, ..., 4; r_{2} = 1, ..., 4); μ_{ic} is the mean logratio; and ${\epsilon}_{ic{r}_{c}}$ is normally distributed with mean 0.
The mixture weights were estimated as part of the model, where they were assigned a Dirichlet prior distribution, (π_{1}, ..., π_{ k }) ~ Dirichlet(1, ..., 1), where k is the number of mixture components. The χ^{2} mixture scale parameters ψ_{ j }were assigned independent Gamma(0.01, 0.01) prior distributions. For each mixture model, three independent chains were run for 50,000 iterations. We discarded the first 20,000 iterations and used a combined sample of the remaining 60,000 iterations for posterior summaries.
Posterior mean (SD) for various variance mixture models
Model  Bayesian Model  EM algorithm  

π_{ j }  $\overline{\sigma}$ _{ j }  π_{ j }  $\overline{\sigma}$ _{ j }  
2 Classes  
1  0.2276 (0.009)  0.7624 (0.008)  0.2286  0.7607 
2  0.7724 (0.0085)  0.3426 (0.0024)  0.7713  0.3422 
BIC  22911.89  22657.89  
3 Classes  
1  0.0677 (0.0039)  1.0220 (0.2627)  0.0688  1.0225 
2  0.4462 (0.0170)  0.4977 (0.0077)  0.4546  0.4970 
3  0.4860 (0.0182)  0.2841 (0.0035)  0.4766  0.2839 
BIC  21819.82  21703.82  
4 Classes  
1  0.2303 (0.0146)  0.6125 (0.0129)  0.2317  0.6125 
2  0.0384 (0.0045)  1.1547 (0.0321)  0.0370  1.1523 
3  0.2318 (0.0218)  0.2413 (0.0053)  0.2285  0.2427 
4  0.4994 (0.0201)  0.3775 (0.0161)  0.5028  0.3782 
BIC  21727.75  21620.86 
where _{m1} and _{m2} are sample sizes for condition 1 and 2, respectively.
Classification of genes under the Bayesian and EM mixture models, with four components
EM model  Bayesian model  Total  

1  2  3  4  
1  1621  3  0  52  1676 
2  18  242  0  0  260 
3  0  0  1848  161  2009 
4  42  0  111  5118  5271 
Total  1681  245  1959  5331  9216 
Top ten ranked genes by different variance models (Genes are listed by their codes)
Genespecific  Constant  2 Classes  3 Classes  4 Classes  Weighted 

4323  2602  1939  2143  2143  2143 
4532  1945  2143  3181  4069  4069 
4069  257  3181  4069  4323  4323 
8076  1400  7161  8903  4532  4532 
4331  1939  2003  4323  7496  8076 
6635  2143  4069  4347  8076  3181 
2026  3181  6731  4532  4586  4586 
2143  7003  6343  6592  5674  8903 
4586  4323  8649  8628  8048  2542 
8892  3151  4323  5746  1631  5674 
Finally, other than setting ${\widehat{\sigma}}_{i\ast}^{2}={\widehat{\sigma}}_{j}^{2}$, which is outright membership, we also considered a weighted variance estimate ${\widehat{\sigma}}_{iw}^{2}={\displaystyle \sum {\pi}_{ij}{\sigma}_{j}^{2}}$ based on a fourcomponent variance mixture model. We ranked genes using this weighted variance estimate to compare the ranking obtained from the other variance estimates (see Table 3). Eight genes coded 2143, 3181, 4069, 4323, 4532, 4586, 8076 and 8903 that appear in the top ten under the weighted variance model, also appear at least once in the top ten of the other variance models. These genes might be interesting, thus requiring further analysis and investigation.
We are aware that different assumptions for the prior distributions for both mixture weights and scales may give different results. We performed some limited sensitivity examination of the results to different specification of the priors. There were slight differences in the results, but the substantive conclusions were not affected.
Conclusion
We have presented a Bayesian variance mixture model for differential gene expression data. This model is a compromise between two extreme models: the too stringent constant gene variance and the overparameterised genespecific variance models, which are both unrealistic assumptions. Our mixture variance model provides a more realistic and flexible estimate for the variance of gene expression data under limited replicates. We believe that in using the (weighted) latent class variances, estimated from a larger number of genes in each derived latent group, the pvalues obtained are more accurate then either using a constant gene variance or genespecific variance estimate.
In our example data, the results obtained from using our model are in close agreement to those obtained using EM algorithm implementation [17], which had been shown, using simulation studies, to be flexible and reliable in both true and false discovery rates. Our results are based on the assumptions of normally distributed logratios and a constant genevariance between the two conditions. We are working on relaxing these conditions, in particular in using a longtailed tdistribution as a robust alternative to allow for the possibility of geneintensity measurement outliers. By varying the degrees of freedom, the tdistribution can also be used to investigate the sensitivity of the posterior results to changes in the prior for the gene intensity measurements.
Methods
Notation and variance models
We assume that the intensity level data is background corrected and normalised according to [23], using an arcsinh transformation based on a model for the dependence of the variance on the mean intensity levels with variance stabalizing properties. These are implemented in a variance stabilizing transformation vsn, the title function of the vsn package, part of the Bioconductor project for R [24]. This results in generalised logtransformed expression intensity values. We develop the methodology for unpaired (indirect) data case, where two samples of interest are each cohybridised with a reference sample. That is, each independent slide is a twocolour microarray experiment. The methodology is easily adapted to paired (direct comparison) data case.
The selection of differentially expressed genes can proceed simply by a test based on log foldchange criterion, ${\overline{\delta}}_{i}={\overline{T}}_{i}{\overline{C}}_{i}$. However, as pointed out in the Standard statistical analyses section, the use of fold changes is limited because the intensities are associated with biological, experimental and measurement errors. Formally, we consider significance testing by assuming that the observed logratios y_{ itr }and y_{ itr' }are modelled by simple linear models:
y_{ itr }= μ_{ iT }+ ε_{ itr }and y_{ icr' }= μ_{ iC }+ ε_{ icr' }
whose usefulness depends on an accurate estimate of the standard error (SE) of $\overline{\delta}$_{ i }. Generally, there will be very few replications; thus the estimation of SE would be very unreliable and unstable. One solution is to discard genes with a small foldchange and very small standard deviations to avoid getting overoptimistic significant results.
Standard models for gene variance
There are two extreme cases to model this variance:

Independent genespecific variance: In this scenario, ${\sigma}_{iT}^{2}={\sigma}_{iC}^{2}={\sigma}_{i}^{2}$ is estimated by:${S}_{i}^{2}=(({M}_{T}1){S}_{iT}^{2}+({m}_{C}1){S}_{iC}^{2}))/{m}_{T}+{m}_{C}2)$(7)
which has ν = m_{ T }+ m_{ C } 2 degrees of freedom. Under the null hypothesis, μ_{ iT }= μ_{ iC }, the standardised foldchange t_{ i }is distributed according to a tdistribution with ν degrees of freedom. This option has low power as it is based on very few replications. This makes estimation of the sample variances very unreliable and unstable, and this results in less powerful ttests.

A constant variance: In this homoscedastic model, all genes are assumed to have the same variance σ^{2}, which is estimated by:${S}^{2}=\frac{1}{g}{\displaystyle \sum _{i=1}^{g}{S}_{i}^{2}}.$(8)
This has a large number of replicates, totalling (m_{ T }+ m_{ C })g. The degrees of freedom of the variance estimate are now ν = (m_{ T }+ m_{ C } 2)g, which makes the statistic t_{ i }behave like a standard Normal (0,1) variate. However, a constant variance over all the genes is too unrealistic an assumption, and it increases the risk of a false positive result for a gene with a larger variance. On the other hand, there is a high risk of missing out on a truly differentiated gene having a small variance, and a large differential effect.
A mixture model for gene variance
where k is the number of latent components (classes); π = (π_{1}, ..., π_{ k }) are the mixture proportions with π_{ j }being the probability that a gene belongs to latent class j (π_{1} + ⋯ + π_{ k }= 1); ψ = (ψ_{1}, ..., ψ_{ k }) are the componentspecific χ^{2} distribution scaling constants and g(w_{ i }ψ_{ j }, ν) are χ^{2} distributions with the scaled constant ψ_{ j }being specific to component j and ν = m_{ T }+ m_{ C } 2 is the common degrees of freedom, common for all components. The variance of the class j is ${\sigma}_{j}^{2}$ = 1/ψ_{ j }, the variance of all the genes in class j.
In mixture modelling, it is convenient to formulate the model using a missing data problem, where each observation w_{ i }is assumed to arise from a specific but unknown component z_{ i }of the mixture. The model can be written in terms of the missing data z = z_{1}, ..., z_{ g }, otherwise known as allocation variables, which are assumed to be independent realizations of discrete random variables Z_{1}, ..., Z_{ g }with
Pr(Z_{ i }= jπ, ψ, ν) = π_{ j }(i = 1, ..., g; j = 1, ..., k).
The allocation of gene i to component j is based on the highest posterior probability π_{ ij }over all the components. Genes in the same latent class j share the same variance 1/ψ_{ j }, the mean variance of all the n_{ j }genes in component j, estimated from a larger number of replicates than purely those per gene (n_{1} + … + n_{ k }= n). Parameters of the model can be estimated by the EM algorithm as shown in [17].
Quantities of interest, such as the posterior classification probabilities in (12), are estimated by pluggingin the point estimates $\widehat{\pi}$ and $\widehat{\psi}$ of π and ψ, respectively. Such plugin estimates do not account for all the variability in estimating the model parameters and, as such, are more likely to underestimate the variance of the model parameters, which might inflate the significance levels. Aside from the problems associated with estimating the variability, the EM algorithm can sometimes have computational problems, not least in finding local maximum of the likelihood surface among several possible local maxima instead of the global maximum. In order to avoid the possibility of such problems, we propose a Bayesian hierarchical structure for the mixture model of the gene differential variance. We use an exchangeable gamma(ς, τ) prior on ψ_{ i }and a dirichlet(α_{1}, ..., α_{ k }) prior on π. The hyperparameters ς and τ can be influential, and therefore in our full Bayesian analysis, these are not fixed, but given vague hyperpriors.
Declarations
Acknowledgements
SOMM gratefully acknowledges the support of the Medical Research Council (MRC) (Grant Ref G0400908). REW was funded by the Engineering and Physical Sciences Research Council (EPSRC) on a postgraduate fellowship and MSG was funded by the Higher Education Funding Council of England (HEFCE). We thank two anonymous referees for their very helpful and constructive comments.
Authors’ Affiliations
References
 Amaratunga D, Cabrera J: Exploration and Analysis of DNA Microarray and Protein Array Data. 2004, New Jersey: Wiley
 Brown PO, Bolstein D: Exploring the new world of the genome with microarrays. Nature Genetics. 1999, 21 (suppl 1): 3337. 10.1038/4462.View ArticlePubMed
 Efron B, Tibshirani R, Storey JD, Tusher V: Empirical Bayes analysis of a microarray experiment. Journal of the American Statistical Association. 2001, 96: 11521160. 10.1198/016214501753382129.View Article
 Lonnstedt I, Speed T: Replicated microarray data. Statist Sin. 2002, 12: 203218.
 Schena M, Heller RA, Theriault TP, Konrad K, Lachenmeier E, Davis RW: Microarrays: Biotechnology's Discovery Platform for Functional Genomics. Trends in Biotechnology. 1998, 16: 301306. 10.1016/S01677799(98)012190.View ArticlePubMed
 Li C, Tseng GC, Wong WH: Modelbased analysis of oligonucleotide arrays and issues in CDNA microarray analysis. Statistical Analysis of Gene Expression Microarray Data. Edited by: Speed T. 2003, Boca Raton, Florida: Chapman & Hall, 134.
 Yang YH, Speed T: Design and analysis of comparative microarray experiments. Statistical Analysis of Gene Expression Microarray Data. Edited by: Speed T. 2003, Boca Raton, Florida: Chapman & Hall, 3591.
 Tusher V, Tibshirani R, Gilbert C: Significance analysis of microarrays applied to the ionizing radiation response. Proceedings of the National Academy of Sciences, USA. 2001, 98: 51165121. 10.1073/pnas.091062498.View Article
 Dean N, Raftery AE: Normal uniform mixture differential gene expression detection for cDNA microarrays. BMC Bioinformatics. 2005, 6: 17310.1186/147121056173.PubMed CentralView ArticlePubMed
 Efron B, Tibshirani R: Empirical Bayes methods and false discovery rates for microarrays. Genetic Epidemiology. 2002, 23: 7086. 10.1002/gepi.1124.View ArticlePubMed
 Broet P, Richardson S, Radvanyi F: Bayesian hierachical model for identifying changes in gene expression from microarray experiments. Journal of Computational Biology. 2002, 9: 671683. 10.1089/106652702760277381.View ArticlePubMed
 Richardson S, Green PJ: On Bayesian analysis of mixtures with an unknown number of components (with discussion). Journal of the Royal Statistical Society, Ser. B. 1997, 59: 731792. 10.1111/14679868.00095.View Article
 Newton MA, Kendziorski CM, Richmond CS, Blattner FR, Tsui KW: On differential variability of expression ratios: Improving statistical inference about gene expression changes from microarray data. Journal of Computational Biology. 2001, 8: 3752. 10.1089/106652701300099074.View ArticlePubMed
 Pan W, Lin J, le Delmar CT: A mixture model approach to detecting differentially expressed genes with microarray data. Fund Integr Genomics. 2003, 3: 1171402. 10.1007/s1014200300857.View Article
 Baldi P, Long AD: A Bayesian framework for the analysis of microarray expression data: regularised ttest and statistical inferences of gene expression changes. Bioinformatics. 2001, 17: 509519. 10.1093/bioinformatics/17.6.509.View ArticlePubMed
 Lewin A, Richardson S, Marshall C, Glazier A, Aitman T: Bayesian modelling of differential gene expression. Biometrics. 2005, 62: 1018. 10.1111/j.15410420.2005.00394.x.View Article
 Delmar P, Robin S, Roux DT, Daudin JJ: Mixture model on the variance for differential analysis of gene expression data. Journal of the Royal Statistical Society, C. 2005, 54: 3150. 10.1111/j.14679876.2005.00468.x.View Article
 Ji Y, Wu C, Liu P, Wang J, Coombes KR: Applications of betamixture models in bioinformatics. Bioinformatics. 2005, 21:
 Spiegelhalter D, Thomas A, Best N, Lunn D: BUGS: Bayesian Inference Using Gibbs Sampling, Version 1.4. 2003, MRC Biostatistics Unit: Cambridge
 Alizadeh AA, Eisen MB, Davis RE, Ma C, Lossos IS, Rosenwald A, Boldrick JC, Sabet H, Tran T, Yu X, Powell JI, Yang L, Marti GE, Moore T, Hudson J, Lu L, Lewis DB, Tibshirani R, Sherlock G, Chan WC, Greiner TC, Weisenburger DD, Armitage JO, Warnke R, Levy R, Wilson W, Grever MR, Byrd JC, Botstein D, Brown PO, Staudt LM: Distinct types of diffuse large Bcell lymphoma identified by gene expression profilling. Nature. 2000, 403: 503511. 10.1038/35000501.View ArticlePubMed
 OliveiraBrochado A, Martins FV: Assessing the number of components in mixture models: a review. ResearchWork In Progress – n. 194. Faculty de Economia. University of Porto
 Brame R, Nagin DS, Wasserman L: Exploring some analytical characterictics of finite mixture modes. Journal of Quantitative Criminology. 2006, 22: 3159. 10.1007/s1094000590018. 2118–2122View Article
 Huber W, von Heydebreck A, Sultmann H, Poustka A, Vingron M: Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics. 2002, 18: S96S104.View ArticlePubMed
 Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JYH, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biology. 2004, 5 (10): R8010.1186/gb2004510r80. Sep 15PubMed CentralView ArticlePubMed
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.