summary.mpl {nlreg} | R Documentation |
Summary Method for ‘mpl’ Objects
Description
Returns a summary list for objects of class mpl
.
Usage
## S3 method for class 'mpl'
summary(object, correlation = FALSE, digits = NULL, ...)
Arguments
object |
a fitted |
correlation |
logical argument. If |
digits |
the number of significant digits to be printed. Defaults to
|
... |
absorbs any additional argument. |
Details
This function is a method for the generic function
summary
for
class mpl
. It can be invoked by calling summary
for an
object of the appropriate class, or directly by calling
summary.mpl
regardless of the class of the object.
Value
A list is returned with the following components:
varPar |
the maximum adjusted profile likelihood estimates of the variance parameters. |
coefficients |
the constrained MLEs of the regression coefficients given the maximum adjusted profile likelihood estimates of the variance parameters. |
offset |
the values passed through the |
varParMLE |
the MLEs of the variance parameters. |
coefMLE |
the MLEs of the regression coefficients. |
varParCov |
the (asymptotic) covariance matrix of the variance parameters, that is, the corresponding block in the inverse of the observed information matrix. |
coefCov |
the (asymptotic) covariance matrix of the regression coefficients, that is, the corresponding block in the inverse of the observed information matrix. |
lmp |
the adjusted profile log likelihood from the fit. |
lp |
the profile log likelihood from the fit. |
stats |
the indicator of which higher order solution was used. |
formula |
the model formula. |
meanFun |
the formula expression of the mean function. |
varFun |
the formula expression of the variance function. |
data |
a list representing a summary of the original data with the following components.
|
nobs |
the number of observations. |
iter |
the number of interations needed for convergence; only if
|
call |
an image of the call to |
ws |
the workspace component of the original
|
See Also
mpl.object
, nlreg.object
,
summary
Examples
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)) )
##
metsulfuron.mpl <- mpl( metsulfuron.nl, trace = TRUE )
summary( metsulfuron.mpl, corr = FALSE )