DUreport: Differential gene expression and differential bin/junction...

Description Usage Arguments Value Author(s) See Also Examples

Description

Estimate differential expression at gene level and differential usage at bin and junction level.

Usage

1
DUreport(counts, targets, pair, group, minGenReads, minBinReads, minRds,ignoreExternal,threshold)

Arguments

counts

An object of class ASpliCounts

targets

A dataframe containing sample, bam and condition columns

pair

vector of length two, either numeric or character, providing the pair of groups to be compared

group

Factorial vector with tags for each sample

minGenReads

Default 10 reads

minBinReads

Default 5 reads

minRds

Default 0.05

ignoreExternal

Ignore Exon Bins at the beggining or end of the transcript. Default TRUE

threshold

Minimun number of junction. Default 5

Value

An ASpliDU object with results at genes, bins and junctions level

Author(s)

Estefania Mancini, Marcelo Yanovsky, Ariel Chernomoretz

See Also

DEXSeq, edgeR Accesors: genesDE, binsDU,junctionsDU Export: writeDU

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
library(RNAseqData.HNRNPC.bam.chr14)
chr14 <- system.file("extdata","chr14.sqlite", package="ASpli")
genome <- loadDb(chr14)
features <- binGenome(genome)
targets <- data.frame(bam=RNAseqData.HNRNPC.bam.chr14_BAMFILES,
                       condition=c(rep("CT",4),rep("KD",4)))
bam <- loadBAM(targets)
counts <- readCounts(features, bam, l=100L, maxISize=50000)
group <- factor(c(rep("CT",4), 
                rep("KD",4)))
pair <- c("CT","KD")  
du <- DUreport(counts, targets, pair, group)
writeDU(du, output.dir="only_du")


Search within the ASpli package
Search all R packages, documentation and source code

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.