hierarchical_fusion_BLR: Standard Hierarchical Monte Carlo Fusion for Bayesian...

In rchan26/BayesLogitFusion: Monte Carlo Fusion for Bayesian Logistic Regression

Description

Standard Hierarchical Monte Carlo Fusion for Bayesian Logistic Regression model

Usage

hierarchical_fusion_BLR(
  N_schedule,
  dim,
  y_split,
  X_split,
  prior_means,
  prior_variances,
  time_schedule,
  m_schedule,
  C,
  power,
  precondition = FALSE,
  L,
  base_samples,
  seed = NULL,
  n_cores = parallel::detectCores()
)

Arguments

N_schedule	vector of length (L-1), where N_schedule[l] is the number of samples per node at level l
dim	dimension of the predictors (= p+1)
y_split	list of length C, where y_split[[c]] is the y responses for sub-posterior c
X_split	list of length C, where X_split[[c]] is the design matrix for sub-posterior c
prior_means	prior for means of predictors
prior_variances	prior for variances of predictors
time_schedule	vector of legnth (L-1), where time_schedule[l] is the time chosen for Fusion at level l
m_schedule	vector of length (L-1), where m_schedule[l] is the number of samples to fuse for level l
C	number of sub-posteriors at the base level
power	exponent of target distribution
precondition	boolean value determining whether or not a preconditioning matrix is to be used
L	total number of levels in the hierarchy
base_samples	list of length C, where samples_to_fuse[c] containg the samples for the c-th node in the level
seed	seed number - default is NULL, meaning there is no seed
n_cores	number of cores to use

Value

A list with components:

samples

list of length (L-1), where samples[[l]][[i]] are the samples for level l, node i

time

list of length (L-1), where time[[l]] is the run time for level l, node i

rho_acc

list of length (L-1), where rho_acc[[l]][i] is the acceptance rate for first fusion step for level l, node i

Q_acc

list of length (L-1), where Q_acc[[l]][i] is the acceptance rate for second fusion step for level l, node i

rhoQ_acc

list of length (L-1), where rhoQ_acc[[l]][i] is the overall acceptance rate for fusion for level l, node i

y_inputs

input y data for each level and node

X_inputs

input X data for each level and node

C_inputs

vector of length (L-'), where C_inputs[l] is the number of sub-posteriors at level l+1 (the input for C to get to level l)

diffusion_times

vector of length (L-1), where diffusion_times[l] are the times for fusion in level l

power

exponent of target distributions in the hierarchy

precondition_matrices

pre-conditioning matricies that were used

rchan26/BayesLogitFusion documentation built on June 13, 2020, 5:03 a.m.