diffCountsBaySeq: Differential Counts wrapper - BaySeq

Description Usage Arguments Value Author(s) References See Also Examples

View source: R/diffCountsBaySeq.R

Description

Most Rcade users will not need to call this function directly. A function that provides a wrapper for the methods in the BaySeq package.

Usage

1
diffCountsBaySeq(counts, targets, annoZones, cl = NULL, getLibsizesArgs = list(estimationType = "quantile", quantile = 0.75), getPriors.NBArgs = list(), getLikelihoods.NBArgs = list(), libsizes)

Arguments

counts

Counts from countReads

targets

Data.frame - Information about the ChIP data files. Mandatory column names are: "fileid", "sampleid", "factor", "filepath".

annoZones

GRanges specifying the bins of interest, with a column in the metadata for the geneID.

cl

cluster from makeCluster in the parallel package.

getLibsizesArgs

List - Arguments to be passed to the getLibsizes function. If a libsizes column is present in the targets file, then these arguments are ignored.

getLibsizesArgs$cD is always ignored.

See getLibsizes for a list of arguments.

getPriors.NBArgs

See getPriors for a list of arguments.

getPriors.NBArgs$cD and getPriors.NBArgs$cl are always ignored.

getLikelihoods.NBArgs

See getLikelihoods for a list of arguments.

getLikelihoods.NBArgs$cD and getLikelihoods.NBArgs$cl are always ignored.

libsizes

Library sizes FIXME

Value

data.frame containing differential count information.

Author(s)

Jonathan Cairns

References

Hardcastle, T. J., & Kelly, K. A. (2010). baySeq: Empirical Bayesian methods for identifying differential expression in sequence count data. BMC Bioinformatics, 11, 422.

See Also

RcadeAnalysis

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
	dir <- file.path(system.file("extdata", package="Rcade"), "STAT1")

	targets <- read.csv(file.path(dir, "targets.csv"), as.is = TRUE)

	anno <- read.csv(file.path(dir, "anno.csv"))
	anno <- anno[order(anno$chromosome_name),]
	colnames(anno) <- c("ENSG","chr","start","end","str")
	ChIPannoZones <- defineBins(anno, zone=c(-1500, 1500), geneID="ENSG")

	counts <- countReads(ChIPannoZones, targets, fileDir = dir)

	x <- diffCountsBaySeq(counts, targets, ChIPannoZones)

jmc200/Rcade documentation built on June 17, 2020, 5:59 a.m.