View source: R/lm_full_Bayes_SR.R
| lm_full_Bayes_SR | R Documentation |
Estimate regression coefficients and scale matrix for noise by using Gibbs MCMC algorithm. The function assumes 1) multivariate t-distribution for noise as a sampling distribution, and 2) noninformative priors for regression coefficients and scale matrix for noise.
lm_full_Bayes_SR(Y, X, dof = Inf, burnincycle = 1000, mcmccycle = 2000)
Y |
An N x K matrix of dependent variables. |
X |
An N x M matrix of regressors. |
dof |
Degree of freedom for multivariate t-distribution.
If |
burnincycle, mcmccycle |
Number of burnin cycles is the number of initially generated sample values to drop. Number of MCMC cycles is the number of generated sample values to compute estimates. |
Consider the multivariate regression:
Y = X \Psi + e, \quad e \sim MVT(0, \nu, \Sigma).
\Psi is a M-by-K matrix of regression coefficients and
\Sigma is a K-by-K scale matrix for multivariate t-distribution for
noise.
Sampling distribution for noise e is multivariate t-distribution with
degree of freedom dof and scale matrix \Sigma: e \sim MVT(0, \nu,
\Sigma).
The priors are noninformative priors: 1) the shrinkage prior for regression
coefficients \Psi, and 2) the reference prior for scale matrix
\Sigma.
The function implements Gibbs MCMC algorithm for estimating regression coefficients Psi and scale matrix Sigma.
A list object with estimated parameters: Psi, Sigma, dof, delta (delta is the reciprocal of lambda), and lambda. Additional components are se.param (standard error of the parameters) and LINEXVARmodel (estimates under LINEX loss).
S. Ni and D. Sun (2005). Bayesian estimates for vector autoregressive models. Journal of Business & Economic Statistics 23(1), 105-117.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.