Learn R Programming

Rfit (version 0.09)

rfit: Rank-based Estimates of Regression Coefficients

Description

Minimizes Jaeckel's dispersion function to obtain a rank-based solution for linear models.

Usage

rfit(formula, data = list(), ...)## S3 method for class 'default':
rfit(formula, data=list(), yhat0 = suppressWarnings(rq(y ~ x)$fitted.values), scores = wscores, symmetric = FALSE, ...)

Arguments

formula
an object of class formula
data
an optional data frame
yhat0
an n by vector of initial fitted values
scores
an object of class 'scores'
symmetric
logical. If 'FALSE' uses median of residuals as estimate of intercept
...
additional arguments to be passed to fitting routines

Value

  • coefficientsestimated regression coefficents with intercept
  • residualsthe residuals, i.e. y-yhat
  • fitted.valuesyhat = x betahat
  • xccentered design matrix
  • tauhatestimated value of the scale parameter tau
  • taushatestimated value of the scale parameter tau_s
  • betahatestimated regression coefficents
  • callCall to the function

Details

Rank-based estimation involves replacing the L2 norm of least squares estimation with a pseudo-norm which is a function of the ranks of the residuals. That is, in rank estimation, the usual notion of Euclidean distance is replaced with another measure of distance which is referred to as Jaeckel's (1972) dispersion function. Jaeckel's dispersion function depends on a score function and a library of commonly used score functions is included. e.g. Wilcoxon and sign score functions.

References

Hettmansperger, T.P. and McKean J.W. (1998), Robust nonparametric statistical methods, London: Arnold.

Jaeckel, L. A. (1972). Estimating regression coefficients by minimizing the dispersion of residuals. Annals of Mathematical Statistics, 43, 1449 - 1458.

Jureckova, J. (1971). Nonparametric estimate of regression coefficients. Annals of Mathematical Statistics, 42, 1328 - 1338.

See Also

summary.rfit

Examples

Run this code
data(baseball)
data(wscores)
fit<-rfit(weight~height,data=baseball)
summary(fit)

Run the code above in your browser using DataLab