estimate_m_rd: Estimate heterozygote miscall rate for different read depth...
In whoa: Evaluation of Genotyping Error in Genotype-by-Sequencing Data

To see how this Rcpp function is applied, see the code in infer_m.

1	estimate_m_rd(Y, R, init_m, num_cats, p_prior, m_prior, num_reps)

`Y`	the 012,-1 matrix that is N x L giving the observed genotypes of the N individuals at L SNPs.
`R`	integer matrix that is N x L giving the read depth categories. These must be indexed from 1 up to num_cats. Missing data should be -1.
`init_m`	starting value for m. Typically you might want to use the m estimated from init_m
`num_cats`	the number of read depth categories.
`p_prior`	two-vector that holds the beta parameters for a prior on the allele frequency for each locus. Typical value is c(0.5, 0.5).
`m_prior`	two-vector that holds the beta parameters for a prior on the heterozygote miscall rate for each locus. Typical value is c(0.5, 0.5).
`num_reps`	the number of MCMC sweeps to do.