SSweibull               package:stats               R Documentation

_W_e_i_b_u_l_l _g_r_o_w_t_h _c_u_r_v_e _m_o_d_e_l

_D_e_s_c_r_i_p_t_i_o_n:

     This 'selfStart' model evaluates the Weibull model for growth
     curve data and its gradient.  It has an 'initial' attribute that
     will evaluate initial estimates of the parameters 'Asym', 'Drop',
     'lrc', and 'pwr' for a given set of data.

_U_s_a_g_e:

     SSweibull(x, Asym, Drop, lrc, pwr)

_A_r_g_u_m_e_n_t_s:

       x: a numeric vector of values at which to evaluate the model.

    Asym: a numeric parameter representing the horizontal asymptote on
          the right side (very small values of 'x').

    Drop: a numeric parameter representing the change from 'Asym' to
          the 'y' intercept.

     lrc: a numeric parameter representing the natural logarithm of the
          rate constant.

     pwr: a numeric parameter representing the power to which 'x' is
          raised.

_D_e_t_a_i_l_s:

     This model is a generalization of the 'SSasymp' model in that it
     reduces to 'SSasymp' when 'pwr' is unity.

_V_a_l_u_e:

     a numeric vector of the same length as 'x'.  It is the value of
     the expression 'Asym-Drop*exp(-exp(lrc)*x^pwr)'.  If all of the
     arguments 'Asym', 'Drop', 'lrc', and 'pwr' are names of objects,
     the gradient matrix with respect to these names is attached as an
     attribute named 'gradient'.

_A_u_t_h_o_r(_s):

     Douglas Bates

_R_e_f_e_r_e_n_c_e_s:

     Ratkowsky, David A. (1983), _Nonlinear Regression Modeling_,
     Dekker. (section 4.4.5)

_S_e_e _A_l_s_o:

     'nls', 'selfStart', 'SSasymp'

_E_x_a_m_p_l_e_s:

     Chick.6 <- subset(ChickWeight, (Chick == 6) & (Time > 0))
     SSweibull(Chick.6$Time, 160, 115, -5.5, 2.5 )  # response only
     Asym <- 160; Drop <- 115; lrc <- -5.5; pwr <- 2.5
     SSweibull(Chick.6$Time, Asym, Drop, lrc, pwr)  # response and gradient
     getInitial(weight ~ SSweibull(Time, Asym, Drop, lrc, pwr), data = Chick.6)
     ## Initial values are in fact the converged values
     fm1 <- nls(weight ~ SSweibull(Time, Asym, Drop, lrc, pwr), data = Chick.6)
     summary(fm1)

