Calculate E_q(U_ijr) and E_q(log U_ijr) when U_ijr follows a Gamma(a_k,g,a_k,g)
qukg_calculation(s_ik = s_ik, t_jg = t_jg, x = x, mu_i = mu_i,
nu_j = nu_j, alpha_c = alpha_c, a = a)A list of 4 elements.
a_tildea_tilde.
b_tildeb_tilde.
exp_utildeexp_utilde.
exp_logutildeexp_logutilde.
s_ik.
t_jg.
a matrix of observations. Columns correspond to biological samples and rows to microorganisms.
mu_i.
nu_j.
alpha_c.
a0.