getPolymerizationRate: Computes the Polymerization Rate for each Alignment

Description Usage Arguments Details Value Examples

Description

'getPolymerizationRate' computes the polymerization rate for each alignment in the cmpH5 file.

Usage

1
getPolymerizationRate(cmpH5, idx = seq.int(1, nrow(cmpH5)))

Arguments

cmpH5

An object of class PacBioCmpH5.

idx

The indices of the alignments to return.

Details

Each alignment in the dataset has one polymerization rate value.

Value

Returns a vector of polymerization rates.

Examples

1
2
3
4
5
6
7
8
  ## Not run: 
  require(pbh5)
  cmpH5 <- PacBioCmpH5(system.file("h5_files",  "aligned_reads.cmp.h5",
  package = "pbh5"))
  pRate <- getPolymerizationRate(cmpH5)
  plot(density(pRate, na.rm = TRUE), xlab = "Polymerization Rate")
 
## End(Not run)

PacificBiosciences/R-pbh5 documentation built on May 7, 2019, 11:54 p.m.