ScanBIC {seqCBS} | R Documentation |
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)
combX |
The number of reads at each unique read position |
combZ |
The number of case/tumor reads at each unique read position |
tauHat |
The change points called |
lik0 |
The null likelihood. Computed in the main routine. |
nTotal |
The total number of reads |
This is meanted to be called as a subrountine of ScanCBS
Returns a numerical value of mBIC for the current model
Jeremy J. Shen