normalizeBySum: Normalization of NucleotideDistr by global number of reads

Description Usage Arguments Value Author(s) See Also Examples

View source: R/NDtransforms.R

Description

normalizeBySum function normalizes the coverage values in NucleotideDistr by dividing all the numbers for all samples by the sum of reads for each sample. The number of reads from each sample may be taken from the database by the function getSumsExp, which is a wrapper for an appropriate SQL procedure. Alternatively, it is passed directly as a vector of numeric values of the same length as the number of samples analyzed. Such simple normalization allows comparisons of the coverage values for samples with different number of reads

Usage

1

Arguments

nd

NucleotideDistr object with raw read counts

r

Vector of numbers. If there is no such parameter, a database procedure summarizing reads is run

Value

NucleotideDistr object

Author(s)

Michal Okoniewski, Anna Lesniewska

See Also

getSumsExp

Examples

1
2
3
4
5
6
7
#  if (xmapConnected())  
#  {
#    rs <- newSeqReads(1,10000,20000,1)
#    nd.cov <- getCoverageFromRS(rs,1:3)
#    nd.norm <- normalizeBySum(nd.cov)
#    nd.norm <- normalizeBySum(nd.cov, r=c(100, 200, 1000))
#  }

rnaSeqMap documentation built on Nov. 8, 2020, 5:50 p.m.