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

Description Usage Arguments Value Author(s) Examples

View source: R/infinite.R

Description

Fit an infinite Bayesian mixture model to IPD and sequence data

Usage

1
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 R@sherrillmix.com

Examples

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

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