split_clusters_3 {rhierbaps}

split_clusters_3

Description

Peform an iteration of the third move in the algorithm. That is split cluster in two and re-allocate one sub-cluster.

Usage

split_clusters_3(
  snp.object,
  partition,
  threshold = 1e-05,
  min.clust.size = 20,
  n.cores = 1
)

Arguments

snp.object

A snp.object containing the processed SNP data.

partition

An integer vector indicating an initial partition of the isolates.

threshold

The increase in marginal log likelihood required to accept a move.

min.clust.size

Clusters smaller than min.clust.size will not be split.

n.cores

The number of cores to use.

Value

The best partition after splitting a cluster and re-allocating as well as a boolean value indicating whether a move increased the marginal likelihood.