diagplotAvgFtd {metaseqR2} | R Documentation

## Create average False (or True) Discovery curves

### Description

This function creates false (or true) discovery curves using a list containing several outputs from `diagplotFtd`.

### Usage

        diagplotAvgFtd(ftdrObj, output = "x11",
            path = NULL, draw = TRUE, ...)

### Arguments

`ftdrObj` | a list with outputs from `diagplotFtd`.
`output` | one or more R plotting device to direct the plot result to. Supported mechanisms: `"x11"` (default), `"png"`, `"jpg"`, `"bmp"`, `"pdf"` or `"ps"`.
`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 `TRUE` (make plots).
`...` | further arguments to be passed to plot devices, such as parameter from `par`.

### Value

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.

### Author(s)

Panagiotis Moulos

### Examples

    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)

[Package metaseqR2 version 1.0.11 Index]