spike_and_slab_logistic_moretrees: Group spike and slab variable selection with Gaussian outcome

Description Usage Arguments Details Value Model Description See Also

View source: R/spike_and_slab_logistic_moretrees.R

Description

Here's a brief description. spike_and_slab_logistic_moretrees performs group variable selection via a spike and slab prior for binary data. The posterior is approximated via variational inference. This function returns the parameters of the variational approximation.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
spike_and_slab_logistic_moretrees(
  dsgn,
  initial_values,
  random_init,
  tol,
  max_iter,
  update_hyper,
  update_hyper_freq,
  hyper_fixed,
  print_freq,
  hyper_random_init,
  vi_random_init
)

Arguments

tol

Convergence tolerance for ELBO.

update_hyper

Update hyperparameters? Default = TRUE.

update_hyper_freq

How frequently to update hyperparameters. Default = every 10 iterations.

y

Integer vector of length n containing outcomes; 1 = success, 0 = failure.

X

Matrix of dimension n x sum(K), where n is the number of units, and K[g] is the number of variables in group g.

W

Matrix of non-sparse regression coefficients of dimension n x m

model

"ss" or "moretrees". If "ss", regular spike and slab variable selection is implemented. If "moretrees", the multi-outcome moretrees model is fitted.

groups

If model = "ss", groups is a list of length G (number of groups), where groups[[g]] is an integer vector specifying the columns of X that belong to group g. NULL if model = "moretrees".

outcomes_units:

If model = "moretrees", outcomes_units is a list of length equal to the number of unique outcomes. Each element of the list is an integer vector indicating which units (entries of y_reord, rows of X_reord) correspond to each outcomes. NULL if model = "ss".

outcomes_nodes:

If model = "moretrees", outcomes_nodes is a list of length equal to the number of unique nodes. Each element of the list is an integer vector indicating which outcomes/leaves are descendants of each node. NULL if model = "ss".

ancestors

If model = "moretrees", ancestors ia list of length equal to the number of unique outcomes. Each element of the list is an integer vector indicating which nodes on the tree (including leaves) are ancestors of the corresponding outcome. NULL if model = "ss".

maxiter

Maximum number of iterations of the VI algorithm.

Details

All the details go here!

Value

A list of variational parameters.

Model Description

Describe group spike and slab prior and all parameters here.

See Also

Other spike and slab functions: compute_betas(), compute_thetas(), ml_by_group(), spike_and_slab_logistic(), spike_and_slab_normal_moretrees(), spike_and_slab_normal(), spike_and_slab()


IQSS/moretrees documentation built on March 20, 2020, 8:44 p.m.