inst/shiny/dropRates/rewards.R

library("shiny")
library("tidyverse")



rewards <- swgoh::rewards %>%
  
  dplyr::mutate(n_sims = ifelse(simulated, attempts, 0)) %>%
  dplyr::group_by(battle, reward, n_sims) %>%
  dplyr::summarize(
    n = sum(attempts) + sum(!simulated),
    y = sum(count),
    .groups = "keep") %>%
  dplyr::ungroup() %>%
  
  # left_join(swgoh::reward_details, by = "reward") %>%
  
  # Poisson model with Ga(0.5, 0.5) prior
  dplyr::mutate(
    # a = 0.5 + y,
    # b = 0.5 + n - y,
    # p  = y / n,
    # 
    # lb = qbeta(.025, a, b),
    # ub = qbeta(.975, a, b),
    # 
    # lb = ifelse(p==0, 0, lb),
    # ub = ifelse(p==0, qbeta(.95, a, b), ub),
    # 
    # lb = ifelse(p==1, qbeta(.05, a, b), lb),
    # ub = ifelse(p==1, 1, ub),
    
    a = 0.5 + y,
    b = 1 + n,
    rate  = y/n,
    
    lb = qgamma(.025, a, b),
    ub = qgamma(.975, a, b),
    
    lb = ifelse(y==0, 0, lb),
    ub = ifelse(y==0, qgamma(.95, a, b), ub),
    
    n_sims = as.integer(n_sims)) 
jarad/swgoh documentation built on March 10, 2021, 9:11 a.m.