R/preprocess.hrd.R

In sztup/scarHRD: Estimating Genomic scars signatures (LST, TAI,LOH)

#' Preprocessing for further analysis
#'
#' @param seg segmentation data
#' @return preprocessed data
preprocess.hrd<-function(seg){
  seg <- seg[!seg[,2] %in% c(paste('chr',c('X','Y','x','y',23,24),sep=''),c('X','Y','x','y',23,24)),]
  seg[,1] <- as.character(seg[,1])

  if(! all(seg[,8] <= seg[,7]) ){
    tmp <- seg
    seg[tmp[,8] > tmp[,7],7]  <- tmp[tmp[,8] > tmp[,7],8]
    seg[tmp[,8] > tmp[,7],8]  <- tmp[tmp[,8] > tmp[,7],7]
  }
  seg <- shrink.seg.ai.wrapper(seg)
  return(seg)

}