segment.baf.phased.multisample: Segment BAF, with the possible inclusion of structural...

View source: R/segmentation.R

segment.baf.phased.multisampleR Documentation

Segment BAF, with the possible inclusion of structural variant breakpoints

Description

This function breaks the genome up into chromosomes, possibly further when SV breakpoints are provided, and runs PCF on each to segment the chromosomes independently.

Usage

segment.baf.phased.multisample(
  samplename,
  inputfile,
  outputfile,
  prior_breakpoints_file = NULL,
  gamma = 10,
  calc_seg_baf_option = 3
)

Arguments

samplename

Name of the sample, which is used to name output figures

inputfile

String that points to the output from the combine.baf.files function. This contains the phased SNPs with their BAF values

outputfile

String where the segmentation output will be written

prior_breakpoints_file

String that points to a file with prior breakpoints (from SVs for example) with chromosome and position columns (Default: NULL)

gamma

The gamma parameter controls the size of the penalty of starting a new segment during segmentation. It is therefore the key parameter for controlling the number of segments (Default 10)

calc_seg_baf_option

Various options to recalculate the BAF of a segment. Options are: 1 - median, 2 - mean, 3 - ifelse median==0 or 1, median, mean. (Default: 3)

Author(s)

jdemeul, sd11


Wedge-Oxford/battenberg documentation built on Aug. 4, 2023, 6:27 p.m.