Learn R Programming

seqCBS (version 1.2.1)

ScanStatNewComp: Main new window scan statistics computation

Description

This is a wrapper function to call the C routines for the scan statistic new candidate segmentation computing from the IGS

Usage

ScanStatNewComp(combZCumSum, combXCumSum, combZPoint, combXPoint, 
	p, nTotal, grid.cur, max.win, statistic)

Arguments

combZCumSum

A cumulative sum of the number of case reads

combXCumSum

A cumulative sum of the number of reads

combZPoint

The number of case reads at the grid points

combXPoint

The number of reads at the grid points

p

The proportion of case reads in the current region

nTotal

The total number of reads in the current region

grid.cur

The current grid to be computed on

max.win

The maximum inter-window to be considered for new scan

statistic

The type of statistic. Can be 'binomial','rabinowitz' or 'normal'

Value

Returns a matrix containing the candidate change points from the new scan

Details

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').

See Also

ScanCBS, ScanIterateGrid