autoplot.curve_params: graph antibody decay curves by antigen isotype

In serocalculator: Estimating Infection Rates from Serological Data

graph antibody decay curves by antigen isotype

Description

graph antibody decay curves by antigen isotype

Usage

## S3 method for class 'curve_params'
autoplot(
  object,
  antigen_isos = unique(object$antigen_iso),
  ncol = min(3, length(antigen_isos)),
  ...
)

Arguments

object	a data.frame() of curve parameters (one or more MCMC samples)
antigen_isos	antigen isotypes to analyze (can subset curve_params)
ncol	how many columns of subfigures to use in panel plot
...	Arguments passed on to plot_curve_params_one_ab verbose verbose output xlim range of x values to graph n_curves how many curves to plot (see details). n_points Number of points to interpolate along the x axis (passed to ggplot2::geom_function()) rows_to_graph which rows of curve_params to plot (overrides n_curves). alpha (passed to ggplot2::geom_function()) how transparent the curves should be: 0 = fully transparent (invisible) 1 = fully opaque log_x should the x-axis be on a logarithmic scale (TRUE) or linear scale (FALSE, default)? log_y should the Y-axis be on a logarithmic scale (default, TRUE) or linear scale (FALSE)?

Details

rows_to_graph

Note that if you directly specify rows_to_graph when calling this function, the row numbers are enumerated separately for each antigen isotype; in other words, for the purposes of this argument, row numbers start over at 1 for each antigen isotype. There is currently no way to specify different row numbers for different antigen isotypes; if you want to do that, you could call plot_curve_params_one_ab() directly for each antigen isotype and combine the resulting panels yourself. Or you could subset curve_params manually, before passing it to this function, and set the n_curves argument to Inf.

Value

a ggplot2::ggplot() object

Examples

library(dplyr)
library(ggplot2)
library(magrittr)

curve <-
  serocalculator_example("example_curve_params.csv") %>%
  read.csv() %>%
  as_curve_params() %>%
  filter(antigen_iso %in% c("HlyE_IgA", "HlyE_IgG")) %>%
  autoplot()

curve

serocalculator documentation built on April 3, 2025, 7:35 p.m.