## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(
  warning = FALSE,
  message = FALSE,
  fig.width = 7.252,
  fig.height = 4,
  comment = "#>",
  fig.retina = 3
)

## -----------------------------------------------------------------------------
library(logitr)
library(dplyr)
library(tidyr)
library(kableExtra) # For tables

numDraws <- unique(runtimes$numDraws)
logitr_time <- runtimes %>%
    filter(package == "logitr") %>%
    rename(time_logitr = time_sec)
time_compare <- runtimes %>%
    left_join(select(logitr_time, -package), by = "numDraws") %>%
    mutate(mult = round(time_sec/ time_logitr, 1)) %>%
    select(-time_logitr)
# Compare raw times
time_compare %>%
    select(-mult) %>%
    pivot_wider(names_from = numDraws, values_from = time_sec) %>% 
    kbl()
# Compare how many times slower compared to logitr
time_compare %>%
    select(-time_sec) %>%
    pivot_wider(names_from = numDraws, values_from = mult) %>% 
    kbl()

## ----eval=FALSE---------------------------------------------------------------
#  library(ggplot2)
#  library(ggrepel)
#  
#  plotColors <- c("black", RColorBrewer::brewer.pal(n = 5, name = "Set1"), "gold")
#  benchmark <- runtimes %>%
#      ggplot(aes(x = numDraws, y = time_sec, color = package)) +
#      geom_line() +
#      geom_point() +
#      geom_text_repel(
#          data = . %>% filter(numDraws == max(numDraws)),
#          aes(label = package),
#          hjust = 0, nudge_x = 40, direction = "y",
#          size = 4.5, segment.size = 0
#      ) +
#      scale_x_continuous(
#          limits = c(0, 1200),
#          breaks = numDraws,
#          labels = scales::comma) +
#      scale_y_continuous(limits = c(0, 300), breaks = seq(0, 300, 100)) +
#      scale_color_manual(values = plotColors) +
#      guides(
#          point = guide_legend(override.aes = list(label = "")),
#          color = guide_legend(override.aes = list(label = ""))) +
#      theme_bw(base_size = 18) +
#      theme(
#          panel.grid.minor = element_blank(),
#          panel.grid.major.x = element_blank(),
#          legend.position = "none",
#          axis.line.x = element_blank(),
#          axis.ticks.x = element_blank()
#      ) +
#      labs(
#          x = "Number of random draws",
#          y = "Computation time (seconds)"
#      )
#  
#  benchmark

## ----eval=FALSE, echo=FALSE---------------------------------------------------
#  ggsave('benchmark.png', benchmark, width = 8.5, height = 6)

## ----probabilities, echo=FALSE------------------------------------------------
knitr::include_graphics('benchmark.png')