canopy.BIC: To get BIC as a model selection criterion

In Canopy: Accessing Intra-Tumor Heterogeneity and Tracking Longitudinal and Spatial Clonal Evolutionary History by Next-Generation Sequencing

Description

To get BIC as a model selection criterion from MCMC sampling results.

Usage

canopy.BIC(sampchain,projectname,K,numchain,burnin,thin,pdf)

Arguments

sampchain	list of sampled trees returned by canopy.sample
projectname	name of project
K	number of subclones (vector)
numchain	number of MCMC chains with random initiations
burnin	burnin of MCMC chains
thin	MCMC chains thinning
pdf	whether a pdf plot of BIC should be generated, default to be TRUE

Value

BIC values (vector) for model selection with plot generated (pdf format).

Author(s)

Yuchao Jiang [e-mail]

Examples

data(MDA231_sampchain)
sampchain = MDA231_sampchain
projectname = 'MD231'
K = 3:6
numchain = 20
burnin = 150
thin = 5
bic = canopy.BIC(sampchain = sampchain, projectname = projectname, K = K,
                 numchain = numchain, burnin = burnin, thin = thin)

Canopy documentation built on May 1, 2019, 7:59 p.m.