invparalogistic             package:VGAM             R Documentation

_I_n_v_e_r_s_e _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  inverse
     paralogistic distribution.

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

     invparalogistic(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 inverse paralogistic distribution is the
     4-parameter generalized beta II distribution with shape parameter
     q=1 and a=p. It is the 3-parameter Dagum distribution with a=p.
     More details can be found in Kleiber and Kotz (2003).

     The inverse paralogistic distribution has density

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

     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(a + 1/a)  gamma(1 - 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:

     'Invparalogistic', 'genbetaII', 'betaII', 'dagum', 'sinmad',
     'fisk', 'invlomax', 'lomax', 'paralogistic'.

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

     y = rinvparalogistic(n=3000, 4, 6)
     fit = vglm(y ~ 1, invparalogistic, trace=TRUE)
     fit = vglm(y ~ 1, invparalogistic(init.a=2.7, init.sc=3.3),
                trace=TRUE, crit="c")
     coef(fit, mat=TRUE)
     Coef(fit)
     summary(fit)

