Description Usage Arguments Details Value References Examples
Bayesian quantile regression based on asymmetric Laplace likelihood with posterior variance adjustment
1 2 |
y |
the response vector |
x |
the design matrix. If the first column of x is not all ones, a column of ones will be added. |
tau |
the quantile level of interest |
niter |
integer: number of iterations to run the chain for. Default 20000. |
burn_in |
integer: discard the first burn_in values. Default 100. |
prop_cov |
covariance matrix giving the covariance of the proposal distribution. This matrix need not be positive definite. If the covariance structure of the target distribution is known (approximately), it can be given here. If not given, the diagonal will be estimated via the Fisher information matrix. |
level |
nominal confidence level for the credible interval |
The function returns the unadjusted and adjusted posterior standard deviation, and unadjusted and adjusted credible intervals for Bayesian quantile regression based on asymmetric Laplace working likelihood.
A list of the following commponents is returned
estpar: posterior mean of the regression coefficient vector
PSD: posterior standard deviation without adjustment
PSD.adj: posterior standard deviation with adjustment
CI.BAL: credible interval without adjustment
CI.BAL.adj: credible interval with adjustment
sig: estimated scale parameter
MCMCsize: effective size of the chain
Yang, Y., Wang, H. and He, X. (2015). Posterior inference in Bayesian quantile regression with asymmetric Laplace likelihood. International Statistical Review, 2015. doi: 10.1111/insr.12114.
1 2 3 4 5 6 7 8 9 10 |
Loading required package: quantreg
Loading required package: SparseM
Attaching package: 'SparseM'
The following object is masked from 'package:base':
backsolve
Loading required package: MHadaptive
Loading required package: MASS
Loading required package: coda
Loading required package: survival
Attaching package: 'survival'
The following object is masked from 'package:quantreg':
untangle.specials
$estpar
[1] -0.06855049 1.99545736 1.97272196
$PSD
[1] 0.08616764 0.07904206 0.09566932
$PSD.adj
[1] 0.09914312 0.07985734 0.11146693
$CI.BAL
LB UB
[1,] -0.2102836 0.07318266
[2,] 1.8654447 2.12546999
[3,] 1.8153599 2.13008398
$CI.BAL.adj
LB UB
[1,] -0.2316264 0.09452543
[2,] 1.8641037 2.12681100
[3,] 1.7893752 2.15606874
$sig
[1] 0.5300029
$MCMCsize
var1 var2 var3
1308.693 1293.016 1290.231
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.