paralogistic              package:VGAM              R Documentation

_P_a_r_a_l_o_g_i_s_t_i_c _D_i_s_t_r_i_b_u_t_i_o_n _F_a_m_i_l_y _F_u_n_c_t_i_o_n

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

     Maximum likelihood estimation of the 2-parameter  paralogistic
     distribution.

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

     paralogistic(link.a = "loge", link.scale = "loge",
                  earg.a=list(), earg.scale=list(),
                  init.a = 1, init.scale = NULL, zero = NULL)

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

link.a, link.scale: Parameter link functions applied to the (positive)
          shape parameter 'a' and (positive) scale parameter 'scale'.
          See 'Links' for more choices.

earg.a, earg.scale: List. Extra argument for each of the links. See
          'earg' in 'Links' for general information.

init.a, init.scale: Optional initial values for 'a' and 'scale'.

    zero: An integer-valued vector specifying which linear/additive
          predictors are modelled as intercepts only. Here, the values
          must be from the set {1,2} which correspond to 'a', 'scale',
          respectively.

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

     The 2-parameter paralogistic distribution is the 4-parameter
     generalized beta II distribution with shape parameter p=1 and a=q.
     It is the 3-parameter Singh-Maddala distribution with a=q. More
     details can be found in Kleiber and Kotz (2003).

     The 2-parameter paralogistic has density

            f(y) = a^2 y^(a-1) / [b^a (1 + (y/b)^a)^(1+a)]

     for a > 0, b > 0, y > 0. Here, b is the scale parameter 'scale',
     and a is the shape parameter. The mean is

         E(Y) = b  gamma(1 + 1/a)  gamma(a - 1/a) /  gamma(a)

     provided a > 1.

_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:

     If the self-starting initial values fail, try experimenting with
     the initial value arguments, especially those whose default value
     is not 'NULL'.

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

     T. W. Yee

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

     Kleiber, C. and Kotz, S. (2003) _Statistical Size Distributions in
     Economics and Actuarial Sciences_, Hoboken, NJ:
     Wiley-Interscience.

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

     'Paralogistic', 'genbetaII', 'betaII', 'dagum', 'fisk',
     'invlomax', 'lomax', 'invparalogistic'.

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

     y = rparalogistic(n=3000, 4, 6)
     fit = vglm(y ~ 1, paralogistic, trace=TRUE)
     fit = vglm(y ~ 1, paralogistic(init.a=2.3, init.sc=5),
                trace=TRUE, crit="c")
     coef(fit, mat=TRUE)
     Coef(fit)
     summary(fit)

