ratio.plot {rama} | R Documentation |
Plot the estimated log2(gamma1/gamma2) against log2(gamma1*gamma2)/2.
ratio.plot(mcmc.obj,col=1,pch=1)
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. |
The graph!
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 98195-4322
fit.model
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)