malacoda_intervals: Plot a variant's posterior intervals
In andrewGhazi/malacoda: Bayesian Analysis of High-Throughput Genomic Assays

Description Usage Arguments Details See Also Examples

Given a result object from the malacoda model, this function plots the posterior intervals for all the modeled parameters. The outer black bars represent the limits of the 95 bar represents the 80 posterior means.

malacoda_intervals(
  sampler_res,
  inner_probs = c(0.1, 0.9),
  outer_probs = c(0.025, 0.975)
)

`sampler_res`	A stanfit object for one variant
`inner_probs`	quantiles of the posterior for the inner bar
`outer_probs`	quantiles of the posterior for the outer bar

The sampler_res object for a given variant is saved into the out_dir argument of fit_mpra_model. The barcode concentrations (i.e. DNA means) are numbered in the same order in which they are given to fit_mpra_model, with poorly represented barcodes removed entirely.

fit_mpra_model get_well_represented