Detect differentially expressed Transcription Start Sites between two conditions (test)

Source: R/AllGenerics.R, R/detectDiffTSS.R

Detect differentially expressed Transcription Start Sites between two conditions (test)

detectDiffTSS(fit, testGroup, contGroup, TSSfile, MAplot_fdr = NA)

# S4 method for DGEGLM
detectDiffTSS(fit, testGroup, contGroup, TSSfile,
  MAplot_fdr = NA)

Arguments

fit	DGEGLM object (output of fitDiffTSS command )
testGroup	Test group name
contGroup	Control group name
TSSfile	The TSS .bed file used for fitDiffTSS command
MAplot_fdr	FDR threshold to mark differentially expressed TSS in MAplot (NA = Don't make an MAplot)

Value

A GRanges object containing p-values of differential expression for each TSS.

Examples

# before running this
# 1. Create a CapSet object
# 2. de-multiplex the fastqs
# 3. map them
# 4. filter duplicate reads from mapped BAM
# 5. detect TSS
# 6. fit the diff TSS model.

# NOT RUN {
# load a previously saved DGEGLM object from step 5
csfit <- load("diffTSS_fit.Rdata")
dir <- system.file("extdata", package = "icetea")
# detect differentially expressed TSS between groups (return MA plot)
detectDiffTSS(csfit, testGroup = "mut", controlGroup = "wt",
               tssFile = file.path(dir, "testTSS_merged.bed"), MAplot_fdr = 0.05)

# }