sampleFastGaussian

Sample from N(mu, Sigma) where Sigma = solve(crossprod(Phi) + solve(D))
and mu = Sigma*crossprod(Phi, alpha):

Monte Carlo sampling algorithms for semiparametric
Bayesian regression analysis. These models feature a nonparametric
(unknown) transformation of the data paired with widely-used
regression models including linear regression, spline regression,
quantile regression, and Gaussian processes. The transformation
enables broader applicability of these key models, including for
real-valued, positive, and compactly-supported data with challenging
distributional features. The samplers prioritize computational
scalability and, for most cases, Monte Carlo (not MCMC) sampling for
greater efficiency. Details of the methods and algorithms are provided
in Kowal and Wu (2024) <doi:10.1080/01621459.2024.2395586>.

Dan Kowal

SeBR

Semiparametric Bayesian Regression Analysis

sampleFastGaussian function

<dl><dt>Phi</dt>
<dd><code>n x p</code> matrix (of predictors)</dd>
<dt>Ddiag</dt>
<dd><code>p x 1</code> vector of diagonal components (of prior variance)</dd>
<dt>alpha</dt>
<dd><code>n x 1</code> vector (of data, scaled by variance)</dd></dl>

Arguments

Sample a Gaussian vector using Bhattacharya et al. (2016) — sampleFastGaussian

<dl>

<dt>Phi</dt>
<dd><code>n x p</code> matrix (of predictors)</dd>


<dt>Ddiag</dt>
<dd><code>p x 1</code> vector of diagonal components (of prior variance)</dd>


<dt>alpha</dt>
<dd><code>n x 1</code> vector (of data, scaled by variance)</dd>

</dl>

sampleFastGaussian: Sample a Gaussian vector using Bhattacharya et al. (2016)

Description

Usage

Value

Arguments

References