P-value bowling

[This article was first published on R on Publishable Stuff, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

I recently went bowling, and you know those weird 3D-animated bowling animations that all bowling alleys seemed to show whenever you made a strike? They are still alive and well! (At least at my local bowling place). And then I thought: Can I get animations like that into my daily data science workflow? With Rstudio’s built-in Viewer tab, I absolutely could! Below you find the code for a much improved t.test function that gives you different animations when you hit a strike ($p < 0.01$), a spare ($p < 0.05$), a “near miss” ($p < 0.1$) and a complete miss ($p > 0.1$).

(If you think this is silly, then I agree. Roughly as silly as using ritualized p-value cutoffs to decide whether an experiment is a “success” or not.)

The bowling.t.test function

# The animations we're going to use
# Source: https://archive.org/details/brunswick-frameworx-ten-pin-animations/
bowling_animation_urls <- c(
  "0.01" = "https://i.imgur.com/Kn8CQbj.gif",
  "0.05" = "https://i.imgur.com/HFTKdDL.gif",
  "0.1" = "https://i.imgur.com/8Sw54Mz.gif",
  "1" = "https://i.imgur.com/kjROdhj.gif"

bowling_animation_html <- sapply(bowling_animation_urls, function(url) {
  # We need to download the animations to the R session's temp directory  
  # to display them in the Rstudio Viewer pane
  animation_path = tempfile(pattern = "bowling_animation", fileext = ".gif")
  download.file(url, animation_path)
  # And we need to create an HTML wrapper around each GIF
  html_path = tempfile(pattern = "bowling_animation", fileext = ".html")
  html <- paste0(
    '<html><body><img src="',
    basename(animation_path) , 
    '" style="width: 100%;height: auto;"/></body></html>'
  writeLines(html, html_path)

# A much improved t.test function
bowling.t.test <- function(...) {
  test_result <- t.test(...)
  selected_animation_html <- bowling_animation_html[
    test_result$p.value < as.numeric(names(bowling_animation_urls)) 

And now trying it out!

y1 = rnorm(30, mean = 2)
y2 = rnorm(30)
# A successful experiment, we're great scientists.
bowling.t.test(y1, y2)

y1 = rnorm(30, mean = 0.7)
y2 = rnorm(30)
# Oh, that was close, but we're still great scientists.
bowling.t.test(y1, y2)

y1 = rnorm(30, mean = 0.65)
y2 = rnorm(30)
# What bad luck! Let's say it's trending towards significance maybe?
bowling.t.test(y1, y2)

y1 = rnorm(30, mean = 0.5)
y2 = rnorm(30)
# What are we even doing with our lives...
bowling.t.test(y1, y2)
To leave a comment for the author, please follow the link and comment on their blog: R on Publishable Stuff.

R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

Never miss an update!
Subscribe to R-bloggers to receive
e-mails with the latest R posts.
(You will not see this message again.)

Click here to close (This popup will not appear again)