getStats: Obtain statistics from a BSseqTstat object
In bsseq: Analyze, manage and store bisulfite sequencing data

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

Essentially an accessor function for the statistics of a BSseqTstat object.

1	getStats(bstat, regions = NULL, ...)

`bstat`	An object of class `BSseqStat` or `BSseqTstat`.
`regions`	An optional `data.frame` or `GenomicRanges` object specifying a number of genomic regions.
`...`	Additional arguments passed to the different backends based on the class of `bstat`; see Details.

Additional argument when the bstat object is of class BSseqTstat:

stat: Which statistics column should be obtained.

An object of class data.frame possible restricted to the regions specified.

Kasper Daniel Hansen khansen@jhsph.edu

BSseqTstat for the BSseqTstat class, and getCoverage and getMeth for similar functions, operating on objects of class BSseq.

if(require(bsseqData)) {
  data(BS.cancer.ex.tstat)
  head(getStats(BS.cancer.ex.tstat))
  reg <- GRanges(seqnames = c("chr22", "chr22"),
     ranges = IRanges(start = c(1, 2*10^7), end = c(2*10^7 +1, 4*10^7)))
  head(getStats(BS.cancer.ex.tstat, regions = reg))
}