R/findBreakPoints.R

In snapCGH: Segmentation, normalisation and processing of aCGH data

findBreakPoints <- function (seg, array) 
{
  states <- seg$state[,array]
  chr <- seg$genes$Chr
    bpoints <- c(1)
    for (i in 2:length(states)) {
        if (states[i] != states[i - 1] || chr[i] != chr[i -1])
            bpoints <- c(bpoints, i - 1, i)
    }
    bpoints <- c(bpoints, length(states))
    bpoints
}