Qtools (version 1.0)

sparsity.rqt: Sparsity Estimation

Description

This function estimates the density and sparsity functions of the residuals from a rq or a rqt object.

Usage

sparsity(object, se = "nid", hs = TRUE)
## S3 method for class 'rq':
sparsity(object, se = "nid", hs = TRUE)
## S3 method for class 'rqs':
sparsity(object, se = "nid", hs = TRUE)
## S3 method for class 'rqt':
sparsity(object, se = "nid", hs = TRUE)

Value

sparsity returns an object of class list that contains three elements:densityestimate of the density of the residuals.sparsityestimate of the sparsity of the residuals.bandwidthbandwidth used for estimation.

Details

This function is based on the code from quantreg::summary.rq and quantreg::bandwidth.rq to estimate the sparsity function for linear quantile regression models (Koenker and Bassett, 1978) and transformation models of Geraci and Jones (2014).

References

Geraci M and Jones MC. Improved transformation-based quantile regression. Canadian Journal of Statistics 2015;43(1):118-132. Koenker R (2013). quantreg: Quantile Regression. R package version 5.05. URL http://CRAN.R-project.org/package=quantreg. Koenker R, Bassett G. Regression quantiles. Econometrica. 1978;46(1):33-50. Powell JL. Estimation of monotonic regression models under quantile restrictions. In: Barnett W, Powell J, Tauchen G, editors. Nonparametric and Semiparametric Methods in Econometrics and Statistics: Proceedings of the Fifth International Symposium on Economic Theory and Econometrics. New York, NY: Cambridge University Press 1991. p. 357-84.

See Also

rq

Examples

data(trees)

# 'rqt' object

fit.rqt <- tsrq(Volume ~ Height, tsf = "bc", symm = FALSE, data = trees,
	lambda = seq(-10, 10, by = 0.01), tau = 0.5)
sparsity(fit.rqt)

# 'rq' object

fit.rq <- rq(Volume ~ Height, data = trees)
sparsity(fit.rq, se = "iid")
sparsity(fit.rq, se = "nid")
sparsity(fit.rq, se = "ker")