VGAM (version 0.9-3)

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", init.pA = NULL)

Arguments

link

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

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)