VGAM (version 0.8-1)

AA.Aa.aa: The AA-Aa-aa Blood Group System

Description

Estimates the parameter of the AA-Aa-aa blood group system.

Usage

AA.Aa.aa(link = "logit", earg=list(), init.pA = NULL)

Arguments

link

Link function applied to pA. See Links for more choices.

earg

List. Extra argument for the link. See earg in Links for general information.

init.pA

Optional initial value for pA.

Value

• An object of class "vglmff" (see vglmff-class). The object is used by modelling functions such as vglm and vgam.

Details

This one parameter model involves a probability called pA. The probability of getting a count in the first column of the input (an AA) is pA*pA.

References

Weir, B. S. (1996) Genetic Data Analysis II: Methods for Discrete Population Genetic Data, Sunderland, MA: Sinauer Associates, Inc.

See Also

AB.Ab.aB.ab, AB.Ab.aB.ab2, ABO, G1G2G3, MNSs.

Examples

y = cbind(53, 95, 38)
fit = vglm(y ~ 1, AA.Aa.aa(link="probit"), trace=TRUE)
rbind(y, sum(y)*fitted(fit))
Coef(fit) # Estimated pA
summary(fit)