G1G2G3                 package:VGAM                 R Documentation

_T_h_e _G_1_G_2_G_3 _B_l_o_o_d _G_r_o_u_p _S_y_s_t_e_m

_D_e_s_c_r_i_p_t_i_o_n:

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

_U_s_a_g_e:

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

_A_r_g_u_m_e_n_t_s:

    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'.

_D_e_t_a_i_l_s:

     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.

_V_a_l_u_e:

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

_N_o_t_e:

     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.

_A_u_t_h_o_r(_s):

     T. W. Yee

_R_e_f_e_r_e_n_c_e_s:

     Lange, K. (2002) _Mathematical and Statistical Methods for Genetic
     Analysis_, 2nd ed. New York: Springer-Verlag.

_S_e_e _A_l_s_o:

     'AA.Aa.aa', 'AB.Ab.aB.ab', 'AB.Ab.aB.ab2', 'ABO', 'MNSs'.

_E_x_a_m_p_l_e_s:

     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, .3, .3, .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)))

