cf_reh_exact: Exact confidence bounds to the relative excess heterozygosity...

In EQUIVNONINF: Testing for Equivalence and Noninferiority

Description

Implementation of the interval estimation procedure described on pp. 305-6 of Wellek S (2010) Testing statistical hypotheses of equivalence and noninferiority. Second edition.

Usage

cf_reh_exact(X1,X2,X3,alpha,SW,TOL,ITMAX)

Arguments

X1	count of homozygotes of the first kind [<-> genotype AA]
X2	count of heterozygotes [<-> genotype AB]
X3	count of homozygotes of the second kind [<-> genotype BB]
alpha	1 - confidence level
SW	width of the search grid for determining an interval covering the parameter point at which the conditional distribution function takes value α and 1-α, respectively
TOL	numerical tolerance to the deviation between the computed confidence limits and their exact values
ITMAX	maximum number of interval-halving steps

Details

The program exploits the structure of the family of all genotype distributions, which is 2-parameter exponential with \log(REH) as one of these parameters.

Value

X1	count of homozygotes of the first kind [<-> genotype AA]
X2	count of heterozygotes [<-> genotype AB]
X3	count of homozygotes of the second kind [<-> genotype BB]
alpha	1 - confidence level
SW	width of the search grid for determining an interval covering the parameter point at which the conditional distribution function takes value α and 1-α, respectively
TOL	numerical tolerance to the deviation between the computed confidence limits and their exact values
ITMAX	maximum number of interval-halving steps
C_l_exact	exact conditional lower (1-α)-confidence bound to REH
C_r_exact	exact conditional upper (1-α)-confidence bound to REH

Author(s)

Stefan Wellek <[e-mail]>
Peter Ziegler <[e-mail]>

References

Wellek S, Goddard KAB, Ziegler A: A confidence-limit-based approach to the assessment of Hardy-Weinberg equilibrium. Biometrical Journal 52 (2010), 253-270.

Wellek S: Testing statistical hypotheses of equivalence and noninferiority. Second edition. Boca Raton: Chapman & Hall/CRC Press, 2010, Par. 9.4.3.

Examples

               
cf_reh_exact(34,118,96,.05,.1,1E-4,25)

EQUIVNONINF documentation built on July 12, 2021, 5:08 p.m.