ratio.plot {rama}R Documentation

Plot the estimated log ratios against the overall intensities

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)

[Package rama version 1.0.1 Index]