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

## ----setup--------------------------------------------------------------------
library(tempodisco)

## -----------------------------------------------------------------------------
data("td_bc_single_ptpt")
mod <- td_bcnm(td_bc_single_ptpt, discount_function = 'hyperbolic')
plot(mod, type = 'summary')

## -----------------------------------------------------------------------------
plot(mod, type = 'summary', verbose = F, log = 'x')

## -----------------------------------------------------------------------------
plot(mod, type = 'summary', verbose = F, log = 'x', p_lines = c(0.1, 0.9))

## -----------------------------------------------------------------------------
data("td_ip_simulated_ptpt")
mod_ip <- td_ipm(td_ip_simulated_ptpt, discount_function = 'hyperbolic')
plot(mod_ip, type = 'summary', log = 'x', verbose = F)

## -----------------------------------------------------------------------------
mod_ip <- kirby_score(td_bc_single_ptpt)
plot(mod_ip, type = 'summary', log = 'x', verbose = F)

## -----------------------------------------------------------------------------
plot(mod, type = 'endpoints')

## -----------------------------------------------------------------------------
del <- sort(unique(mod$data$del))[2]
plot(mod, type = 'endpoints', del = del, val_del = 1000)

## -----------------------------------------------------------------------------
plot(mod, type = 'link')

## -----------------------------------------------------------------------------
ddm <- td_ddm(td_bc_single_ptpt, discount_function = 'exponential',
              gamma_par_starts = 0.01,
              beta_par_starts = 0.5,
              alpha_par_starts = 3.5,
              tau_par_starts = 0.9)
plot(ddm, type = 'rt', q_lines = c(0.05, 0.95), legend = T)