filterSegments {AneuFinder} | R Documentation |
filterSegments
filters out segments below a specified minimal segment size. This can be useful to get rid of boundary effects from the Hidden Markov approach.
filterSegments(segments, min.seg.width)
segments |
A |
min.seg.width |
The minimum segment width in base-pairs. |
The input model
with adjusted segments.
Aaron Taudt
## Load an HMM file <- list.files(system.file("extdata", "primary-lung", "hmms", package="AneuFinderData"), full.names=TRUE) hmm <- loadFromFiles(file)[[1]] ## Check number of segments before and after filtering length(hmm$segments) hmm$segments <- filterSegments(hmm$segments, min.seg.width=2*width(hmm$bins)[1]) length(hmm$segments)