# A Bayesian Approach for Decision Making on the Identification of Genes with Different Expression Levels: An Application toEscherichia coli Bacterium Data

Authors: Erlandson F. Saraiva; Francisco Louzada; Luís A. Milan; Silvana Meira; Juliana Cobre
Journal: Computational and Mathematical Methods in Medicine
 (2012)

---

## Abstract

A common interest in gene expression data analysis is to identify from a large pool of candidate genes the genes that present significant changes in expression levels between a treatment and a control biological condition. Usually, it is done using a statistic value and a cutoff value that are used to separate the genes differentially and nondifferentially expressed. In this paper, we propose a Bayesian approach to identify genes differentially expressed calculating sequentially credibility intervals from predictive densities which are constructed using the sampled mean treatment effect from all genes in study excluding the treatment effect of genes previously identified with statistical evidence for difference. We compare our Bayesian approach with the standard ones based on the use of thet-test and modified t-tests via a simulation study, using small sample sizes which are common in gene expression data analysis. Results obtained report evidence that the proposed approach performs better than standard ones, especially for cases with mean differences and increases in treatment variance in relation to control variance. We also apply the methodologies to a well-known publicly available data set on Escherichia coli bacterium.


---

## Body

## 1. Introduction

The DNA array technology is capable of providing gene expression levels measurements for thousands of genes simultaneously under different biological experimental conditions. In these experiments, total RNA is reverse-transcribed to create either radioactive or fluorescent-labeled cDNA which is hybridized with a large DNA library of gene fragments attached to a glass or membrane support [1]. After this, a scanner of high resolution is used to obtain the color intensity from each spot. So, the color intensities are normalized in order to obtain the expression level of genes. For further discussion and additional references on DNA array technology see [2–9].Obtaining the expression levels, a common objective is to identify genes that present significant changes in gene expression levels between treatment and control experimental condition. The identification of these genes is important because it may bring to light new biological discoveries, such as which genes may be involved in the origin and/or evolution of the same disease of genetic origin or which genes react to a drug stimulus. Thus, we aim to establish the use of these experiments as tools in medicine [10].As the observed expression levels incorporate different sources of variability present in the process of obtaining fluorescent intensity measurements [11], statistical methods are important to identify the genes differentially expressed. One of the first approaches proposed to identify genes differentially expressed was the fold-change approach [2, 3]. In this approach, a gene is considered differentially expressed if the average of the logarithm of the observed expression levels in treatment and control varies more than a cutoff point, Rc, which is previously prefixed. This approach however is not adequate to yield good results, once a cutoff value Rc may have different significance for different observed expression levels. Besides, this approach does not consider the variability of the observed expression levels from treatment and control.Another method commonly used for gene expression data analysis is the so-called two-samplet-test (TT) for the log transformed data [1, 8]. The problem with the application of TT to this kind of data is the usual small size of treatment and control samples in genetic studies, which may lead to underestimated variances and small power of test. To avoid such limitations, some TT modifications were proposed, such as the Cyber-t (CT) proposed by [1] and the Bayesian t-test (BTT) proposed by [12]. Basically, the main idea is to consider modifications of the standard error estimate of the two-sample difference present in the denominator of the standard t-statistics.In this paper, we propose a Bayesian approach to identify genes differentially expressed by calculating sequentially credibility intervals from predictive densities which are constructed using all treatment effects excluding the treatment effect of genes previously identified with statistical evidence for difference. This procedure avoids the small sample size problem, usual in gene expression data analysis, and allows us to use the normality assumption for observed data [9].In order to verify the performance of the proposed approach and compare it with the conventional ones based on the use of thet-test and modified t-tests, we present a simulation study. The comparison is done in terms of the true positive rate, false positive rate, and true discovery rate. Results obtained report evidence that our proposed approach performs better than t-test and modified t-tests, especially for cases with mean differences and increases in treatment variance in relation to control variance. We also apply the methods to a real dataset, obtained from the experiment carried through with Escherichia coli bacterium, described in [5].The paper is organized as follows. In Section2, we develop our Bayesian approach constructing the predictive density and describing our criteria to identify the genes differentially expressed. In Section 3, the method is compared with the t-test and modified t-tests using simulated datasets and a real dataset. In Section 4, we conclude the paper with final remarks on the proposed method.


## 2. Predictive Modeling for Gene Expression

Consider a DNA array experiment withn genes and two experimental conditions which we name by control (c) and treatment (t). Suppose that control and treatment are replicated nc and nt times, respectively. Denote by xigh the ith observed expression level (or its logarithm) for gene g in experimental condition h,h∈{c,t}, and g=1,…,n. Let xgh={x1gh,…,xnhgh} be realizations of independent random variables Xgh={X1gh,…,Xnhgh}, for g=1,…,n and h∈{c,t}.Consider that(1)Yg=1nt∑i=1ntXigt-1nc∑i=1ncXigc
is the sampled mean treatment effect for gene g,g=1,…,n, and Y={Y1,…,Yn} is the set of all sampled mean treatment effects.Thus, consideringY, we can determine the predictive density for a new observation Yn+1, given Y, and build a 100(1-α)% credibility interval for Yn+1. In order to develop our idea and as often found in gene expression data analysis [1, 8, 9, 11–13], we assume that Y is an independent sample generated from a normal distribution with mean μ and variance σ2,(2)Y1,Y2,…,Yn~N(μ,σ2).The likelihood function is given by(3)L(μ,σ2)∝(σ2)-n/2exp⁡{-n(y¯-μ)2+(n-1)s22σ2},
where y¯=(1/n)∑g=1nyg and s2=(1/(n-1))∑g=1n(yg-y¯)2 are the sample mean and variance of y, respectively.Since parametersμ and σ2 have a direct interpretation in the context of the gene expression data analysis, so we may express expert opinions in terms of prior distributions for parameters. In order to explore the fully conjugation, consider that joint prior distribution for parameters μ and σ2 is given by(4)μ∣σ2~N(μ0,σ2λ),σ2~IG(τ2,β2),
where μ0,λ,τ, and β are known hyperparameters, and ℐ𝒢(·) represent the inverse gamma distribution with mean (β/2)/(τ/2)-1.Updating the prior distributions in (4) via likelihood function in (3), the joint posterior distribution for (μ,σ2) is given by(5)μ,σ2∣y~N(μ*,σ2λ+n)IG(τ*2,β*2),
where μ*=(n/(n+λ)y¯)+(λ/(n+λ))μ0, τ*=τ+n+1, and β*=β+(n-1)·s2+nλ(y¯-μ0)2/(n+λ), for g=1,…,n.Considering now thatYn+1 is a new observation, independent from Y, the predictive distribution of Yn+1∣y is given by(6)Yn+1∣y~tτ*(μ*,β*(n+λ+1)τ*(n+λ)),
where tτ* represents the Student's t-distribution with location parameter μ*, scale β*(n+λ+1)/τ*(n+λ), and τ* degrees of freedom.From (6), the variance of Yn+1∣y is given by(7)Var⁡(Yn+1∣y)=(τ*τ*-2)⋅(β*(n+λ+1)τ*(n+λ)),
and a 100(1-α)% credibility interval for Yn+1∣y is given by(8)I(1-α)(Yn+1∣y)=(μ*-t1-α/2,τ*Var⁡(Yn+1∣y),μ*+t1-α/2,τ*Var⁡(Yn+1∣y)),
where t1-α/2,τ* denotes the quantile 1-α/2 of the standard t-student distribution with τ* degrees of freedom.


### 2.1. Predictive Approach Criterion

Letyord={y(1),y(2),…,y(n)}, the set y in increasing numerical order, y(1)<y(2)<⋯<y(n). Assuming that y(g) is a future observation in relation to the set composite by all observed treatment effect except the gth treatment effect, y(-g)ord={y(1),…,y(g-1),y(g+1),…,y(n)}, so the distribution of Y(g) is given by (6) with posterior parameter calculated using y(-g)ord and I1-α(Y(g)∣y(-g)ord) is a 100(1-α)% credibility interval for Y(g), g=1,2,…,n.In order to identify the genes differentially expressed, we fixE(μ∣y(-g)ord)=μ*=0 and set up I(g)threshold=t1-α/2,τ*Var({y(-g)ord}∖{ydif}) as a threshold, where {ydif} is the set that will be composite by treatment effect from genes identified with evidences for difference, and Var({y(-g)ord}∖{ydif}) is calculated according to (7) for the set {y(-g)ord} excluding the set {ydif}. For g=1,…,n, the identification of the genes differentially expressed is given by the following steps:  (i)
CalculateI(g)threshold;(ii)
if|y(g)|≤I(g)threshold, then gene g does not presents statistical evidence for differences;(iii)
if|y(g)|>I(g)threshold, then gene g present statistical evidence for differences. Do ydif=ydif∪y(g).


## 2.1. Predictive Approach Criterion

Letyord={y(1),y(2),…,y(n)}, the set y in increasing numerical order, y(1)<y(2)<⋯<y(n). Assuming that y(g) is a future observation in relation to the set composite by all observed treatment effect except the gth treatment effect, y(-g)ord={y(1),…,y(g-1),y(g+1),…,y(n)}, so the distribution of Y(g) is given by (6) with posterior parameter calculated using y(-g)ord and I1-α(Y(g)∣y(-g)ord) is a 100(1-α)% credibility interval for Y(g), g=1,2,…,n.In order to identify the genes differentially expressed, we fixE(μ∣y(-g)ord)=μ*=0 and set up I(g)threshold=t1-α/2,τ*Var({y(-g)ord}∖{ydif}) as a threshold, where {ydif} is the set that will be composite by treatment effect from genes identified with evidences for difference, and Var({y(-g)ord}∖{ydif}) is calculated according to (7) for the set {y(-g)ord} excluding the set {ydif}. For g=1,…,n, the identification of the genes differentially expressed is given by the following steps:  (i)
CalculateI(g)threshold;(ii)
if|y(g)|≤I(g)threshold, then gene g does not presents statistical evidence for differences;(iii)
if|y(g)|>I(g)threshold, then gene g present statistical evidence for differences. Do ydif=ydif∪y(g).


## 3. Data Analysis

In this section, we illustrate the predictive approach (PA) applied to artificial and real datasets. The real data set was extracted from the site (http://www.jbc.org/) and refers to an experiment realized with Escherichia Coli bacterium using nylon membranes, described by [5].Moreover, we compare the PA results with the results obtained by considering three well-known methods to identify differentially expressed genes: the two-samplet-test (TT) and Cyber-t test (CT) proposed by [1] and with the Bayesian t-test (BTT) proposed by [12].In the TT, the hypothesis test is based on the statistics(9)tg=x¯gt-x¯gcsgt2/nt+sgc2/nc,
which follows a Student's t-distribution with df=[sgc2/nc+sgt2/nt]2/[(sgc2/nc)2/(nc-1)+(sgt2/nt)2/(nt-1)], degrees of freedom, where x¯gh and sgh2 are the sample mean and variance for gene g in experimental condition h={c,t}. Fixing a significance level α, if |tg| is greater than a threshold t1-α/2,df (quantile 1-α/2 of Student's t distribution with df degrees of freedom), then the test conclude for difference of expression levels.Reference [1] proposed a two-sample t-test replacing the denominator of (1) by a pooled variance estimated via a Bayesian approach. So, the authors implement the Cyber-t software using the statistics(10)tg=x¯gt-x¯gcσ̃gt2/ngt+σ̃gc2/ngc
and the degrees of freedom df=ν0+ngc+ngt-2, where σ̃gh2=ν0σ02+(nh-1)sgh2/(ν0+nh-2), for h∈{c,t}, where ν0 and σ02 are hyperparameters. The authors assume that k>2 points are needed to properly estimate the standard deviation and keep ng+ν0=k, where ng=ngc+ngt. They suggest to fix k=10 and so ν0=10-ng. To fix a value for σ02, the authors say “one could use the standard deviation of the entire set of observations or, depending on the situation, of particularly categories of genes.” Using only the information from observations of the gene g, we fix σ02=((ng-1)/ng)sg2, as suggested by [1], where sg2 is the sample variance of the set xg={xgc,xgt}.Based on [1, 12], develop a Bayesian approach and show that(11)Δμ-Δx¯σn1/ngt+1/ngc|xgc,xgt~tνn,
where Δx¯=x¯gt-x¯gc,νnσn2=ν0σ02+(ngc-1)sgc2+(ngt-1)sgt2, νn=ν0+ngc+ngt-2, and tνn represent the Student's t distribution with νn degrees of freedom, for g=1,…,n. As suggested by authors, we fix ν0=ng and σ02=sg2.


### 3.1. Artificial Data

To generate the artificial data sets, we consider that observations from control group are generated from a normal distribution with meanμc and variance σc2, Xigc~𝒩(μc,σc2), for i=1,…,nc and g=1,…,n. We fix μgc=-14 and σgc2=0.8. These values are the average of the observed mean and variance of the expression levels (log transformed) from control group of the Escherichia coli bacterium dataset. We fix n=1.000, and the sample sizes nc and nt were fixed at 4 and 8.To generate the observations from treatment group, we follow the steps:(i)
From index{1,…,n}, we choose randomly p% of these index to indicate the cases generated with difference, p∈{5,10,20};(ii)
if the indexg∈{1,…,n} is chosen, then we consider an indicator variable 𝕀g=1 and generate Xigt~𝒩(μt,σt2), for i=1,…,nt. In order to verify how the method behaves when (μt,σt2) moves away from (μc,σc2), we simulate its values using(12)μgt=μgc±δσc,σt=γσc,
for δ={0,0.25,0.50,0.75,1,1.25,1.50,1.75,2} and γ={1,2,3}. The differential cases represent the situation overexpressed, defined by signal + in expression μt, and under expressed, defined by the signal − in expression of μt. We use p=pover+punder, for pover={3,7,5} and punder={2,3,15}, respectively. For example, p=5 is composite by pover=3 plus punder=2(iii)
if the indexg∈{1,…,n} is not chosen, then set up 𝕀g=0 and generate Xigt~𝒩(μc,σc2), for i=1,…,nc.For PA application, we fix the hyperparameters in order to have weak informative priors. We set up (i)τ and β in a way that E[σ2]=(β/2)/((τ/2)-1)=R2, where R=max(y)-min(y) is the length of the interval of variation of the observed data y. Thus, we obtain β=(τ-2)·R2. So, we fix τ=3, μ0=0, and λ=10-2. We also set up α=0.05 to calculate the credibility intervals and for the t-tests.To record the cases identified with difference by PA, we consider an indicator variable𝕀gPA=1 for cases, so that |y(g)|>I(g)threshold. Otherwise, 𝕀gPA=0. Analogously, for TT, CT, and BTT, we consider 𝕀gmethod=1 (method={TT,CT,BTT}) for cases with Pvalueg<0.05 and 𝕀gmethod=0, otherwise.In order to compare the performance of methods, we calculate the true positive rate given by(13)Pmethod=∑g=1nIg⋅Igmethod∑g=1nIg,
where method={PA,TT,CT,BTT}.We calculate the true positive rate forS=100 different datasets generated according to steps (i) to (iii) described above, and we present the results using the mean of the true positive rate, that is given by P¯method=∑s=1SPmethod(s)/S, where Pmethod(s) is the true positive rate calculated for sth generated dataset for method={PA,TT,CT,BTT}.Tables1, 2, and 3 present the P¯method value for nc=nt=4 and Tables 4, 5, and 6 present the P¯method value for nc=nt=8, for method={PA,TT,CT,BTT}.Table 1
True positive rate,nc=nt=4, p=5% (3% over and 2% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0670.0720.1200.1890.3080.4100.5320.6530.738TT0.0400.0460.0720.1220.1980.2900.3860.5070.601CT0.0310.0370.0560.0920.1600.2430.3230.4420.533BTT0.0410.0460.0730.1220.2000.2950.3880.5070.6042PA0.2200.2290.2580.3080.3630.4250.4920.5600.630TT0.0500.0490.0630.0890.1070.1360.1820.2210.277CT0.0380.0420.0500.0730.0920.1180.1550.1990.250BTT0.0530.0540.0640.0920.1140.1450.1940.2390.2993PA0.3570.3700.3670.3970.4240.4540.4860.5300.584TT0.0490.0570.0590.0670.0810.1020.1120.1450.170CT0.0450.0510.0500.0580.0740.0940.1060.1380.162BTT0.0600.0630.0680.0770.0940.1160.1330.1640.195Table 2
True positive rate,nc=nt=4, p=10% (7% over and 3% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0630.0760.1100.1810.2570.3590.4470.5380.631TT0.0430.0540.0740.1210.1910.2860.3790.4930.613CT0.0320.0400.0540.0960.1540.2330.3210.4260.545BTT0.0430.0550.0740.1230.1930.2860.3840.4970.6162PA0.2080.2110.2370.2750.3160.3670.4340.4970.556TT0.0450.0520.0640.0850.1030.1410.1840.2330.286CT0.0360.0420.0500.0700.0860.1200.1560.2050.259BTT0.0480.0540.0670.0900.1080.1500.1950.2530.3123PA0.3170.3140.3370.3460.3650.4070.4340.4730.510TT0.0490.0580.0560.0710.0800.0990.1180.1440.173CT0.0420.0510.0490.0600.0720.0900.1100.1340.165BTT0.0560.0650.0650.0800.0910.1110.1360.1660.200Table 3
True positive rate,nc=nt=4, p=20% (5% over and 15% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0660.0770.1170.1760.2410.3220.4110.5100.610TT0.0430.0490.0790.1230.1970.2850.3840.4940.606CT0.0310.0350.0600.0980.1570.2320.3220.4300.540BTT0.0420.0490.0800.1240.1970.2870.3870.5010.6132PA0.1790.1870.2230.2600.3130.3690.4250.4870.545TT0.0490.0510.0630.0830.1030.1420.1840.2360.286CT0.0390.0400.0510.0690.0840.1220.1600.2080.254BTT0.0510.0550.0670.0880.1100.1530.1970.2530.3073PA0.2740.2760.3030.3280.3560.3830.4280.4640.511TT0.0520.0560.0630.0700.0800.0980.1220.1420.175CT0.0450.0490.0550.0630.0730.0900.1130.1350.167BTT0.0590.0640.0710.0790.0930.1130.1380.1650.201Table 4
True positive rate,nc=nt=8, p=5% (3% over and 2% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0610.0840.1790.3230.4890.6360.7870.8790.939TT0.0470.0750.1450.2820.4600.6300.7980.8990.957CT0.0380.0570.1200.2420.4150.5800.7650.8730.944BTT0.0470.0730.1440.2830.4590.6300.7980.8970.9582PA0.2180.2370.3040.3790.4770.5830.6630.7510.819TT0.0510.0590.0880.1360.2110.3010.4030.5100.621CT0.0430.0510.0790.1230.1880.2740.3770.4870.597BTT0.0530.0620.0940.1430.2190.3140.4220.5320.6403PA0.3520.3650.3870.4330.4810.5470.6070.6710.731TT0.0480.0590.0720.0910.1290.1740.2240.2900.361CT0.0430.0530.0660.0860.1240.1660.2180.2800.354BTT0.0550.0680.0770.1040.1450.1920.2470.3180.393Table 5
True positive rate,nc=nt=8, p=10% (7% over and 3% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0670.0880.1680.2870.4070.5360.6640.7730.858TT0.0520.0710.1460.2820.4430.6260.7850.8940.957CT0.0410.0560.1250.2480.3990.5810.7510.8700.945BTT0.0510.0700.1460.2810.4430.6250.7850.8940.9572PA0.2040.2240.2670.3320.4160.4940.5790.6480.728TT0.0500.0590.0910.1350.2140.3040.4020.5080.628CT0.0420.0500.0800.1220.1930.2800.3780.4850.604BTT0.0520.0630.0950.1430.2230.3160.4190.5250.6473PA0.3320.3250.3400.3830.4180.4780.5310.5800.634TT0.0500.0560.0710.0930.1270.1710.2230.2850.356CT0.0460.0520.0660.0880.1210.1630.2150.2770.348BTT0.0570.0630.0790.1030.1420.1930.2460.3130.387Table 6
True positive rate,nc=nt=8, p=20% (5 over and 15 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0660.0920.1620.2630.3810.5130.6580.7680.864TT0.0520.0730.1500.2780.4540.6320.7900.8950.958CT0.0400.0600.1250.2430.4070.5890.7540.8720.945BTT0.0510.0720.1490.2770.4530.6320.7900.8960.9572PA0.1770.2060.2530.3180.4080.4990.5750.6520.725TT0.0510.0610.0880.1390.2120.3030.4040.5070.627CT0.0420.0540.0780.1240.1940.2810.3800.4820.605BTT0.0540.0640.0930.1450.2220.3170.4220.5280.6473PA0.2710.2940.3220.3700.4170.4680.5260.5870.630TT0.0520.0570.0670.0950.1250.1750.2270.2820.350CT0.0490.0530.0620.0910.1190.1690.2200.2750.343BTT0.0600.0660.0770.1080.1400.1940.2520.3100.383As we move from the left to the right side of the tables, in each line, we have the distances between control and treatment means, which are increasing. As we move from top to down in columns of the tables, we have the distance between the treatment and control variances, which are increasing. Increasing the sample sizes from 4 to 8, all methods increase its performance.Fornc=4 and all values of δ and γ used, the PA present better performance than TT, CT, and BTT. Moving away the treatment distribution from the control distribution (increasing δ and γ), the true positive rate obtained by PA is greater than t-tests.Fornc=8 and γ=1 fixed the TT, CT, and BTT present better performance than PA for the same values of δ used. For example, for p=5 and δ={1.5,1.75,2} the t tests present greater true positive rate than PA. But increasing the value of γ, γ={2,3}, the PA presents greater true positive rate.Besides, note that TT, CT, and BTT present similar results with a slight advantage for BTT, that is, greater true positive rate than TT and CT. Also note that true positive rate obtained by CT is smaller than TT for all cases simulated. It happens because we use only the information from observations from geneg to fix the hyperparameter σ02. In order to obtain better results, [1, 12] suggest to fix σ0 as the standard deviation estimated by pooling together all the neighboring genes contained in a window of size w. But, the authors do not discuss how to define a good value w to lead to satisfactory results.We also compare the performance of the methods using the mean of the false positive rate and the mean of the true discovery rate. The mean of the false positive rates is presented in Tables7–9 and in Tables 10–12. All methods present a small false positive rate.Table 7
False positive rate,nc=nt=4, p=5% (3 over and 2 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0620.0620.0580.0540.0500.0450.0390.0330.028TT0.0400.0410.0420.0410.0400.0420.0420.0410.040CT0.0300.0300.0310.0290.0290.0300.0310.0300.029BTT0.0410.0410.0420.0410.0400.0420.0420.0410.0402PA0.0510.0510.0500.0450.0410.0370.0310.0270.023TT0.0410.0420.0420.0410.0410.0410.0410.0400.041CT0.0300.0310.0300.0300.0290.0300.0300.0290.030BTT0.0410.0420.0420.0410.0410.0420.0420.0400.0413PA0.0400.0370.0380.0330.0300.0250.0230.0180.016TT0.0400.0410.0400.0410.0430.0400.0410.0420.040CT0.0290.0300.0290.0300.0310.0290.0300.0310.030BTT0.0400.0410.0400.0410.0430.0400.0410.0420.041Table 8
False positive rate,nc=nt=4, p=10% (7 over and 3 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0610.0600.0560.0480.0390.0310.0230.0170.012TT0.0400.0400.0410.0420.0420.0410.0410.0400.041CT0.0290.0300.0300.0300.0300.0300.0300.0300.030BTT0.0400.0410.0410.0420.0410.0410.0410.0400.0412PA0.0460.0440.0390.0330.0270.0210.0170.0120.009TT0.0410.0400.0420.0420.0400.0410.0410.0410.041CT0.0300.0290.0310.0300.0290.0300.0300.0300.030BTT0.0410.0400.0410.0410.0400.0410.0410.0410.0413PA0.0290.0260.0230.0200.0160.0130.0100.0070.005TT0.0410.0420.0400.0420.0410.0410.0400.0410.041CT0.0300.0300.0290.0300.0290.0310.0290.0300.030BTT0.0410.0420.0400.0420.0410.0420.0410.0410.042Table 9
False positive rate,nc=nt=4, p=20% (5 over and 15 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0620.0620.0550.0460.0360.0270.0200.0140.011TT0.0400.0400.0410.0410.0410.0410.0420.0420.041CT0.0300.0290.0290.0290.0300.0300.0300.0300.030BTT0.0400.0410.0400.0410.0410.0410.0410.0420.0422PA0.0340.0350.0340.0290.0270.0200.0160.0110.008TT0.0410.0410.0420.0410.0410.0400.0400.0420.041CT0.0300.0300.0310.0290.0290.0290.0290.0300.029BTT0.0410.0410.0420.0410.0410.0400.0400.0420.0413PA0.0160.0160.0150.0150.0140.0110.0100.0070.005TT0.0420.0410.0420.0400.0410.0410.0420.0400.042CT0.0300.0290.0300.0290.0300.0300.0310.0290.031BTT0.0420.0410.0410.0400.0410.0410.0420.0400.042Table 10
False positive rate,nc=nt=8, p=5% (3% over and 2% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0630.0600.0560.0480.0420.0340.0260.0190.014TT0.0470.0480.0470.0470.0470.0470.0500.0480.047CT0.0370.0370.0360.0370.0370.0370.0390.0380.037BTT0.0470.0470.0460.0470.0470.0470.0500.0480.0472PA0.0530.0510.0460.0410.0320.0260.0200.0150.011TT0.0470.0470.0480.0500.0480.0480.0480.0470.048CT0.0370.0370.0360.0370.0370.0370.0390.0380.037BTT0.0470.0470.0480.0490.0470.0470.0480.0470.0483PA0.0400.0390.0340.0290.0230.0200.0140.0110.008TT0.0470.0480.0480.0490.0480.0480.0480.0490.048CT0.0370.0370.0380.0390.0380.0380.0370.0380.038BTT0.0470.0470.0480.0490.0470.0480.0480.0490.048Table 11
False positive rate,nc=nt=8, p=10% (7% over and 3% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0640.0590.0490.0390.0270.0170.0110.0060.003TT0.0490.0480.0480.0490.0480.0480.0490.0470.048CT0.0390.0370.0380.0380.0380.0380.0380.0370.037BTT0.0490.0480.0480.0490.0480.0480.0480.0470.0482PA0.0450.0410.0340.0260.0190.0120.0070.0040.002TT0.0480.0470.0470.0470.0470.0490.0480.0480.046CT0.0370.0370.0360.0370.0370.0390.0370.0370.037BTT0.0480.0470.0470.0470.0470.0490.0470.0480.0463PA0.0270.0240.0200.0150.0100.0070.0040.0020.001TT0.0480.0480.0480.0460.0470.0500.0470.0480.048CT0.0370.0380.0370.0370.0360.0390.0360.0370.038BTT0.0480.0480.0480.0470.0460.0500.0470.0480.048Table 12
False positive rate,nc=nt=8, p=20% (5 over and 15 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0630.0580.0470.0320.0210.0140.0110.0060.004TT0.0480.0480.0480.0480.0490.0480.0480.0490.049CT0.0380.0370.0370.0380.0380.0380.0380.0380.038BTT0.0480.0470.0480.0480.0480.0480.0480.0490.0482PA0.0360.0360.0300.0230.0170.0120.0070.0040.002TT0.0470.0480.0480.0480.0480.0480.0480.0470.048CT0.0370.0380.0370.0370.0370.0370.0370.0370.038BTT0.0470.0480.0480.0480.0470.0480.0470.0470.0483PA0.0140.0170.0160.0120.0100.0070.0050.0020.001TT0.0480.0480.0470.0480.0470.0480.0470.0470.048CT0.0370.0380.0360.0380.0370.0380.0370.0360.038BTT0.0480.0470.0470.0480.0470.0480.0470.0470.048The mean of the true discovery rates is presented in Tables13–15 and in Tables 16–18. The PA presents greater true discovery rate than t-tests for all values of δ and γ used. Besides, note that increasing the value of δ and γ, the true discovery rate increases in both directions for PA. But the same does not happen with the t-tests, in which the proportion of identification increases only as the value of δ increases, that is, when the mean of the treatment distribution moves away from the mean of the control distribution. Increasing the variance of treatment (increasing the value of γ), the t-tests present a reduction of its performance, in opposite to PA which presents an improvement in its performance.Table 13
True discovery rate,nc=nt=4, p=5% (3 over and 2 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0530.0570.0980.1550.2460.3260.4200.5130.586TT0.0500.0560.0830.1370.2080.2690.3290.3950.444CT0.0530.0600.0890.1440.2240.3000.3590.4390.497BTT0.0520.0550.0840.1370.2080.2720.3310.3940.4422PA0.1850.1920.2170.2670.3210.3800.4630.5260.596TT0.0590.0580.0740.1020.1220.1480.1880.2240.262CT0.0620.0680.0790.1150.1450.1710.2140.2680.304BTT0.0630.0630.0750.1060.1290.1540.1970.2400.2773PA0.3250.3510.3410.3940.4300.4930.5330.6100.662TT0.0610.0690.0710.0790.0900.1180.1260.1520.182CT0.0740.0820.0830.0910.1110.1460.1570.1910.223BTT0.0730.0750.0830.0900.1040.1310.1460.1700.202Table 14
True discovery rate,nc=nt=4, p=10% (7 over and 3 under).
γMethodδ00.250.500.7511.251.51.7521PA0.1020.1240.1800.2950.4240.5640.6820.7760.857TT0.1060.1310.1680.2420.3380.4380.5080.5770.623CT0.1070.1300.1700.2610.3650.4670.5470.6160.673BTT0.1060.1320.1680.2460.3430.4390.5100.5790.6272PA0.3370.3500.4060.4830.5690.6610.7470.8220.877TT0.1090.1270.1460.1850.2240.2760.3350.3890.436CT0.1150.1380.1560.2070.2490.3040.3700.4360.490BTT0.1130.1320.1540.1940.2310.2870.3470.4080.4583PA0.5550.5750.6190.6660.7220.7850.8360.8820.922TT0.1170.1340.1340.1590.1790.2090.2450.2810.318CT0.1350.1570.1550.1810.2140.2440.2950.3300.379BTT0.1310.1470.1530.1760.2000.2280.2700.3090.347Table 15
True discovery rate,nc=nt=4, p=20% (5 over and 15 under).
γMethodδ00.250.500.7511.251.51.7521PA0.2090.2350.3500.4880.6290.7490.8410.9030.934TT0.2090.2320.3280.4320.5440.6340.6990.7470.786CT0.2030.2310.3360.4570.5680.6610.7300.7830.817BTT0.2070.2300.3310.4340.5460.6360.7040.7500.7862PA0.5720.5750.6260.6920.7430.8220.8740.9190.944TT0.2290.2410.2750.3400.3850.4700.5360.5840.638CT0.2430.2560.2890.3670.4140.5140.5830.6320.684BTT0.2350.2530.2880.3500.4020.4910.5520.6040.6543PA0.8150.8170.8340.8450.8700.9010.9180.9470.961TT0.2370.2560.2750.3030.3300.3730.4230.4670.511CT0.2690.2970.3180.3520.3810.4310.4800.5330.577BTT0.2570.2840.3010.3310.3650.4090.4520.5060.547Table 16
True discovery rate,nc=nt=8, p=5% (3% over and 2% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0480.0680.1450.2620.3850.5030.6190.7170.790TT0.0490.0760.1410.2390.3410.4140.4580.4990.517CT0.0520.0750.1500.2580.3710.4520.5120.5490.573BTT0.0500.0750.1410.2420.3400.4150.4590.4980.5192PA0.1790.2000.2610.3310.4430.5460.6380.7270.800TT0.0540.0610.0870.1270.1890.2490.3090.3640.407CT0.0580.0670.0980.1440.2120.2780.3540.4120.457BTT0.0560.0640.0930.1320.1970.2590.3190.3750.4163PA0.3240.3360.3810.4500.5290.6000.6990.7620.825TT0.0510.0600.0730.0890.1240.1620.1980.2380.286CT0.0580.0690.0840.1040.1480.1890.2360.2780.332BTT0.0580.0700.0780.1010.1380.1760.2150.2560.304Table 17
True discovery rate,nc=nt=8, p=10% (7% over and 3% under).
γMethodδ00.250.500.7511.251.51.7521PA0.1040.1410.2760.4520.6290.7770.8680.9330.970TT0.1050.1400.2540.3920.5070.5950.6430.6780.690CT0.1060.1440.2690.4200.5400.6340.6880.7240.740BTT0.1030.1400.2540.3920.5080.5950.6440.6800.6912PA0.3360.3820.4670.5910.7160.8210.9010.9460.973TT0.1040.1220.1770.2420.3340.4090.4850.5420.603CT0.1120.1310.1970.2680.3660.4480.5320.5940.649BTT0.1080.1290.1850.2540.3440.4190.4970.5520.6113PA0.5820.6030.6570.7430.8180.8820.9290.9630.984TT0.1040.1120.1410.1820.2320.2760.3460.3980.451CT0.1210.1310.1640.2100.2700.3170.3950.4560.507BTT0.1170.1260.1550.1980.2530.3020.3690.4230.473Table 18
True discovery rate,nc=nt=8, p=20% (5 over and 15 under).
γMethodδ00.250.500.7511.251.51.7521PA0.2070.2830.4650.6710.8200.9030.9390.9690.983TT0.2130.2790.4380.5910.7010.7660.8050.8220.832CT0.2110.2880.4560.6160.7280.7930.8340.8510.863BTT0.2120.2780.4380.5920.7030.7660.8060.8220.8322PA0.5560.5910.6790.7760.8560.9100.9530.9750.988TT0.2110.2420.3130.4200.5270.6130.6790.7290.767CT0.2210.2640.3410.4540.5700.6530.7190.7650.801BTT0.2230.2530.3260.4320.5400.6230.6910.7390.7733PA0.8310.8170.8400.8870.9150.9420.9670.9840.992TT0.2130.2300.2660.3340.3970.4760.5470.6010.645CT0.2470.2580.3020.3770.4470.5270.5990.6550.693BTT0.2370.2550.2920.3610.4260.5020.5740.6250.667Results show a better performance of the PA in relation to TT, CT, and BTT, specially, when the difference refers to variance of the variable involved. From the biological practical point of view, it shows us that PA may identify with differences genes which are not identified by TT, CT, and BTT, specially, genes with differences in means and variances.


### 3.2. Escherichia coli Data

In this section consider the gene expression data set onEscherichia coli bacterium, composed by n=4290 genes [5]. Figure 1 shows the treatment and control observed means and variances for all genes of this dataset.Treatment and control observed means and variances.
(a)
Observed means(b)
Observed variancesResults for PA are presented in Figure2. Results for TT, CT, and BTT are presented in Figure 3. These figures show the observed treatment and control means and variances of genes identified with evidence for difference by considering PA, TT, CT, and BTT, respectively. The PA identifies 340 genes with evidences for difference, while TT identifies 222, CT 219, and BTT 288 genes.Treatment and control observed means and variances and genes identified with evidence for difference by PA.
(a)
Genes identified by PA(b)
Genes identified by PATreatment and control observed means and variances and genes identified with evidence for difference by TT, CT, and BTT.
(a)
Genes identified by TT(b)
Genes identified by TT(c)
Genes identified by CT(d)
Genes identified by CT(e)
Genes identified by BTT(f)
Genes identified by BTTNote that genes with means well apart are better identified by PA than by the other methods. Moreover, genes with mean and variances well apart are identified by PA and not identified by TT, CT, and BTT, as can be noted in Figure2. Examples are genes 2766 (b1326(f262)) and 3254 (dbpA) that are highlighted in Figures 2(a) and 2(b). Genes with means well apart and similar variances are however identified by TT, CT, and BTT. An example is the gene 10 (hdeB) that is highlighted in Figures 2(a) and 2(b). One possible reason for this is the low performance of TT, CT and BTT in situations with differences in means and variances, as observed in the artificial data sets. Besides, note that PA is capable to identify differentially expressed genes which are not identified by TT, CT, and BTT, specially, genes with differences in means and variances.


## 3.1. Artificial Data

To generate the artificial data sets, we consider that observations from control group are generated from a normal distribution with meanμc and variance σc2, Xigc~𝒩(μc,σc2), for i=1,…,nc and g=1,…,n. We fix μgc=-14 and σgc2=0.8. These values are the average of the observed mean and variance of the expression levels (log transformed) from control group of the Escherichia coli bacterium dataset. We fix n=1.000, and the sample sizes nc and nt were fixed at 4 and 8.To generate the observations from treatment group, we follow the steps:(i)
From index{1,…,n}, we choose randomly p% of these index to indicate the cases generated with difference, p∈{5,10,20};(ii)
if the indexg∈{1,…,n} is chosen, then we consider an indicator variable 𝕀g=1 and generate Xigt~𝒩(μt,σt2), for i=1,…,nt. In order to verify how the method behaves when (μt,σt2) moves away from (μc,σc2), we simulate its values using(12)μgt=μgc±δσc,σt=γσc,
for δ={0,0.25,0.50,0.75,1,1.25,1.50,1.75,2} and γ={1,2,3}. The differential cases represent the situation overexpressed, defined by signal + in expression μt, and under expressed, defined by the signal − in expression of μt. We use p=pover+punder, for pover={3,7,5} and punder={2,3,15}, respectively. For example, p=5 is composite by pover=3 plus punder=2(iii)
if the indexg∈{1,…,n} is not chosen, then set up 𝕀g=0 and generate Xigt~𝒩(μc,σc2), for i=1,…,nc.For PA application, we fix the hyperparameters in order to have weak informative priors. We set up (i)τ and β in a way that E[σ2]=(β/2)/((τ/2)-1)=R2, where R=max(y)-min(y) is the length of the interval of variation of the observed data y. Thus, we obtain β=(τ-2)·R2. So, we fix τ=3, μ0=0, and λ=10-2. We also set up α=0.05 to calculate the credibility intervals and for the t-tests.To record the cases identified with difference by PA, we consider an indicator variable𝕀gPA=1 for cases, so that |y(g)|>I(g)threshold. Otherwise, 𝕀gPA=0. Analogously, for TT, CT, and BTT, we consider 𝕀gmethod=1 (method={TT,CT,BTT}) for cases with Pvalueg<0.05 and 𝕀gmethod=0, otherwise.In order to compare the performance of methods, we calculate the true positive rate given by(13)Pmethod=∑g=1nIg⋅Igmethod∑g=1nIg,
where method={PA,TT,CT,BTT}.We calculate the true positive rate forS=100 different datasets generated according to steps (i) to (iii) described above, and we present the results using the mean of the true positive rate, that is given by P¯method=∑s=1SPmethod(s)/S, where Pmethod(s) is the true positive rate calculated for sth generated dataset for method={PA,TT,CT,BTT}.Tables1, 2, and 3 present the P¯method value for nc=nt=4 and Tables 4, 5, and 6 present the P¯method value for nc=nt=8, for method={PA,TT,CT,BTT}.Table 1
True positive rate,nc=nt=4, p=5% (3% over and 2% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0670.0720.1200.1890.3080.4100.5320.6530.738TT0.0400.0460.0720.1220.1980.2900.3860.5070.601CT0.0310.0370.0560.0920.1600.2430.3230.4420.533BTT0.0410.0460.0730.1220.2000.2950.3880.5070.6042PA0.2200.2290.2580.3080.3630.4250.4920.5600.630TT0.0500.0490.0630.0890.1070.1360.1820.2210.277CT0.0380.0420.0500.0730.0920.1180.1550.1990.250BTT0.0530.0540.0640.0920.1140.1450.1940.2390.2993PA0.3570.3700.3670.3970.4240.4540.4860.5300.584TT0.0490.0570.0590.0670.0810.1020.1120.1450.170CT0.0450.0510.0500.0580.0740.0940.1060.1380.162BTT0.0600.0630.0680.0770.0940.1160.1330.1640.195Table 2
True positive rate,nc=nt=4, p=10% (7% over and 3% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0630.0760.1100.1810.2570.3590.4470.5380.631TT0.0430.0540.0740.1210.1910.2860.3790.4930.613CT0.0320.0400.0540.0960.1540.2330.3210.4260.545BTT0.0430.0550.0740.1230.1930.2860.3840.4970.6162PA0.2080.2110.2370.2750.3160.3670.4340.4970.556TT0.0450.0520.0640.0850.1030.1410.1840.2330.286CT0.0360.0420.0500.0700.0860.1200.1560.2050.259BTT0.0480.0540.0670.0900.1080.1500.1950.2530.3123PA0.3170.3140.3370.3460.3650.4070.4340.4730.510TT0.0490.0580.0560.0710.0800.0990.1180.1440.173CT0.0420.0510.0490.0600.0720.0900.1100.1340.165BTT0.0560.0650.0650.0800.0910.1110.1360.1660.200Table 3
True positive rate,nc=nt=4, p=20% (5% over and 15% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0660.0770.1170.1760.2410.3220.4110.5100.610TT0.0430.0490.0790.1230.1970.2850.3840.4940.606CT0.0310.0350.0600.0980.1570.2320.3220.4300.540BTT0.0420.0490.0800.1240.1970.2870.3870.5010.6132PA0.1790.1870.2230.2600.3130.3690.4250.4870.545TT0.0490.0510.0630.0830.1030.1420.1840.2360.286CT0.0390.0400.0510.0690.0840.1220.1600.2080.254BTT0.0510.0550.0670.0880.1100.1530.1970.2530.3073PA0.2740.2760.3030.3280.3560.3830.4280.4640.511TT0.0520.0560.0630.0700.0800.0980.1220.1420.175CT0.0450.0490.0550.0630.0730.0900.1130.1350.167BTT0.0590.0640.0710.0790.0930.1130.1380.1650.201Table 4
True positive rate,nc=nt=8, p=5% (3% over and 2% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0610.0840.1790.3230.4890.6360.7870.8790.939TT0.0470.0750.1450.2820.4600.6300.7980.8990.957CT0.0380.0570.1200.2420.4150.5800.7650.8730.944BTT0.0470.0730.1440.2830.4590.6300.7980.8970.9582PA0.2180.2370.3040.3790.4770.5830.6630.7510.819TT0.0510.0590.0880.1360.2110.3010.4030.5100.621CT0.0430.0510.0790.1230.1880.2740.3770.4870.597BTT0.0530.0620.0940.1430.2190.3140.4220.5320.6403PA0.3520.3650.3870.4330.4810.5470.6070.6710.731TT0.0480.0590.0720.0910.1290.1740.2240.2900.361CT0.0430.0530.0660.0860.1240.1660.2180.2800.354BTT0.0550.0680.0770.1040.1450.1920.2470.3180.393Table 5
True positive rate,nc=nt=8, p=10% (7% over and 3% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0670.0880.1680.2870.4070.5360.6640.7730.858TT0.0520.0710.1460.2820.4430.6260.7850.8940.957CT0.0410.0560.1250.2480.3990.5810.7510.8700.945BTT0.0510.0700.1460.2810.4430.6250.7850.8940.9572PA0.2040.2240.2670.3320.4160.4940.5790.6480.728TT0.0500.0590.0910.1350.2140.3040.4020.5080.628CT0.0420.0500.0800.1220.1930.2800.3780.4850.604BTT0.0520.0630.0950.1430.2230.3160.4190.5250.6473PA0.3320.3250.3400.3830.4180.4780.5310.5800.634TT0.0500.0560.0710.0930.1270.1710.2230.2850.356CT0.0460.0520.0660.0880.1210.1630.2150.2770.348BTT0.0570.0630.0790.1030.1420.1930.2460.3130.387Table 6
True positive rate,nc=nt=8, p=20% (5 over and 15 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0660.0920.1620.2630.3810.5130.6580.7680.864TT0.0520.0730.1500.2780.4540.6320.7900.8950.958CT0.0400.0600.1250.2430.4070.5890.7540.8720.945BTT0.0510.0720.1490.2770.4530.6320.7900.8960.9572PA0.1770.2060.2530.3180.4080.4990.5750.6520.725TT0.0510.0610.0880.1390.2120.3030.4040.5070.627CT0.0420.0540.0780.1240.1940.2810.3800.4820.605BTT0.0540.0640.0930.1450.2220.3170.4220.5280.6473PA0.2710.2940.3220.3700.4170.4680.5260.5870.630TT0.0520.0570.0670.0950.1250.1750.2270.2820.350CT0.0490.0530.0620.0910.1190.1690.2200.2750.343BTT0.0600.0660.0770.1080.1400.1940.2520.3100.383As we move from the left to the right side of the tables, in each line, we have the distances between control and treatment means, which are increasing. As we move from top to down in columns of the tables, we have the distance between the treatment and control variances, which are increasing. Increasing the sample sizes from 4 to 8, all methods increase its performance.Fornc=4 and all values of δ and γ used, the PA present better performance than TT, CT, and BTT. Moving away the treatment distribution from the control distribution (increasing δ and γ), the true positive rate obtained by PA is greater than t-tests.Fornc=8 and γ=1 fixed the TT, CT, and BTT present better performance than PA for the same values of δ used. For example, for p=5 and δ={1.5,1.75,2} the t tests present greater true positive rate than PA. But increasing the value of γ, γ={2,3}, the PA presents greater true positive rate.Besides, note that TT, CT, and BTT present similar results with a slight advantage for BTT, that is, greater true positive rate than TT and CT. Also note that true positive rate obtained by CT is smaller than TT for all cases simulated. It happens because we use only the information from observations from geneg to fix the hyperparameter σ02. In order to obtain better results, [1, 12] suggest to fix σ0 as the standard deviation estimated by pooling together all the neighboring genes contained in a window of size w. But, the authors do not discuss how to define a good value w to lead to satisfactory results.We also compare the performance of the methods using the mean of the false positive rate and the mean of the true discovery rate. The mean of the false positive rates is presented in Tables7–9 and in Tables 10–12. All methods present a small false positive rate.Table 7
False positive rate,nc=nt=4, p=5% (3 over and 2 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0620.0620.0580.0540.0500.0450.0390.0330.028TT0.0400.0410.0420.0410.0400.0420.0420.0410.040CT0.0300.0300.0310.0290.0290.0300.0310.0300.029BTT0.0410.0410.0420.0410.0400.0420.0420.0410.0402PA0.0510.0510.0500.0450.0410.0370.0310.0270.023TT0.0410.0420.0420.0410.0410.0410.0410.0400.041CT0.0300.0310.0300.0300.0290.0300.0300.0290.030BTT0.0410.0420.0420.0410.0410.0420.0420.0400.0413PA0.0400.0370.0380.0330.0300.0250.0230.0180.016TT0.0400.0410.0400.0410.0430.0400.0410.0420.040CT0.0290.0300.0290.0300.0310.0290.0300.0310.030BTT0.0400.0410.0400.0410.0430.0400.0410.0420.041Table 8
False positive rate,nc=nt=4, p=10% (7 over and 3 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0610.0600.0560.0480.0390.0310.0230.0170.012TT0.0400.0400.0410.0420.0420.0410.0410.0400.041CT0.0290.0300.0300.0300.0300.0300.0300.0300.030BTT0.0400.0410.0410.0420.0410.0410.0410.0400.0412PA0.0460.0440.0390.0330.0270.0210.0170.0120.009TT0.0410.0400.0420.0420.0400.0410.0410.0410.041CT0.0300.0290.0310.0300.0290.0300.0300.0300.030BTT0.0410.0400.0410.0410.0400.0410.0410.0410.0413PA0.0290.0260.0230.0200.0160.0130.0100.0070.005TT0.0410.0420.0400.0420.0410.0410.0400.0410.041CT0.0300.0300.0290.0300.0290.0310.0290.0300.030BTT0.0410.0420.0400.0420.0410.0420.0410.0410.042Table 9
False positive rate,nc=nt=4, p=20% (5 over and 15 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0620.0620.0550.0460.0360.0270.0200.0140.011TT0.0400.0400.0410.0410.0410.0410.0420.0420.041CT0.0300.0290.0290.0290.0300.0300.0300.0300.030BTT0.0400.0410.0400.0410.0410.0410.0410.0420.0422PA0.0340.0350.0340.0290.0270.0200.0160.0110.008TT0.0410.0410.0420.0410.0410.0400.0400.0420.041CT0.0300.0300.0310.0290.0290.0290.0290.0300.029BTT0.0410.0410.0420.0410.0410.0400.0400.0420.0413PA0.0160.0160.0150.0150.0140.0110.0100.0070.005TT0.0420.0410.0420.0400.0410.0410.0420.0400.042CT0.0300.0290.0300.0290.0300.0300.0310.0290.031BTT0.0420.0410.0410.0400.0410.0410.0420.0400.042Table 10
False positive rate,nc=nt=8, p=5% (3% over and 2% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0630.0600.0560.0480.0420.0340.0260.0190.014TT0.0470.0480.0470.0470.0470.0470.0500.0480.047CT0.0370.0370.0360.0370.0370.0370.0390.0380.037BTT0.0470.0470.0460.0470.0470.0470.0500.0480.0472PA0.0530.0510.0460.0410.0320.0260.0200.0150.011TT0.0470.0470.0480.0500.0480.0480.0480.0470.048CT0.0370.0370.0360.0370.0370.0370.0390.0380.037BTT0.0470.0470.0480.0490.0470.0470.0480.0470.0483PA0.0400.0390.0340.0290.0230.0200.0140.0110.008TT0.0470.0480.0480.0490.0480.0480.0480.0490.048CT0.0370.0370.0380.0390.0380.0380.0370.0380.038BTT0.0470.0470.0480.0490.0470.0480.0480.0490.048Table 11
False positive rate,nc=nt=8, p=10% (7% over and 3% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0640.0590.0490.0390.0270.0170.0110.0060.003TT0.0490.0480.0480.0490.0480.0480.0490.0470.048CT0.0390.0370.0380.0380.0380.0380.0380.0370.037BTT0.0490.0480.0480.0490.0480.0480.0480.0470.0482PA0.0450.0410.0340.0260.0190.0120.0070.0040.002TT0.0480.0470.0470.0470.0470.0490.0480.0480.046CT0.0370.0370.0360.0370.0370.0390.0370.0370.037BTT0.0480.0470.0470.0470.0470.0490.0470.0480.0463PA0.0270.0240.0200.0150.0100.0070.0040.0020.001TT0.0480.0480.0480.0460.0470.0500.0470.0480.048CT0.0370.0380.0370.0370.0360.0390.0360.0370.038BTT0.0480.0480.0480.0470.0460.0500.0470.0480.048Table 12
False positive rate,nc=nt=8, p=20% (5 over and 15 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0630.0580.0470.0320.0210.0140.0110.0060.004TT0.0480.0480.0480.0480.0490.0480.0480.0490.049CT0.0380.0370.0370.0380.0380.0380.0380.0380.038BTT0.0480.0470.0480.0480.0480.0480.0480.0490.0482PA0.0360.0360.0300.0230.0170.0120.0070.0040.002TT0.0470.0480.0480.0480.0480.0480.0480.0470.048CT0.0370.0380.0370.0370.0370.0370.0370.0370.038BTT0.0470.0480.0480.0480.0470.0480.0470.0470.0483PA0.0140.0170.0160.0120.0100.0070.0050.0020.001TT0.0480.0480.0470.0480.0470.0480.0470.0470.048CT0.0370.0380.0360.0380.0370.0380.0370.0360.038BTT0.0480.0470.0470.0480.0470.0480.0470.0470.048The mean of the true discovery rates is presented in Tables13–15 and in Tables 16–18. The PA presents greater true discovery rate than t-tests for all values of δ and γ used. Besides, note that increasing the value of δ and γ, the true discovery rate increases in both directions for PA. But the same does not happen with the t-tests, in which the proportion of identification increases only as the value of δ increases, that is, when the mean of the treatment distribution moves away from the mean of the control distribution. Increasing the variance of treatment (increasing the value of γ), the t-tests present a reduction of its performance, in opposite to PA which presents an improvement in its performance.Table 13
True discovery rate,nc=nt=4, p=5% (3 over and 2 under).
γMethodδ00.250.500.7511.251.51.7521PA0.0530.0570.0980.1550.2460.3260.4200.5130.586TT0.0500.0560.0830.1370.2080.2690.3290.3950.444CT0.0530.0600.0890.1440.2240.3000.3590.4390.497BTT0.0520.0550.0840.1370.2080.2720.3310.3940.4422PA0.1850.1920.2170.2670.3210.3800.4630.5260.596TT0.0590.0580.0740.1020.1220.1480.1880.2240.262CT0.0620.0680.0790.1150.1450.1710.2140.2680.304BTT0.0630.0630.0750.1060.1290.1540.1970.2400.2773PA0.3250.3510.3410.3940.4300.4930.5330.6100.662TT0.0610.0690.0710.0790.0900.1180.1260.1520.182CT0.0740.0820.0830.0910.1110.1460.1570.1910.223BTT0.0730.0750.0830.0900.1040.1310.1460.1700.202Table 14
True discovery rate,nc=nt=4, p=10% (7 over and 3 under).
γMethodδ00.250.500.7511.251.51.7521PA0.1020.1240.1800.2950.4240.5640.6820.7760.857TT0.1060.1310.1680.2420.3380.4380.5080.5770.623CT0.1070.1300.1700.2610.3650.4670.5470.6160.673BTT0.1060.1320.1680.2460.3430.4390.5100.5790.6272PA0.3370.3500.4060.4830.5690.6610.7470.8220.877TT0.1090.1270.1460.1850.2240.2760.3350.3890.436CT0.1150.1380.1560.2070.2490.3040.3700.4360.490BTT0.1130.1320.1540.1940.2310.2870.3470.4080.4583PA0.5550.5750.6190.6660.7220.7850.8360.8820.922TT0.1170.1340.1340.1590.1790.2090.2450.2810.318CT0.1350.1570.1550.1810.2140.2440.2950.3300.379BTT0.1310.1470.1530.1760.2000.2280.2700.3090.347Table 15
True discovery rate,nc=nt=4, p=20% (5 over and 15 under).
γMethodδ00.250.500.7511.251.51.7521PA0.2090.2350.3500.4880.6290.7490.8410.9030.934TT0.2090.2320.3280.4320.5440.6340.6990.7470.786CT0.2030.2310.3360.4570.5680.6610.7300.7830.817BTT0.2070.2300.3310.4340.5460.6360.7040.7500.7862PA0.5720.5750.6260.6920.7430.8220.8740.9190.944TT0.2290.2410.2750.3400.3850.4700.5360.5840.638CT0.2430.2560.2890.3670.4140.5140.5830.6320.684BTT0.2350.2530.2880.3500.4020.4910.5520.6040.6543PA0.8150.8170.8340.8450.8700.9010.9180.9470.961TT0.2370.2560.2750.3030.3300.3730.4230.4670.511CT0.2690.2970.3180.3520.3810.4310.4800.5330.577BTT0.2570.2840.3010.3310.3650.4090.4520.5060.547Table 16
True discovery rate,nc=nt=8, p=5% (3% over and 2% under).
γMethodδ00.250.500.7511.251.51.7521PA0.0480.0680.1450.2620.3850.5030.6190.7170.790TT0.0490.0760.1410.2390.3410.4140.4580.4990.517CT0.0520.0750.1500.2580.3710.4520.5120.5490.573BTT0.0500.0750.1410.2420.3400.4150.4590.4980.5192PA0.1790.2000.2610.3310.4430.5460.6380.7270.800TT0.0540.0610.0870.1270.1890.2490.3090.3640.407CT0.0580.0670.0980.1440.2120.2780.3540.4120.457BTT0.0560.0640.0930.1320.1970.2590.3190.3750.4163PA0.3240.3360.3810.4500.5290.6000.6990.7620.825TT0.0510.0600.0730.0890.1240.1620.1980.2380.286CT0.0580.0690.0840.1040.1480.1890.2360.2780.332BTT0.0580.0700.0780.1010.1380.1760.2150.2560.304Table 17
True discovery rate,nc=nt=8, p=10% (7% over and 3% under).
γMethodδ00.250.500.7511.251.51.7521PA0.1040.1410.2760.4520.6290.7770.8680.9330.970TT0.1050.1400.2540.3920.5070.5950.6430.6780.690CT0.1060.1440.2690.4200.5400.6340.6880.7240.740BTT0.1030.1400.2540.3920.5080.5950.6440.6800.6912PA0.3360.3820.4670.5910.7160.8210.9010.9460.973TT0.1040.1220.1770.2420.3340.4090.4850.5420.603CT0.1120.1310.1970.2680.3660.4480.5320.5940.649BTT0.1080.1290.1850.2540.3440.4190.4970.5520.6113PA0.5820.6030.6570.7430.8180.8820.9290.9630.984TT0.1040.1120.1410.1820.2320.2760.3460.3980.451CT0.1210.1310.1640.2100.2700.3170.3950.4560.507BTT0.1170.1260.1550.1980.2530.3020.3690.4230.473Table 18
True discovery rate,nc=nt=8, p=20% (5 over and 15 under).
γMethodδ00.250.500.7511.251.51.7521PA0.2070.2830.4650.6710.8200.9030.9390.9690.983TT0.2130.2790.4380.5910.7010.7660.8050.8220.832CT0.2110.2880.4560.6160.7280.7930.8340.8510.863BTT0.2120.2780.4380.5920.7030.7660.8060.8220.8322PA0.5560.5910.6790.7760.8560.9100.9530.9750.988TT0.2110.2420.3130.4200.5270.6130.6790.7290.767CT0.2210.2640.3410.4540.5700.6530.7190.7650.801BTT0.2230.2530.3260.4320.5400.6230.6910.7390.7733PA0.8310.8170.8400.8870.9150.9420.9670.9840.992TT0.2130.2300.2660.3340.3970.4760.5470.6010.645CT0.2470.2580.3020.3770.4470.5270.5990.6550.693BTT0.2370.2550.2920.3610.4260.5020.5740.6250.667Results show a better performance of the PA in relation to TT, CT, and BTT, specially, when the difference refers to variance of the variable involved. From the biological practical point of view, it shows us that PA may identify with differences genes which are not identified by TT, CT, and BTT, specially, genes with differences in means and variances.


## 3.2. Escherichia coli Data

In this section consider the gene expression data set onEscherichia coli bacterium, composed by n=4290 genes [5]. Figure 1 shows the treatment and control observed means and variances for all genes of this dataset.Treatment and control observed means and variances.
(a)
Observed means(b)
Observed variancesResults for PA are presented in Figure2. Results for TT, CT, and BTT are presented in Figure 3. These figures show the observed treatment and control means and variances of genes identified with evidence for difference by considering PA, TT, CT, and BTT, respectively. The PA identifies 340 genes with evidences for difference, while TT identifies 222, CT 219, and BTT 288 genes.Treatment and control observed means and variances and genes identified with evidence for difference by PA.
(a)
Genes identified by PA(b)
Genes identified by PATreatment and control observed means and variances and genes identified with evidence for difference by TT, CT, and BTT.
(a)
Genes identified by TT(b)
Genes identified by TT(c)
Genes identified by CT(d)
Genes identified by CT(e)
Genes identified by BTT(f)
Genes identified by BTTNote that genes with means well apart are better identified by PA than by the other methods. Moreover, genes with mean and variances well apart are identified by PA and not identified by TT, CT, and BTT, as can be noted in Figure2. Examples are genes 2766 (b1326(f262)) and 3254 (dbpA) that are highlighted in Figures 2(a) and 2(b). Genes with means well apart and similar variances are however identified by TT, CT, and BTT. An example is the gene 10 (hdeB) that is highlighted in Figures 2(a) and 2(b). One possible reason for this is the low performance of TT, CT and BTT in situations with differences in means and variances, as observed in the artificial data sets. Besides, note that PA is capable to identify differentially expressed genes which are not identified by TT, CT, and BTT, specially, genes with differences in means and variances.


## 4. Discussion

Identifying genes with difference, in what concerns gene expression, may help biologists to study and understand some function of genes and infer possible relationships among genes and proteins. In this paper we propose a Bayesian approach to identify differentially expressed genes based on predictive density.In order to verify the performance of the PA approach and compare it with TT, CT, and BTT, we considered artificial and real datasets. Results show a better performance of PA in relation to thet-tests in identifying difference, mainly, in presence of different variances. The main advantage of the proposed method is that it is easy to use like a usual two-sample t-test but presents better performance in situations with small sample size.The biological interest in this fact is that PA may bring to light genes that are not identified when we use only the TT or the modifiedt-test ones. Moreover, the PA can be easily implemented in usual softwares such as the software R (the Comprehensive R Archive Network, http://cran.r-project.org/). The source code used for the data analysis was implemented in software R and can be obtained by emailing the authors.According to [1], gene expression data can be analyzed in at least three levels of increasing complexity. In the first level, each gene is analyzed separately, where the objective is to verify whether the observed expression in treatment experimental condition is significantly different from observed expression in control experimental condition. In the second level, clusters of genes are analyzed in terms of common functionalities and interactions. In the third level, the objective is to infer and understand the relationship among genes. As it should be clear by now, in this paper, we focus on the first level of analysis. However, as future work we intent to present the adjustment of the proposed method to control the false discovery rate for multiple testing hypotheses, when thousands of hypotheses are realized simultaneously, as proposed by [14], as well as to make a systematic comparison with his methodology. Besides, we also will development a multivariate approach in order to consider dependence among genes.



---

*Source: 953086-2012-03-05.xml*