Learn R Programming

Canopy (version 1.0.0)

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

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).

Examples

Run this code
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)

Run the code above in your browser using DataLab