ThetaMater.M2: ThetaMater.M2: (Bayesian Model 2) MCMC function that uses the...

In radamsRHA/ThetaMater: ThetaMater: Bayesian estimation of the population size parameter theta from genomic data

ThetaMater.M2: (Bayesian Model 2) MCMC function that uses the MCMCmetrop1R function from MCMCpack with a discrete gamma rate variation model

Description

This function returns a list of each step in the MCMC sampling chain

Usage

ThetaMater.M2(k.vec, l.vec, n.vec, c.vec, K, alpha.param, ngens, burnin, thin,
  theta.shape, theta.scale)

Arguments

k.vec	Vector of mutation counts
l.vec	Vector of locus lengths
n.vec	Vector of sample numbers
c.vec	Vector of data pattern counts
K	Number of classese to approximate the gamma distribution
alpha.param	Shape of the gamma distribution for describing the amount of among-locus rate variation
ngens	Number of generations to run the MCMC
burnin	Number of generations to discard from MCMC chain
thin	Number of generations to thin in the MCMC chain
theta.shape	Shape parameter of the gamma distribution for setting the prior on theta
theta.scale	Scale parameter of the gamma distribution for setting the prior on theta

Examples

library(Rcpp)
library(ThetaMater)
library(MCMCpack)

radamsRHA/ThetaMater documentation built on April 25, 2024, 10:11 a.m.