The QUANTLIFE Procedure

Confidence Interval

Direct computation of the covariance of the parameter estimators involves a complicated density estimation. Instead, the QUANTLIFE procedure computes confidence intervals for the quantile regression parameters beta left-parenthesis tau right-parenthesis by using resampling methods. The QUANTLIFE procedure implements two different methods, the exponentially weighted method and the pairwise resampling method.

Exponentially Weighted Method

This method samples weights w Subscript i Baseline comma i equals 1 comma ellipsis comma n comma from a standard exponential distribution that has mean 1 and variance 1. Then it computes the censored quantile regression estimators ModifyingAbove beta With caret left-parenthesis tau right-parenthesis based on the observed data left-parenthesis x Subscript i Baseline comma upper Y Subscript i Baseline comma normal upper Delta Subscript i Baseline right-parenthesis with the weights w Subscript i. These steps are repeated B times (where B is the value of the NREP= option in the PROC QUANTLIFE statement). The confidence intervals can be obtained from these B estimates. You can specify this method by using the CI=EW option in the PROC QUANTLIFE statement.

Pairwise Method

This method samples left-parenthesis x Subscript i Baseline comma upper Y Subscript i Baseline comma normal upper Delta Subscript i Baseline right-parenthesis with replacement and computes the quantile regression estimators ModifyingAbove beta With caret left-parenthesis tau right-parenthesis based on the resampled data. These steps are repeated B times (where B is the value of the NREP= option in the PROC QUANTLIFE statement). The confidence intervals can be obtained from these B estimates. You can specify this method by using the CI=PW option in the PROC QUANTLIFE statement.

Testing Effects of Covariates

Consider the linear model

y Subscript i Baseline equals bold x prime Subscript 1 i Baseline bold-italic beta 1 plus bold x prime Subscript 2 i Baseline bold-italic beta 2 plus epsilon Subscript i

where bold-italic beta 1 and bold-italic beta 2 are p-dimensional and q-dimensional parameters, respectively, and epsilon Subscript i, i equals 1 comma ellipsis comma n, are errors. Denote bold x prime Subscript i Baseline equals left-parenthesis bold x prime Subscript 1 i Baseline comma bold x prime Subscript 2 i right-parenthesis, and let ModifyingAbove bold-italic beta With caret Subscript 1 Baseline left-parenthesis tau right-parenthesis and ModifyingAbove bold-italic beta With caret Subscript 2 Baseline left-parenthesis tau right-parenthesis be the parameter estimates for bold-italic beta 1 and bold-italic beta 2, respectively, at the tauth quantile.

The QUANTLIFE procedure implements the Wald test for the null hypothesis:

upper H 0 colon beta 2 left-parenthesis tau right-parenthesis equals 0

The Wald test statistic, which is based on the estimated coefficients ModifyingAbove beta With caret Subscript 2 from the unrestricted fitted model, is given by

upper T Subscript upper W Baseline left-parenthesis tau right-parenthesis equals ModifyingAbove beta With caret prime Subscript 2 Baseline left-parenthesis tau right-parenthesis ModifyingAbove normal upper Sigma With caret left-parenthesis tau right-parenthesis Superscript negative 1 Baseline ModifyingAbove beta With caret Subscript 2 Baseline left-parenthesis tau right-parenthesis

where ModifyingAbove normal upper Sigma With caret left-parenthesis tau right-parenthesis is an estimator of the covariance of ModifyingAbove beta With caret Subscript 2 Baseline left-parenthesis tau right-parenthesis, which is obtained by using resampling methods.

Last updated: December 09, 2022