Description Usage Arguments Value Author(s) See Also Examples
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
1 | normalizeBySum(nd,r=NULL)
|
nd |
|
r |
Vector of numbers. If there is no such parameter, a database procedure summarizing reads is run |
NucleotideDistr object
Michal Okoniewski, Anna Lesniewska
getSumsExp
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))
# }
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.