VGAM (version 0.9-6)

A1A2A3: The A1A2A3 Blood Group System

Description

Estimates the three independent parameters of the the A1A2A3 blood group system.

Usage

A1A2A3(link = "logit", inbreeding = TRUE, ip1 = NULL, ip2 = NULL, iF = NULL)

Arguments

link

Link function applied to p1, p2 and f. See Links for more choices.

inbreeding

Logical. Is the inbreeding coefficient $f$ equal to 0?

ip1, ip2, iF

Optional initial value for p1, p2 and f.

Value

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

Details

The parameters p1 and p2 are probabilities, so that p3=1-p1-p2 is the third probability. The parameter f is the third independent parameter if inbreeding = TRUE. If inbreeding = FALSE then $f = 0$.

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, ABO, MNSs.

Examples

ymat <- cbind(108, 196, 429, 143, 513, 559)
fit <- vglm(ymat ~ 1, A1A2A3(link = probit), trace = TRUE, crit = "coef")
fit <- vglm(ymat ~ 1, A1A2A3(link = logit, ip1 = 0.3, ip2 = 0.3, iF = 0.02),
            trace = TRUE, crit = "coef")
fit <- vglm(ymat ~ 1, A1A2A3(link = "identitylink"), trace = TRUE)
Coef(fit)  # Estimated p1, p2 and f
rbind(ymat, sum(ymat) * fitted(fit))
sqrt(diag(vcov(fit)))