emcncf: EM estimate of copy number and cellular fraction of segment...

In rptashkin/facets2n: Cellular Faction and Copy Numbers from Tumor Sequencing

EM estimate of copy number and cellular fraction of segment clusters

Description

Uses genotype mixture model to estimate the cluster specific copy number and cellular fraction. Uses estimates based on the cnlr.median and mafR as initial values for the EM iteration.

Usage

emcncf(
  x,
  trace = FALSE,
  unif = FALSE,
  min.nhet = 15,
  maxiter = 10,
  eps = 0.001
)

Arguments

x	(list) the output from procSample. This function uses the elements jointseg, out and dipLogR from the output.
trace	(logical) flag to print the EM criteria at every step
unif	(logical) random EM start values of cellular fractions instead of clusteredcncf values
min.nhet	(numeric) minimum number of heterozygote snps in a segment used to call minor cn
maxiter	(numeric) maximum number of EM iterations
eps	the (numeric) convergence threshold

Value

A list containing:

loglik	loglikelihood value of the fitted model
purity	fraction tumor cells in the tumor sample
ploidy	average total copy number of the tumor cells
dipLogR	estimated logR value of diploid segments
cncf	dataframe consisting of the columns of segmentation output as well as cellular fraction (cf), total (tcn) and lesser (lcn) copy number of each segment and their em counterpart (with .em suffix)

rptashkin/facets2n documentation built on May 11, 2022, 1:34 p.m.