R/ratio.plot.R

Defines functions ratio.plot

Documented in ratio.plot

ratio.plot<-function(mcmc.obj,col=1,pch=1)
  {
    if(class(mcmc.obj)!="mcmc")
      stop("'mcmc.obj' should be of type 'mcmc'" , call. = TRUE)
    
    
    ### Estimate the gene effects in sample 1
    x1<-mat.mean(mcmc.obj$gamma1)[,1]
    x2<-mat.mean(mcmc.obj$gamma2)[,1]

    plot((x1+x2)/2,x1-x2,pch=pch,col=col,xlab="Overall intensity (log2)",ylab="Log ratio (log2)")
    
  }

Try the rama package in your browser

Any scripts or data that you put into this service are public.

rama documentation built on May 6, 2019, 2:55 a.m.