Description Usage Arguments Value Author(s) See Also Examples
View source: R/merging.segments.R
It is an option to merge adjacent segments, for which the median values in either or both log2ratio and log2mBAF dimensions are not substantially different. For WGS and SNP array, it is recommended to do so.
1 2 3 4 5 6 7 |
data |
a data frame containing log2ratio and log2mBAF data generated
by |
segs.stat |
a data frame containing segment locations and summary statistics
resulting from |
use.null.data |
logical. If only data for probes located in normal copy
segments to be used for bootstrapping. Default is |
N |
the number of replicates drawn by bootstrap. |
maxL |
integer. The maximum length in terms of number of probes a bootstrapped segment
may span. Default is |
merge.pvalue.cutoff |
a p-value cut-off for merging. If the empirical p-value is greater than the cut-off value, the two adjacent segments under consideration will be merged. |
do.manual.baseline |
logical. If baseline adjustment to be done manually. Default is |
log2mBAF.left, log2mBAF.right, log2ratio.bottom, log2ratio.up |
left, right, bottom and up boundaries to be specified manually by a visual inspectio of
2-D diagnosis plot generated by |
seed |
integer. Random seed can be set for reproducibility of results. |
verbose |
logical. If more details to be output. Default is |
A data frame with the same columns as the one generated by
joint.segmentation
.
Zhongyang Zhang <zhongyang.zhang@mssm.edu>
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | data(seq.data)
data(seq.segs)
## Not run:
seq.segs.merge <- merging.segments(data=seq.data, segs.stat=seq.segs,
use.null.data=TRUE,
N=1000, maxL=2000,
merge.pvalue.cutoff=0.05, verbose=TRUE)
## End(Not run)
## how the results look like
data(seq.segs.merge)
head(seq.segs.merge)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.