G1G2G3 {VGAM} | R Documentation |
Estimates the three independent parameters of the the G1G2G3 blood group system.
G1G2G3(link = "logit", earg=list(), ip1 = NULL, ip2 = NULL, iF = NULL)
link |
Link function applied to p1 , p2 and f .
See Links for more choices.
|
earg |
List. Extra argument for the link.
See earg in Links for general information.
|
ip1, ip2, iF |
Optional initial value for p1 , p2 and f .
|
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.
An object of class "vglmff"
(see vglmff-class
).
The object is used by modelling functions such as vglm
and vgam
.
The input can be a 6-column matrix of counts,
with columns corresponding to
G_1/G_1
,
G_1/G_2
,
G_1/G_3
,
G_2/G_2
,
G_2/G_3
,
G_3/G_3
(in order).
Alternatively, the input can be a 6-column matrix of
proportions (so each row adds to 1) and the weights
argument is used to specify the total number of counts for each row.
T. W. Yee
Lange, K. (2002) Mathematical and Statistical Methods for Genetic Analysis, 2nd ed. New York: Springer-Verlag.
AA.Aa.aa
,
AB.Ab.aB.ab
,
AB.Ab.aB.ab2
,
ABO
,
MNSs
.
y = cbind(108, 196, 429, 143, 513, 559) fit = vglm(y ~ 1, G1G2G3(link=probit), trace=TRUE, crit="coef") fit = vglm(y ~ 1, G1G2G3(link=logit, ip1=.3, ip2=.3, iF=.02), trace=TRUE, crit="coef") fit = vglm(y ~ 1, G1G2G3(link="identity"), trace=TRUE) Coef(fit) # Estimated p1, p2 and f rbind(y, sum(y)*fitted(fit)) sqrt(diag(vcov(fit)))