getCumulativeAdvanceTime: Compute the CumulativeAdvanceTime
In PacificBiosciences/R-pbh5: Interface to Pacific Biosciences HDF5 Files

Description Usage Arguments Details Value See Also Examples

'getCumulativeAdvanceTime' computes the cumulative advance time for each idx. Cumulative advance time is cumsum(pw + ipd).

1	getCumulativeAdvanceTime(cmpH5, idx, unit = c("seconds", "frames"))

`cmpH5`	An object of class `PacBioCmpH5`.
`idx`	The indices of the alignments to retrieve.
`unit`	Whether you want seconds or frames.

Inter-pulse distances are given in seconds.

A list of vectors of inter-pulse distances.

PacBioCmpH5-class

  require(pbh5)
  cmpH5 <- PacBioCmpH5(system.file("h5_files",  "aligned_reads.cmp.h5",
                       package = "pbh5"))
  ipd <- getCumulativeAdvanceTime(cmpH5)
  logCumulativeAdvanceTime <- lapply(ipd, function(x) log(x))
  plot(density(logCumulativeAdvanceTime[[1]], na.rm = TRUE),
       main = "CumulativeAdvanceTimes")