https://github.com/cran/quantreg
Tip revision: ff94aa39c643db8e130ab6d9f2f3bff5d7f97cc8 authored by Roger Koenker on 13 April 2012, 00:00:00 UTC
version 4.79
version 4.79
Tip revision: ff94aa3
nlrq.Rd
\name{nlrq}
\alias{nlrq}
\alias{nlrqModel}
\alias{print.nlrq}
\alias{summary.nlrq}
\alias{deviance.nlrq}
\alias{formula.nlrq}
\alias{coef.nlrq}
\alias{fitted.nlrq}
\alias{logLik.nlrq}
\alias{AIC.nlrq}
\alias{extractAIC.nlrq}
\alias{predict.nlrq}
\alias{print.summary.nlrq}
\alias{tau.nlrq}
\title{ Function to compute nonlinear quantile regression estimates}
\description{
This function implements an R version of an interior point method
for computing the solution to quantile regression problems which
are nonlinear in the parameters. The algorithm is based on interior
point ideas described in Koenker and Park (1994).
}
\usage{
nlrq(formula, data=parent.frame(), start, tau=0.5,
control, trace=FALSE,method="L-BFGS-B")
\method{summary}{nlrq}(object, ...)
\method{print}{summary.nlrq}(x, digits = max(5, .Options$digits - 2), ...)
}
\arguments{
\item{formula}{ formula for model in nls format; accept self-starting models }
\item{data}{ an optional data frame in which to evaluate the variables in
`formula' }
\item{start}{a named list or named numeric vector of starting estimates }
\item{tau}{ a vector of quantiles to be estimated}
\item{control}{ an optional list of control settings. See `nlrq.control' for
the names of the settable control values and their effect.}
\item{trace}{ logical value indicating if a trace of the iteration progress
should be printed. Default is `FALSE'. If `TRUE' intermediary results
are printed at the end of each iteration. }
\item{method}{ method passed to optim for line search, default is "L-BFGS-B"
but for some problems "BFGS" may be preferable. See \code{\link{optim}} for
further details. Note that the algorithm wants to pass
upper and lower bounds for the line search to optim, which is fine for
the L-BFGS-B method. Use of other methods will produce warnings about
these arguments -- so users should proceed at their own risk.}
\item{object}{an object of class nlrq needing summary.}
\item{x}{an object of class summary.nlrq needing printing.}
\item{digits}{Significant digits reported in the printed table.}
\item{...}{Optional arguments passed to printing function.}
}
\details{An `nlrq' object is a type of fitted model object. It has methods
for the generic functions `coef' (parameters estimation at best solution),
`formula' (model used), `deviance' (value of the objective function at best
solution), `print', `summary', `fitted' (vector of fitted variable according
to the model), `predict' (vector of data points predicted by the model, using
a different matrix for the independent variables) and also for the function
`tau' (quantile used for fitting the model, as the tau argument of the
function). Further help is also available for the method `residuals'.
The summary method for nlrq uses a bootstrap approach based on the final
linearization of the model evaluated at the estimated parameters.
}
\value{
A list consisting of:
\item{m}{an `nlrqModel' object similar to an `nlsModel' in package nls}
\item{data }{the expression that was passed to `nlrq' as the data argument.
The actual data values are present in the environment of the
`m' component. }
}
\author{Based on S code by Roger Koenker modified for R and to accept models
as specified by nls by Philippe Grosjean.}
\references{ Koenker, R. and Park, B.J. (1994). An Interior Point Algorithm for
Nonlinear Quantile Regression, Journal of Econometrics, 71(1-2): 265-283.
}
\seealso{ \code{\link{nlrq.control}} , \code{\link{residuals.nlrq}} }
\examples{
# build artificial data with multiplicative error
Dat <- NULL; Dat$x <- rep(1:25, 20)
set.seed(1)
Dat$y <- SSlogis(Dat$x, 10, 12, 2)*rnorm(500, 1, 0.1)
plot(Dat)
# fit first a nonlinear least-square regression
Dat.nls <- nls(y ~ SSlogis(x, Asym, mid, scal), data=Dat); Dat.nls
lines(1:25, predict(Dat.nls, newdata=list(x=1:25)), col=1)
# then fit the median using nlrq
Dat.nlrq <- nlrq(y ~ SSlogis(x, Asym, mid, scal), data=Dat, tau=0.5, trace=TRUE)
lines(1:25, predict(Dat.nlrq, newdata=list(x=1:25)), col=2)
# the 1st and 3rd quartiles regressions
Dat.nlrq <- nlrq(y ~ SSlogis(x, Asym, mid, scal), data=Dat, tau=0.25, trace=TRUE)
lines(1:25, predict(Dat.nlrq, newdata=list(x=1:25)), col=3)
Dat.nlrq <- nlrq(y ~ SSlogis(x, Asym, mid, scal), data=Dat, tau=0.75, trace=TRUE)
lines(1:25, predict(Dat.nlrq, newdata=list(x=1:25)), col=3)
# and finally "external envelopes" holding 95 percent of the data
Dat.nlrq <- nlrq(y ~ SSlogis(x, Asym, mid, scal), data=Dat, tau=0.025, trace=TRUE)
lines(1:25, predict(Dat.nlrq, newdata=list(x=1:25)), col=4)
Dat.nlrq <- nlrq(y ~ SSlogis(x, Asym, mid, scal), data=Dat, tau=0.975, trace=TRUE)
lines(1:25, predict(Dat.nlrq, newdata=list(x=1:25)), col=4)
leg <- c("least squares","median (0.5)","quartiles (0.25/0.75)",".95 band (0.025/0.975)")
legend(1, 12.5, legend=leg, lty=1, col=1:4)
}
\keyword{models}
\keyword{regression}
\keyword{nonlinear}
\keyword{robust}