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

## ----setup, output=FALSE, message=FALSE, warning=FALSE------------------------
library(serosv)
library(dplyr)

## -----------------------------------------------------------------------------
pl <- parvob19_be_2001_2003 %>% 
  penalized_spline_model(status_col = "seropositive", s = "tp", framework = "pl") 
pl

## ----fig.width=7, fig.height=4------------------------------------------------
plot(pl)

## -----------------------------------------------------------------------------
glmm <- parvob19_be_2001_2003 %>% 
  penalized_spline_model(status_col = "seropositive", s = "tp", framework = "glmm") 
glmm

## ----fig.width=7, fig.height=4------------------------------------------------
plot(glmm)