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

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="logit"), trace=TRUE)
fit = vglm(y ~ 1, AA.Aa.aa(link="probit"), trace=TRUE)
fit = vglm(y ~ 1, AA.Aa.aa(link="cloglog", init.p=0.9), trace=TRUE)
fit = vglm(y ~ 1, AA.Aa.aa(link="identity"), trace=TRUE)
rbind(y, sum(y)*fitted(fit))
Coef(fit) # Estimated pA
summary(fit)