estimate_m: heterozygotes miscall rate

In thierrygosselin/radiator: RADseq Data Exploration, Manipulation and Visualization using R

heterozygotes miscall rate

Description

estimate the heterozygotes miscall rate from a simple model

Usage

estimate_m(
  data,
  nreps = 200,
  m_init = stats::runif(1),
  a0 = 0.5,
  a1 = 0.5,
  sm = 0.005
)

Arguments

data	A tidy data frame (biallelic only).
nreps	number of MCMC sweeps to do. Default: nreps = 200.
m_init	initial starting value for m must be between 0 and 1
a0	beta parameter for reference alleles
a1	beta parameter for alternate alleles
sm	standard devation of proposal distribution for m

Author(s)

Eric Anderson [e-mail]

thierrygosselin/radiator documentation built on Nov. 7, 2024, 1:30 p.m.