This is a wrapper function to call the C routines for the scan statistic new candidate segmentation computing from the IGS
ScanStatNewComp(combZCumSum, combXCumSum, combZPoint, combXPoint,
p, nTotal, grid.cur, max.win, statistic)
A cumulative sum of the number of case reads
A cumulative sum of the number of reads
The number of case reads at the grid points
The number of reads at the grid points
The proportion of case reads in the current region
The total number of reads in the current region
The current grid to be computed on
The maximum inter-window to be considered for new scan
The type of statistic. Can be 'binomial','rabinowitz' or 'normal'
Returns a matrix containing the candidate change points from the new scan
The computations are done in C for speed. The three types of statistics are by the use of exact binomial likelihood ('binomial'), score statistic ('rabinowitz') or using normal approximation to the binomial ('normal').