VGAM (version 0.9-1)

G1G2G3: The G1G2G3 Blood Group System

Description

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

Usage

G1G2G3(link = "logit", ip1 = NULL, ip2 = NULL, iF = NULL)

Arguments

link

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

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.

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

Examples

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