runWF {signatureSearch} | R Documentation

## Run the entire GESS/FEA workflow

### Description

This function runs the entire GESS/FEA workflow when providing the query drug and cell type, as well as selecting the reference database (e.g. "cmap" or "lincs"), defining the specific GESS and FEA methods. In this case, the query GES is drawn from the reference database. The N (defined by the 'N_gess_drugs' argument) top ranking hits in the GESS tables were then used for FEA where three different annotation systems were used: GO Molecular Function (GO MF), GO Biological Process (GO BP) and KEGG pathways. The GESS/FEA results will be stored in a list object in R session. A working environment named by the use case will be created under users current working directory or under other directory defined by users. This environment contains a `results` folder where the GESS/FEA result tables were written to. The working environment also contains a template Rmd vignette as well as a rended HTML report, users could make modifications on the Rmd vignette as they need and re-render it to generate their HTML report.

### Usage

    runWF(
      drug,
      cell,
      refdb,
      gess_method,
      fea_method,
      N_gess_drugs = 100,
      env_dir = ".",
      tau = TRUE,
      Nup = 150,
      Ndown = 150,
      higher = 1,
      lower = -1,
      method = "spearman",
      pvalueCutoff = 1,
      qvalueCutoff = 1,
      minGSSize = 5,
      maxGSSize = 500
    )

### Arguments

`drug` | character(1) representing query drug name (e.g. vorinostat). This query drug should be included in the `refdb`
`cell` | character(1) indicating the cell type that the query drug treated in. Details about cell type options in LINCS database can be found in the `cell_info` table after load the 'signatureSearch' package and running 'data("cell_info")'
`refdb` | character(1), one of "lincs", "lincs_expr", "cmap", "cmap_expr", or path to the HDF5 file built from `build_custom_db` function
`gess_method` | character(1), one of "LINCS", "CORsub", "CORall", "Fisher", "CMAP", "gCMAP". When `gess_method` is "CORsub" or "CORall", only "lincs_expr" or "cmap_expr" databases are supported.
`fea_method` | character(1), one of "dup_hyperG", "mGSEA", "mabs", "hyperG", "GSEA"
`N_gess_drugs` | number of unique drugs in GESS result used as input of FEA
`env_dir` | character(1), directory under which the result environment located. The default is users current working directory in R session, can be checked via `getwd()` command in R
`tau` | TRUE or FALSE indicating whether to compute Tau scores if `gess_method` is set as 'LINCS'
`Nup` | integer(1). Number of most up-regulated genes to be subsetted for GESS query when `gess_method` is CMAP, LINCS or CORsub
`Ndown` | integer(1). Number of most down-regulated genes to be subsetted for GESS query when `gess_method` is CMAP, LINCS or CORsub
`higher` | numeric(1), it is defined when gess_method argument is 'gCMAP' or 'Fisher' representing the 'upper' threshold of subsetting genes with a score larger than 'higher'
`lower` | numeric(1), it is defined when gess_method argument is 'gCMAP' or 'Fisher' representing the 'lower' threshold of subsetting genes
`method` | One of 'spearman' (default), 'kendall', or 'pearson', indicating which correlation coefficient to use
`pvalueCutoff` | double, p-value cutoff for FEA result
`qvalueCutoff` | double, qvalue cutoff for FEA result
`minGSSize` | integer, minimum size of each gene set in annotation system
`maxGSSize` | integer, maximum size of each gene set in annotation system

### Value

list object containing GESS/FEA result tables

### Examples

    drug="vorinostat"; cell="SKB"
    refdb <- system.file("extdata", "sample_db.h5", package="signatureSearch")
    env_dir <- tempdir()
    wf_list <- runWF(drug, cell, refdb, gess_method="LINCS",
        fea_method="dup_hyperG", N_gess_drugs=10, env_dir=env_dir, tau=FALSE)

[Package signatureSearch version 1.1.5 Index]