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

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

ThetaMater.M1: (Bayesian Model 2) MCMC function that uses the MCMCmetrop1R function from MCMCpack

Description

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

Usage

ThetaMater.M1(k.vec, l.vec, n.vec, c.vec, 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
ngens	Number of generations to run the MCMC
burnin	Number of generations to discard as burnin
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.