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

In mskcc/facets: Cellular Fraction and Copy Numbers from Tumor Sequencing

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=1e-3)

Arguments

x	the output from procSample. This function uses the elements jointseg, out and dipLogR from the output.
trace	flag to print the EM criteria at every step
unif	random EM start values of cellular fractions instead of clusteredcncf values
min.nhet	minimum number of heterozygote snps in a segment used to call minor cn
maxiter	maximum number of EM iterations
eps	the 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)

mskcc/facets documentation built on Oct. 15, 2021, 3:12 p.m.