mpl                  package:nlreg                  R Documentation

_M_a_x_i_m_u_m _A_d_j_u_s_t_e_d _P_r_o_f_i_l_e _L_i_k_e_l_i_h_o_o_d _E_s_t_i_m_a_t_i_o_n - _G_e_n_e_r_i_c _F_u_n_c_t_i_o_n

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

     Calculates the maximum adjusted profile likelihood estimates.

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

     mpl(fitted, ...)

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

  fitted: any fitted model object for which the maximum adjusted
          profile likelihood estimates can be calculated. 

     ...: absorbs any additional argument. 

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

     This function is generic (see 'methods'); method  functions can be
     written to handle specific classes of data.  Classes which already
     have methods for this function include: 'nlreg'.

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

     the maximum adjusted profile likelihood estimates for all
     parameters of a regression model or for a subset of them.

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

     'mpl.nlreg', 'nlreg.object', 'methods'

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

     data(metsulfuron)
     metsulfuron.nl <- 
         nlreg( formula = log(area) ~ log( b1+(b2-b1) / (1+(dose/b4)^b3) ), 
                weights = ~ ( 1+dose^exp(g) )^2, data = metsulfuron, hoa = TRUE, 
                start = c(b1 = 138, b2 = 2470, b3 = 2, b4 = 0.07, g = log(0.3)) )
     mpl( metsulfuron.nl, trace = TRUE )
     ##
     options( object.size = 10000000 )
     data(chlorsulfuron)
     chlorsulfuron.nl <- 
         nlreg( log(area) ~ log( b1+(b2-b1) / (1+(dose/b4)^b3) ), 
             weights = ~ ( 1+k*dose^g*(b2-b1)^2/(1+(dose/b4)^b3)^4*b3^2*dose^(2*b3-2)/
                         b4^(2*b3)/(b1+(b2-b1)/(1+(dose/b4)^b3))^2 ),
             start = c(b1 = 2.2, b2 = 1700, b3 = 2.8, b4 = 0.28, g = 2.7, k = 1), 
             data = chlorsulfuron, hoa = TRUE, trace = TRUE,  
             control = list(x.tol = 10^-12, rel.tol = 10^-12, step.min = 10^-12) )
     mpl( chlorsulfuron.nl, trace = TRUE )

