Load a cross-sectional antibody survey data set

Description

Load a cross-sectional antibody survey data set

Usage

as_pop_data(
  data,
  antigen_isos = NULL,
  age = "Age",
  value = "result",
  id = "index_id",
  standardize = TRUE
)

Arguments

data a data.frame() or tibble::tbl_df
antigen_isos character() vector of antigen isotypes to be used in analyses
age a character() identifying the age column
value a character() identifying the value column
id a character() identifying the person ID column. The same ID must be used for every biomarker observation from the same person.
standardize a logical() to determine standardization of columns

Value

a pop_data object (a tibble::tbl_df with extra attribute antigen_isos)

Examples

Code
library("serocalculator")

library(magrittr)
xs_data <-
  serocalculator_example("example_pop_data.csv") |>
  read.csv() |>
  as_pop_data()

print(xs_data)
# A tibble: 200 × 8
   id    Country  cluster               catchment   age ageCat antigen_iso value
   <chr> <chr>    <chr>                 <chr>     <dbl> <chr>  <chr>       <dbl>
 1 P1    Pakistan Lyari Other           kgh        13.2 5-15   HlyE_IgA    0.568
 2 P2    Pakistan Gillani Railway Stat… aku        18   16+    HlyE_IgA    5.69 
 3 P3    Pakistan Gillani Railway Stat… aku         7.3 5-15   HlyE_IgA    1.23 
 4 P4    Pakistan Civic Centre          aku         2.6 <5     HlyE_IgA    1.08 
 5 P5    Pakistan Gillani Railway Stat… aku         3.9 <5     HlyE_IgA    1.43 
 6 P6    Pakistan Civic Centre          aku        13   5-15   HlyE_IgA    3.06 
 7 P7    Pakistan Machar Colony         kgh        11   5-15   HlyE_IgA    0.779
 8 P8    Pakistan Machar Colony         kgh        12   5-15   HlyE_IgA    1.90 
 9 P9    Pakistan Machar Colony         kgh        16   16+    HlyE_IgA    1.41 
10 P10   Pakistan Gillani Railway Stat… aku        14.4 5-15   HlyE_IgA    3.10 
# ℹ 190 more rows