# get_est_cpp: Summarize the posterior mean, sd and confidence interval... In zhenkewu/lotR: Latent class analysis with Observed Trees in R (lotR)

 get_est_cpp R Documentation

## 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.