infiniteGaussian: Fit an infinite Bayesian mixture model to IPD and sequence...

In sherrillmix/infiniter: Bayesian infinite mixture models

Description

Fit an infinite Bayesian mixture model to IPD and sequence data

Usage

infiniteGaussian(ipds, seqs, nIter = 100, maxK = 100, mc.cores = 4)

Arguments

ipds	a matrix of IPD ratios with a row for each position of interest and a column for each surrounding base
seqs	a matrix of single characters or vector of strings giving sequences of interest
nIter	number of JAGS iterations
maxK	the maximum number of clusters expected (smaller is less computationally expensive)
mc.cores	number of cores (and chains) to use

Value

list with matrix giving simulated data and stored JAGS model (can be used for further iterations)

Author(s)

Scott Sherrill-Mix [e-mail]

Examples

fakeSeqs<-generateFakeClusters()
infiniteGaussian(fakeSeqs$ipds,fakeSeqs$seqs,maxK=10)

sherrillmix/infiniter documentation built on May 29, 2019, 9:24 p.m.