VGAM (version 0.7-5)

ABO: The ABO Blood Group System

Description

Estimates the two independent parameters of the the ABO blood group system.

Usage

ABO(link = "logit", earg=list(), ir = NULL, ip = NULL)

Value

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

Details

The parameters p and q are probabilities, so that r=1-p-q is the third probability. The probabilities p and r correspond to A and O respectively, so that q is the probability for B. It is easier to make use of initial values for r than for q.

References

Lange, K. (2002) Mathematical and Statistical Methods for Genetic Analysis, 2nd ed. New York: Springer-Verlag.

See Also

AA.Aa.aa, AB.Ab.aB.ab, AB.Ab.aB.ab2, G1G2G3, MNSs.

Examples

y = cbind(A=725, B=258, AB=72, O=1073) # Order matters, not the name
fit = vglm(y ~ 1, ABO(link=logit), trace=TRUE, cri="coef")
fit = vglm(y ~ 1, ABO(link=identity), trace=TRUE, cri="coef")
coef(fit, matrix=TRUE)
Coef(fit) # Estimated p and q
rbind(y, sum(y)*fitted(fit))
sqrt(diag(vcov(fit)))