## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----setup--------------------------------------------------------------------
library(serosv)
library(magrittr)

## -----------------------------------------------------------------------------
linelisting <- parvob19_fi_1997_1998[order(parvob19_fi_1997_1998$age), ]
aggregated <- hav_bg_1964

# View the 2 different data format
head(linelisting)
head(aggregated)

# fit with aggregated data
model1 <- polynomial_model(aggregated, k=1)
plot(model1)
# fit with linelisting data
model2 <- polynomial_model(linelisting, k=1, status_col = "seropositive")
plot(model2)

## -----------------------------------------------------------------------------
transform_data(
  linelisting,
  stratum_col="age", 
  status_col="seropositive") %>% 
  polynomial_model(k=1) %>% 
  plot()