Description Usage Arguments Details Value Author(s) References See Also Examples
Test for differentially expressed genes in a three sample setup. This code can be used with both cDNA microarrays or Affymetrix chip.
1 
sample1 
The matrix of intensity from the sample 1. Each row corresponds to a different gene. 
sample2 
The matrix of intensity from the sample 2. Each row corresponds to a different gene. 
sample3 
The matrix of intensity from the sample 3. Each row corresponds to a different gene. 
B 
The number of iteration used in the MCMC algorithm. 
min.iter 
The length of the burnin period in the MCMC algorithm. 
batch 
The thinning value to be used in the MCMC. Only every 
mcmc.obj 
An object of type bridge2, as returned by

all.out 
A logical value indicating if all the parameters
should be output. If 
verbose 
A logical value indicating if the current MCMC iteration number should be printed out. 
log 
A logical value indicating if the data are log transformed. 
robust 
A logical value indicating if a t model ( 
This code fits a robust Bayesian hierarchical model for testing for differential expression. Outliers are modeled explicitly using a $t$distribution. The model includes an exchangeable prior for the variances which allow different variances for the genes but still shrink extreme empirical variances. This function DO NOT perform normalization. The data should be normalized before hands such as centering the mean expression of each experiment. More details can be found in the references below.
An object of type bridge3
containing the sampled values from the
posterior distribution.
gamma1 
A matrix, each row contains the sampled values from the corresponding gene effect in sample 1. 
gamma2 
A matrix, each row contains the sampled values from the corresponding gene effect in sample 2. 
gamma3 
A matrix, each row contains the sampled values from the corresponding gene effect in sample 3. 
lambda.gamma1 
A vector containing the sampled values for the precision of the gene effect prior for the component corresponding to sample 1. 
lambda.gamma2 
A vector containing the sampled values for the precision of the gene effect prior for the component corresponding to sample 2. 
lambda.gamma3 
A vector containing the sampled values for the precision of the gene effect prior for the component corresponding to sample 3. 
lambda.gamma12 
A vector containing the sampled values for the precision of the gene effect prior for the component where sample 1 and sample 2 are combined. 
lambda.gamma13 
A vector containing the sampled values for the precision of the gene effect prior for the component where sample 1 and sample 3 are combined. 
lambda.gamma23 
A vector containing the sampled values for the precision of the gene effect prior for the component where sample 2 and sample 3 are combined. 
lambda.gamma123 
A vector containing the sampled values for the precision of the gene effect prior for the component where all the samples are combined. 
lambda_eps1 
A matrix, each row contains the sampled values from the corresponding gene precision in sample 1. 
lambda_eps2 
A matrix, each row contains the sampled values from the corresponding gene precision in sample 2. 
lambda_eps3 
A matrix, each row contains the sampled values from the corresponding gene precision in sample 3. 
a.eps1 
A vector containing the sampled values for the mean of the prior of the genes precision in sample 1. 
b.eps1 
A vector containing the sampled values for the variance of the prior of the genes precision in sample 1. 
a.eps2 
A vector containing the sampled values for the mean of the prior of the genes precision in sample 2. 
b.eps2 
A vector containing the sampled values for the variance of the prior of the genes precision in sample 2. 
a.eps3 
A vector containing the sampled values for the mean of the prior of the genes precision in sample 3. 
b.eps3 
A vector containing the sampled values for the variance of the prior of the genes precision in sample 3. 
w1 
A matrix, each element (i,j) correspond to the posterior mean of the sampled weights of replicate j in gene i and sample 1. To save memory, we only store the posterior means of the weigths. 
w2 
A matrix, each element (i,j) correspond to the posterior mean of the sampled weights of replicate j in gene i and sample 2. To save memory, we only store the posterior means of the weigths. 
w3 
A matrix, each element (i,j) correspond to the posterior mean of the sampled weights of replicate j in gene i and sample 3. To save memory, we only store the posterior means of the weigths. 
w1 
A matrix, each element (i,j) correspond to the posterior mean of the sampled weights of replicate j in gene i and sample 1. To save memory, we only store the posterior means of the weigths. 
nu1 
A matrix containing the sampled degrees of freedom in sample 1. 
nu2 
A matrix containing the sampled degrees of freedom in sample 2. 
nu3 
A matrix containing the sampled degrees of freedom in sample 3. 
w.mix 
The posterior mixing proportions in the mixture component. 
prop.model 
The posterior proportions of each component for each each gene. 
move 
The proportion of moves between components. This should be used as a diagnostic tool. 
Raphael Gottardo
Robust Estimation of cDNA Microarray Intensities with Replicates Raphael Gottardo, Adrian E. Raftery, Ka Yee Yeung, and Roger Bumgarner Department of Statistics, University of Washington, Box 354322, Seattle, WA 981954322
1 2 3 4 5 
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.