This computes mBIC for the current change point model. We then use this to determine the appropriate model complexity.
ScanBIC(combX, combZ, tauHat, lik0, nTotal)
The number of reads at each unique read position
The number of case/tumor reads at each unique read position
The change points called
The null likelihood. Computed in the main routine.
The total number of reads
Returns a numerical value of mBIC for the current model
This is meanted to be called as a subrountine of ScanCBS