VGAM (version 0.9-4)

AB.Ab.aB.ab2: The AB-Ab-aB-ab2 Blood Group System

Description

Estimates the parameter of the the AB-Ab-aB-ab2 blood group system.

Usage

AB.Ab.aB.ab2(link = "logit", init.p = NULL)

Arguments

link

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

init.p

Optional initial value for p.

Value

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

Warning

There may be a bug in the deriv and weight slot of the family function.

Details

This one parameter model involves a probability called p.

References

Elandt-Johnson, R. C. (1971) Probability Models and Statistical Methods in Genetics, New York: Wiley.

See Also

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

Examples

ymat <- cbind(68, 11, 13, 21)  # See Elandt-Johnson, pp.430,427
fit <- vglm(ymat ~ 1, AB.Ab.aB.ab2(link = cloglog), trace = TRUE, crit = "coef")
Coef(fit)  # Estimated p
rbind(ymat, sum(ymat) * fitted(fit))
sqrt(diag(vcov(fit)))  # Estimated variance is approx 0.0021