[BioC] Tuning segment() in DNAcopy package to get more segments on each chromosome
littleduck24 at gmail.com
Mon Nov 14 21:13:11 CET 2011
I am trying to tune segment() in DNAcopy package to generate more
segments. Currently I have few segments in each chromosome (min 3
I want to have a least 20 segments in each chromosome. min 5 markers
in each segment. also the number of markers in moving window =1000.
My 1st question is how the parameter kmax and nmin work in the cbs
algorithm? I read the literature, but I couldn't figure it out.
My 2nd question is in the segment function min.width=5 for min 5
markers in each segment. So how to set the window size = 1000?
kmax is the maximum width of smaller segment for permutation in the
nmin is the minimum length of data for which the approximation of
maximum statistic is used under the hybrid method. should be larger
segment(x, alpha = 0.01, nperm = 10000, p.method =
"hybrid", min.width=5, kmax=25, nmin=200,
undo.splits = "prune",undo.prune=0.05,
More information about the Bioconductor