SSweibull                package:nls                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:

     data(ChickWeight)
     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)

