binary_preds.R

library(tidyverse)
expand.grid(day = 1:7, week = 1:13) %>%
    mutate(mu = 3*sin(day*2*pi/7)) %>%
    mutate(p = boot::inv.logit(mu)) %>%
    mutate(obs  = rbinom(n(), size = 1, p = p),
           pred = rbinom(n(), size = 1, p = p)) %>%
    gather(key, value, obs, pred) %>%
    mutate(value = factor(value)) %>%
    ggplot(data =., aes(x = day, y = key)) +
    geom_tile(aes(fill = value)) +
    coord_equal() +
    facet_grid(week ~ .)