Is there a relationship between circumcision and erectile dysfunction?

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.

```
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")
```

```
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.