Circumcision

Is there a relationship between circumcision and erectile dysfunction?

Jonatan Pallesen
02-03-2018

Is there a relationship between circumcision and erectile dysfunction? Inspired by Seth Stephens-Davidowitz I explore this question using Google search frequencies. Specifically I look at the search frequencies of “erectile dysfunction” by state, and look at the correlation with circumcision rates by state.

read data


library(pacman)

p_load(tidyverse, magrittr, pander, janitor)

source('../../src/extra.R', echo = F, encoding="utf-8")



circumcision <- read_csv("data/circumcision_rates.csv")

ed <- read_csv("data/erectile_searches.csv")

df <- inner_join(circumcision, ed, by = "state")

plot


df %>% plot_lm("erectile_dysfunction_searches", "circumcision_rate", alpha = 1)


We see quite a strong correlation in this data. However, such state level analysis can be confounded in a multitude of ways. I won’t draw any conclusions from this, but I think it’s an interesting exploratory analysis.