sigoraWrappR: sigora Wrapper function

Description Usage Arguments Value Examples

View source: R/sigoraWrappR.R

Description

Provides a wrapper for sigora::sigora() combined with sigora::ora() and returns a list used for the generation of an .Rmd report.

Usage

1
2
3
4
5
6
7
8
9
sigoraWrappR(
  data,
  threshold = 0.5,
  score_col = "",
  GPSrepos = sigora::load_data("kegH"),
  db = "",
  greater_than = TRUE,
  idmap = sigora::idmap
)

Arguments

data

input data.frame (at least two columns, first column containing IDs, other columns numerical ranks, i.e. fold changes)

threshold

fold change threshold above which (in absolute terms) a protein is considered differentially regulated

score_col

Name of the fold change column, in case the input file contains multiple contrasts

GPSrepos

GPS repository used as background, can be generated via sigoraWrappR

db

database used for the generation of the GPS repository

greater_than

Logical. Whether fc_threshold should be applied as greater_than (default is TRUE) or strictly less than (FALSE)

idmap

id mapping table, dafault sigora::idmap

Value

Returns a list containing the following elements:

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
library(prora)
df <- get_UniprotID_from_fasta_header(prora::exampleContrastData)
myGPSrepo <- makeGPS_wrappR(ids = df$UniprotID,target = "KEGG")
res <- sigoraWrappR(df,score_col = "estimate", GPSrepos = myGPSrepo$gps,
                     threshold = 0.5)

myGPSrepoGO <- makeGPS_wrappR(df$UniprotID, target = "GO")
res <- sigoraWrappR(df,score_col = "estimate", GPSrepos = myGPSrepoGO$gps,
                     threshold = 0.5)
## Not run: 
  res$ora
res$sigora$summary_results
res$data
GPSrepo <-  myGPSrepoGO$gps

## End(Not run)

protViz/fgczgseaora documentation built on Dec. 14, 2021, 9:22 p.m.