DiSTATIS plot of posterior or bootstrap samples

plot_distatis(
  obj,
  ndraws = 1000,
  ncores = 1,
  show_donors = TRUE,
  show_markers = TRUE,
  repel = TRUE
)

Arguments

obj: Object of class cytoeffect_poisson or cytoeffect_poisson_mcle computed by poisson_lognormal or poisson_lognormal_mcle
ndraws: Number of posterior or bootstrap samples
ncores: Number of cores
show_donors: Include donor random effect
show_markers: Include markers
repel: Repel marker names

Value

ggplot2 object

Examples

if (FALSE) {
set.seed(1)
df = simulate_data()
str(df)
fit = poisson_lognormal(df,
                        protein_names = names(df)[3:ncol(df)],
                        condition = "condition",
                        group = "donor",
                        r_donor = 2,
                        warmup = 200, iter = 325, adapt_delta = 0.95,
                        num_chains = 1)
plot_distatis(fit, ndraws = 125)
fit = poisson_lognormal_mcle(df,
                             protein_names = names(df)[3:ncol(df)],
                             condition = "condition",
                             group = "donor",
                             ncores = 1)
plot_distatis(fit, ndraws = 125)
}