Getting Started: Compute Doodle match probability • doodleproblem

library(doodleproblem)
library(tidyverse)
library(tictoc)

Parameters

n_colleagues <- 3
p <- 3
o <- 10
dep <- 0

Doodle match probability

prob_doodle_match(n_colleagues = 5,
                  o = 10,
                  p = 3,
                  dep = 0)
#> [1] 0.023

We can also a somewhat more detailled output by setting return_prob to FALSE:

prob_doodle_match(n_colleagues = 5,
                  o = 10,
                  p = 3,
                  dep = 0,
                  return_prob = FALSE)
#> # A tibble: 2 × 3
#>    Freq n_matches prop_matches
#>   <int>     <int>        <dbl>
#> 1     5        26        0.026
#> 2    NA       974        0.974

Doodle samples 1

smples1 <- doodle_samples(n_colleagues, o, p)
glimpse(smples1)
#> Rows: 1,018
#> Columns: 4
#> $ slots_chosen <fct> 6, 7, NA, NA, NA, 2, NA, NA, NA, NA, 4, 10, 8, NA, NA, NA…
#> $ Freq         <int> 3, 3, NA, NA, NA, 3, NA, NA, NA, NA, 3, 3, 3, NA, NA, NA,…
#> $ is_match     <lgl> TRUE, TRUE, FALSE, FALSE, FALSE, TRUE, FALSE, FALSE, FALS…
#> $ id           <dbl> 1, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 10, 11, 12, 13, 14, 15,…

The parameters given to the function call can be found using attr afterwards:

attr(smples1, "n_colleagues")
#> [1] 3
attr(smples1, "o")
#> [1] 10

A similar dataset is available as data:

data(smples1)

Doodle samples 2

Let’s set dep to .5:

smples2 <- doodle_samples(n_colleagues, o, p, dep = .5)
glimpse(smples2)
#> Rows: 1,300
#> Columns: 4
#> $ slots_chosen <fct> 4, NA, 6, 7, 1, 6, 8, 10, 5, 9, NA, 1, 4, 7, NA, NA, 2, 7…
#> $ Freq         <int> 3, NA, 3, 3, 3, 3, 3, 3, 3, 3, NA, 3, 3, 3, NA, NA, 3, 3,…
#> $ is_match     <lgl> TRUE, FALSE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TRUE, TR…
#> $ id           <dbl> 1, 2, 3, 3, 4, 5, 5, 5, 6, 6, 7, 8, 8, 9, 10, 11, 12, 12,…

A similar dataset is available as data:

data(smples2)

Summarize samples

summarize_smples is what is finally called in prob_doodle_match.

summarize_samples(smples1)
#> # A tibble: 2 × 3
#>    Freq n_matches prop_matches
#>   <int>     <int>        <dbl>
#> 1     3       276        0.271
#> 2    NA       744        0.729