getFC: Detect fold changes

Description Usage Arguments Value Author(s) Examples

View source: R/RNAsense.R

Description

For each gene and for each time point, RNA-seq count data is analyzed for fold changes between two experimental conditions. This functions bases on functions from the R package NBPSeq package for fold change analysis

Usage

1
2
getFC(dataset = mydata, myanalyzeConditions = analyzeConditions,
  cores = 1, mytimes = times, pValueFC = 0.01)

Arguments

dataset

Object of class SummarizedExperiment, output of SummarizedExperiment, as assays use a numeric matrix with your RNAseq count data, rows correspond to different genes, columns correspond to different experiments, as rowData provide a DataFrame with columns name (geneID) and genename (the gene names), as colData provide a DataFrame with columns condition, time and replicate

myanalyzeConditions

Character vector, Name of experimental conditions

cores

Numeric, Number of cores for parallelization, default 1 for no parallelization

mytimes

Numeric vector, Time points of the time-resolved RNA-seq data

pValueFC

Numeric, p-value for fold change detection

Value

Data.frame containing gene names, log fold change and p-values calculated from NBPSeq, each gene appears as often as available time points

Author(s)

Marcus Rosenblatt, marcus.rosenblatt@fdm.uni-freiburg.de

Examples

1
2
3
4
5
6
7
data(MZsox)
mydata <- MZsox[seq(1,nrow(MZsox), by=20),]
resultFC <- getFC(dataset = mydata,
myanalyzeConditions = c("WT", "MZsox"),
cores = 1,
mytimes = c(2.5,3,3.5,4,4.5,5,5.5,6),
pValueFC=0.01)

marcusrosenblatt/RNAsense documentation built on Dec. 7, 2020, 5:26 p.m.