get_est_cpp: Summarize the posterior mean, sd and confidence interval...

In zhenkewu/lotR: Latent class analysis with Observed Trees in R (lotR)

Summarize the posterior mean, sd and confidence interval (grouped or individual leaf nodes)

Description

Summarize the posterior mean, sd and confidence interval (grouped or individual leaf nodes)

Usage

get_est_cpp(
  prob,
  mu_gamma,
  sigma_gamma,
  mu_alpha,
  Sigma_alpha,
  anc,
  cardanc,
  z
)

Arguments

prob	a vector of variational probability. Length = p. At the extremes, it can also be a vector of zeros and ones, indicating which nodes are selected based on variational probability (prob > 0.5).
mu_gamma	variational Gaussian means (for s_u=1 component) for J*K logit(class-specific response probabilities); (J,K,p) array; In R, we used a list of p (J,K) matrices
sigma_gamma	variational Gaussian variances (for s_u=1 component) for J*K logit(class-specific response probabilities); (J,K,p) array
mu_alpha	variational Gaussian mean vectors (for s_u=1 component) - this is a p by K-1 matrix; in R, we used a list of p vectors (each of length K-1)
Sigma_alpha	variational Gaussian variances (for s_u=1 component) this is an array of dimension (K-1, K-1, p); in R, we used a list of p matrices, each of dimension K-1 by K-1.
anc	a list of pL vectors, each vector has the node ids of the ancestors; lengths may differ. The ancestors include the node concerned.
cardanc	a numeric vector of length pL; integers. The number of ancestors for each leaf node
z	= ci_level+(1-ci_level)/2

Value

a list

Named("beta_est")=beta_est, Named("beta_sd_est")=beta_sd_est, Named("beta_cil")=beta_cil, Named("beta_ciu")=beta_ciu,

Named("eta_est")=eta_est, Named("eta_var_est")=eta_var_est, Named("eta_cil")=eta_cil, Named("eta_ciu")=eta_ciu

zhenkewu/lotR documentation built on April 24, 2022, 2:36 a.m.