A1A2A3(link = "logit", inbreeding = FALSE, ip1 = NULL, ip2 = NULL, iF = NULL)
p1
, p2
and f
.
See Links
for more choices.p1
, p2
and f
."vglmff"
(see vglmff-class
).
The object is used by modelling functions such as vglm
and vgam
.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$ and Hardy-Weinberg
Equilibrium (HWE) is assumed.AA.Aa.aa
,
AB.Ab.aB.ab
,
ABO
,
MNSs
.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")
Coef(fit) # Estimated p1 and p2
rbind(ymat, sum(ymat) * fitted(fit))
sqrt(diag(vcov(fit)))
Run the code above in your browser using DataLab