diagplotAvgFtd {metaseqR2} | R Documentation |
This function creates false (or true) discovery
curves using a list containing several outputs
from diagplotFtd
.
diagplotAvgFtd(ftdrObj, output = "x11", path = NULL, draw = TRUE, ...)
ftdrObj |
a list with outputs from
|
output |
one or more R plotting
device to direct the plot result to.
Supported mechanisms: |
path |
the path to create output files. |
draw |
boolean to determine whether
to plot the curves or just return the
calculated values (in cases where the user
wants the output for later averaging
for example). Defaults to |
... |
further arguments to be passed to
plot devices, such as parameter from
|
A named list with two members: the first member
(avgFtdr
) contains a list with the
means and the standard deviations of the averaged
ftdrObj
and are used to create the plot.
The second member (path
) contains the
path to the created figure graphic.
Panagiotis Moulos
p11 <- 0.001*matrix(runif(300),100,3) p12 <- matrix(runif(300),100,3) p21 <- 0.001*matrix(runif(300),100,3) p22 <- matrix(runif(300),100,3) p31 <- 0.001*matrix(runif(300),100,3) p32 <- matrix(runif(300),100,3) p1 <- rbind(p11,p21) p2 <- rbind(p12,p22) p3 <- rbind(p31,p32) rownames(p1) <- rownames(p2) <- rownames(p3) <- paste("gene",1:200,sep="_") colnames(p1) <- colnames(p2) <- colnames(p3) <- paste("method",1:3,sep="_") truth <- c(rep(1,40),rep(-1,40),rep(0,20), rep(1,10),rep(2,10),rep(0,80)) names(truth) <- rownames(p1) ftdObj1 <- diagplotFtd(truth,p1,N=100,draw=FALSE) ftdObj2 <- diagplotFtd(truth,p2,N=100,draw=FALSE) ftdObj3 <- diagplotFtd(truth,p3,N=100,draw=FALSE) ftdObj <- list(ftdObj1,ftdObj2,ftdObj3) avgFtdObj <- diagplotAvgFtd(ftdObj)