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

## ----setup--------------------------------------------------------------------
library(delaydiscount)

## -----------------------------------------------------------------------------
remedi_rc_result <- jb_rule_check(remedi)
head(remedi_rc_result)

## -----------------------------------------------------------------------------
remedi_with_rc <- merge(remedi, remedi_rc_result)
remedi_rc_pass <- dplyr::filter(remedi_with_rc, 
                                C1 == TRUE, C2 == TRUE)

## -----------------------------------------------------------------------------
prep_remedi <- prepare_data_frame(remedi_rc_pass)

## -----------------------------------------------------------------------------
ln_k_ests <- get_subj_est_ln_k(prep_remedi)
head(ln_k_ests)

## -----------------------------------------------------------------------------
model_fit <- dd_hyperbolic_model(prep_remedi)

## -----------------------------------------------------------------------------
model_fit$ln_k_mean

## -----------------------------------------------------------------------------
model_fit$var

## -----------------------------------------------------------------------------
model_fit$model_test

## -----------------------------------------------------------------------------
model_fit$pairwise_f_tests