posterior_predictive_log_lambda: Posterior predictive of Poisson log-normal mixed model

Description Usage Arguments Value Examples

View source: R/posterior_predictive_log_lambda.R

Description

Posterior predictive of Poisson log-normal mixed model

Usage

1
posterior_predictive_log_lambda(obj, k = 1, show_donors = TRUE)

Arguments

obj

Object of class cytoeffect_poisson computed using poisson_lognormal

k

Draw from HMC chain

show_donors

Include donor random effect

Value

tibble object

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
set.seed(1)
df = simulate_data(n_cells = 10)
str(df)
fit = poisson_lognormal(df,
                        protein_names = names(df)[3:ncol(df)],
                        condition = "condition",
                        group = "donor",
                        r_donor = 2,
                        warmup = 200, iter = 325,
                        num_chains = 1)
posterior_predictive_log_lambda(fit, k = 1)

ChristofSeiler/cytoeffect documentation built on Nov. 24, 2020, 12:25 a.m.