gibbs_2sampleDAG0: Implementation of Gibbs step for two-sample DAG0

In junsoukchoi/BayesDAG0: Bayesian Differential Zero-Inflated Negative Binomial Directed Acyclic Graphs

Implementation of Gibbs step for two-sample DAG0

Description

Implementation of Gibbs step for two-sample DAG0

Usage

gibbs_2sampleDAG0(param, tuning, priors, x0, x1, temp = 1)

Arguments

param	a list with each tag corresponding to a parameter name. Valid tags are 'E0', 'D', 'E1', 'U', 'zeta', 'eta' 'alpha0', 'alpha1', 'beta0', 'beta1', 'delta0', 'delta1', 'gamma0', 'gamma1', 'psi0', 'psi1', 'tau', and 'rho'. The value portion of each tag is the parameter values to be updated.
tuning	a list with each tag corresponding to a parameter name. Valid tags are 'E', 'E_rev', 'U', 'zeta', 'eta', 'alpha0', 'alpha1', 'beta0', 'beta1', 'delta0', 'delta1', 'gamma0', 'gamma1', 'psi0', and 'psi1'. The value portion of each tag defines the standard deviations of Normal proposal distributions for Metropolis sampler for each parameter. The tag 'E' is for the birth or death update schemes for E = {E_0, E_1}, while 'E_rev' is for the reversal schemes. The value portion of both tags should consist of the standard deviations of Gaussian proposals for alpha, beta, delta, gamma, zeta, and eta for an update of E. Also, the value portion of the tag 'U' should include the standard deviations of Gaussian proposals for zeta and eta for an update of U.
priors	a list with each tag corresponding to a parameter name. Valid tags are 'psi', 'tau', and 'rho'. The value portion of each tag defines the hyperparameters for two-sample DAG0.
x0	a matrix containing data for the control group.
x1	a matrix containing data for the case/treatment group.
temp	a temperature.

Value

a list with the tags 'param' and 'acceptance'. The value portion of the tag 'param' gives parameter values after performing Gibbs step for two-sample DAG0. The value portion of the tag 'acceptance' shows the acceptance indicators for each parameter.

junsoukchoi/BayesDAG0 documentation built on Feb. 26, 2025, 2:37 p.m.