Main function

Main function

zinb_w_cov(Y_mat, X_mat, z_vec, s_vec, tau_mukj = 1, S = 20000L,
  burn_rate = 0.5, tau_mu0 = 1, tau_phi = 1, tau_beta = 1,
  a_omega = 1, b_omega = 1, a_pi = 1, b_pi = 1, a_p = 1,
  b_p = 1, a_mu = 2, b_mu = 10, a_phi = 1, b_phi = 0.01,
  a_beta = 2, b_beta = 10, mu0_mean = 5, beta_lim = 5)

Arguments

Y_mat,	abundance matrix
X_mat,	covariate matrix
z_vec,	group vector
s_vec,	scaling factor
tau_mukj	default = 1,
S	default = 20000,
burn_rate	default = 0.5,
tau_mu0	= 1,
tau_phi	= 1,
tau_beta	= 1,
a_omega	default = 0.2, uniform prior
b_omega	defeault = 1.8, uniform prior
a_pi	default = 1, control the sparsity of R matrix
b_pi	default = 0.1, control the sparsity of R matrix
a_p	default = 0.4, control the sparsity of Delta matrix
b_p	default = 1.6, control the sparsity of Delta matrix
a_mu	default = 2, spike and slab for mu_kj
b_mu	default = 15, spike and slab for mu_kj
a_phi	default = 10, dispersion parameter
b_phi	default = 1, dispersion parameter
a_beta	default = 2, spike and slab for beta_rj
b_beta	default = 15, spike and slab for beta_rj
mu0_mean,	mean
beta_lim	= 5, limit of beta