SinglefPrior: Samples from the single f prior.

In HWEBayes: Bayesian investigation of Hardy-Weinberg Equilibrium via estimation and testing.

Description

Function to sample from the single f prior, that is the Dirichlet and normal on λ, where λ=\log((f-f_{\min})/(1-f)).

Usage

SinglefPrior(nsim, alpha, lambdamu, lambdasd)

Arguments

nsim	number of simulations from prior
alpha	vector of k parameters for the Dirichlet prior on the k allele frequencies.
lambdamu	mean of the normal prior on λ.
lambdasd	standard deviation of the normal prior on λ.

Value

p	sample for vector of k allele frequencies
f	sample of f parameters
lgts	samples for logits of baseline logits
lambda	samples for λ

Author(s)

Jon Wakefield ([e-mail])

References

Wakefield, J. (2010). Bayesian methods for examining Hardy-Weinberg equilibrium. Biometrics; Vol 66:257-65

See Also

SinglefReject, HWEsimdat

Examples

SinglefSamp <- SinglefPrior(nsim=1000,alpha=c(1,1,1,1),
  lambdamu=-2.95,lambdasd=1.07)

HWEBayes documentation built on May 2, 2019, 11:07 a.m.