SegAnnotBasesC: SegAnnotBasesC
In SegAnnot: Annotation-aware segmentation via dynamic programing

Description Usage Arguments Author(s) Examples

C implementation of mapping bases to indices.

1	SegAnnotBasesC(y, base, starts, ends)

`y`	The vector of observations to segment.
`base`	Position of each observation in base pairs.
`starts`	Starts of the 1-annotated regions in base pairs.
`ends`	Ends of the 1-annotated regions in base pairs.

Guillem Rigaill, Toby Dylan Hocking

library(SegAnnot)
data(profiles,package="SegAnnot")
if(interactive()){
  for(info in profiles){
    ann <- subset(info$ann, annotation=="1breakpoint")
    pro <- info$pro
    fit <- SegAnnotBases(pro$log, pro$pos, ann$min, ann$max)
    fitc <- SegAnnotBasesC(pro$log, pro$pos, ann$min, ann$max)
    for(must.be.equal in c("sR","eR","change","cost")){
      stopifnot(all(fit[[must.be.equal]] == fitc[[must.be.equal]]))
    }
  }
}