HWLindley: Calculate a posteriori density for Lindley's alpha

In HardyWeinberg: Statistical Tests and Graphics for Hardy-Weinberg Equilibrium

Calculate a posteriori density for Lindley's alpha

Description

Function HWLindley calculates the posterior density for disequilibrium measure alpha, as defined by Lindley (1988).

Usage

HWLindley(alphaseq = seq(-3, 3, by = 0.01), x)

Arguments

alphaseq	a single value or a sequence of values for alpha
x	the genotype count vector in format (AA,AB,BB)

Details

Numerical integration is used to compute the density.

Value

a vector with values of the density for each value in alphaseq

Author(s)

Jan Graffelman [e-mail]

References

Lindley, D.V. (1988) Statistical Inference Concerning Hardy-Weinberg Equilibrium. In: Bernardo, J.M., DeGroot, M.H., Lindley, D.V. and Smith, A.F.M. Bayesian Statistics, 3, pp. 307-326. Oxford University Press.

See Also

HWPosterior

Examples

x <- c(MM=298,MN=489,NN=213)
post.dens <- HWLindley(seq(-2,2,by=0.01),x)
## Not run: 
  plot(seq(-2,2,by=0.01),post.dens,type="l")

## End(Not run)

HardyWeinberg documentation built on May 7, 2022, 5:05 p.m.