getFCFromND {rnaSeqMap}R Documentation

getFCFromND - calculating fold change of coverages

Description

This function calculates the fold change of two sample coverages from a NucleotideDistr objects. The coverages are assumed to be after logarithmic transformation, so the function basically subtracts the value and generates new NucleotideDistr object with a single vector of fold changes.

Usage

getFCFromND(nd, exps)

Arguments

nd

NucleotideDistr object with coverages

exps

a pair of numbers of samples in the experiment

Value

NucleotideDistr object of type "FC" with a single vector of fold changes

Author(s)

Michal Okoniewski, Anna Lesniewska

Examples

#  if (xmapConnected())  
#  {
#    rs <- newSeqReads(1,1,20000,1)
#    rs <- addExperimentsToReadset(rs,1:3)
#    nd.cov <- getCoverageFromRS(rs,1:3)
#    nd.fc <- getFCFromND(nd.cov,c(1,3))
#  }

[Package rnaSeqMap version 2.48.0 Index]