ratio.plot: Plot the estimated log ratios against the overall intensities

In rama: Robust Analysis of MicroArrays

Description

Plot the estimated log2(gamma1/gamma2) against log2(gamma1*gamma2)/2.

Usage

ratio.plot(mcmc.obj,col=1,pch=1)

Arguments

mcmc.obj	An object of class mcmc as returned by fit.model
col	The color to be used for the symbols
pch	The type of symbols to be used.

Value

The graph!

Author(s)

Raphael Gottardo

References

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 98195-4322

See Also

fit.model

Examples

data(hiv)
### Initialize the proposals
mcmc.hiv<-fit.model(hiv[1:10,c(1:4)],hiv[1:10,c(5:8)],B=2000,min.iter=000,batch=1,shift=30,mcmc.obj=NULL,dye.swap=TRUE,nb.col1=2)
ratio.plot(mcmc.hiv,col=1,pch=1)

Example output

rama documentation built on Nov. 8, 2020, 8:02 p.m.