sample_params_lm_gprior: Sample the parameters for a linear regression

In drkowal/rSTAR: MCMC and EM algorithms for Simultaneous Transformation and Rounding (STAR) Models

Sample the parameters for a linear regression

Description

Sample the parameters for a linear regression model assuming a g-prior for the coefficients.

Usage

sample_params_lm_gprior(y, X, params, psi = NULL, XtX = NULL)

Arguments

y	n x 1 vector of data
X	n x p matrix of predictors
params	the named list of parameters containing mu: vector of conditional means (fitted values) sigma: the conditional standard deviation coefficients: a named list of parameters that determine mu
psi	the prior variance for the g-prior
XtX	the p x p matrix of crossprod(X) (one-time cost); if NULL, compute within the function

Value

The updated named list params with draws from the full conditional distributions of sigma and coefficients (and updated mu).

Note

The parameters in coefficients are:

• beta: the p x 1 vector of regression coefficients components of beta

Examples

# Simulate data for illustration:
sim_dat = simulate_nb_lm(n = 100, p = 5)
y = sim_dat$y; X = sim_dat$X

# Initialize:
params = init_params_lm_gprior(y = y, X = X)

# Sample:
params = sample_params_lm_gprior(y = y, X = X, params = params)
names(params)
names(params$coefficients)

drkowal/rSTAR documentation built on July 5, 2023, 2:18 p.m.