outputGeneTables: Output gene tables

Description Usage Arguments Value Author(s) Examples

View source: R/RNAsense.R

Description

Output information on switching genes (up/down) in tabular format (gene identifier/gene name) are created as .txt file and written to the specified working directory. Two of the generated files (geneNamelist) contain gene lists with gene name for genes that switch up and down respectively. The other two (genelist) contain exactly the same output but with gene identifiers instead of gene names depending on what you prefer for further analysis. Each column corresponds to a combination of switch time point, fold change direction and time point of fold change. All genes for which fold change was detected at the indicated time point and switch was detected at the indicated time point are listed in the corresponding column. Note that a single gene may appear multiple times. The fifth .txt file (switchList) contains information on detected switches in a different format. The output consists of table with six columns with each row corresponding to one gene. Detected switches are indicated by 1, -1 and 0 for switch up, switch down and no switch, respectively. If a switch was detected, the column timepoint indicated the corresponding time point of switch detection.

Usage

1
2
outputGeneTables(myresultCombined = resultCombined, mytimes = times,
  myanalyzeConditions = analyzeConditions, mywd = NULL)

Arguments

myresultCombined

data.frame, output of combineResults

mytimes

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

myanalyzeConditions

character vector, the conditions that were analyzed

mywd

character, working directory to which results will be written, if NULL the current working directory is used

Value

Working directory where results have been written to

Author(s)

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

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
library(ggplot2)
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))
resultSwitch <- getSwitch(dataset = mydata,
experimentStepDetection = "WT",
cores = 1,
mytimes = c(2.5,3,3.5,4,4.5,5,5.5,6))
resultCombined <- combineResults(resultSwitch, resultFC)
outputGeneTables(resultCombined,
mytimes = c(2.5,3,3.5,4,4.5,5,5.5,6),
myanalyzeConditions = c("WT", "MZsox"))

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