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

## -----------------------------------------------------------------------------
# Load required packages
library(dplyr)
library(tidyr)
library(ggplot2)
library(joinpointR)

# Load example data
data(hiv_data)

## -----------------------------------------------------------------------------
mod1 <- model_jp(
  data = hiv_data,
  value = hiv_rate,
  time = year,
  group = "sex",
  k = 3,
  step = TRUE,
  test = TRUE
)

# Show the model output
mod1

## -----------------------------------------------------------------------------
mod2 <- model_jp(
  data = hiv_data,
  value = hiv_rate,
  time = year,
  group = "sex",
  k = 1,
  step = FALSE,
  test = FALSE
)

# Show the model output
mod2

## -----------------------------------------------------------------------------
mod3 <- model_jp(
  data = hiv_data,
  value = hiv_rate,
  time = year,
  group = c("region", "sex"),
  step = TRUE
)

## -----------------------------------------------------------------------------
get_apc(mod2)

## -----------------------------------------------------------------------------
get_apc(mod1, digits = 1)

## -----------------------------------------------------------------------------
get_aapc(mod1, digits = 1, show_ci = TRUE)

## -----------------------------------------------------------------------------
get_aapc(mod1, digits = 1, show_ci = FALSE)

## -----------------------------------------------------------------------------
summary_jp(mod1)

## -----------------------------------------------------------------------------
# id: tab-2
mod2 |>
  summary_jp(digits = 1) |>
  jp_to_ft()

## -----------------------------------------------------------------------------
mod3 |>
  summary_jp(digits = 1) |>
  jp_to_ft()

## -----------------------------------------------------------------------------
mod3 |>
  summary_jp(digits = 1) |>
  jp_to_ft(lan = "es")

## -----------------------------------------------------------------------------
gg_jpoint(mod1)

## -----------------------------------------------------------------------------
gg_jpoint(mod1, facets = "none")

## -----------------------------------------------------------------------------
gg_jpoint(mod3)

## -----------------------------------------------------------------------------
gg_jpoint(mod3, facets = "grid")

## -----------------------------------------------------------------------------
gg_jpoint(mod3, facets = "grid", psize = 5)

## -----------------------------------------------------------------------------
gg_jpoint(mod3, facets = "grid", ptr = 1)

## -----------------------------------------------------------------------------
gg_jpoint(mod3, facets = "grid", obs = FALSE)

## -----------------------------------------------------------------------------
gg_jpoint(mod3, facets = "grid", jp = FALSE)

## -----------------------------------------------------------------------------
gg_jpoint(mod3, facets = "grid", jp = FALSE, cbpal = "blue_fluoride")

## -----------------------------------------------------------------------------
gg_jpoint(mod3, facets = "grid", cb = FALSE)