Quarantine Sexual Activity

library(ggplot2)
library(gganimate)
library(viridis)
library(hrbrthemes)

days = rep(seq(1, 21), 2)
f.o = c(3, 6, 5, 5, 0, 0, 0, 0, 0, 0, 0, 5, 7, 0, 7, 5, 2, 0, 11, 18, 2)
m.o = c(0, 0, 0, 1, 2, 0, 1, 0, 1, 0, 2, 0, 1, 1, 2, 1, 2, 0, 1, 0, 1)
gender = c(rep("Female", 21), rep("Male", 21))

q.o = data.frame(days, o = c(f.o, m.o), gender)

q.o %>%
  ggplot(aes(x=days, y = o, colour = gender, size = gender)) +
    geom_line() +
    scale_color_manual(values = c("pink", "blue")) +
    scale_size_manual(values = c(2, 2)) +
    ggtitle("Number of Orgasms Between Partners During 3 Weeks of Quarantine") +
    theme_ipsum() +
    ylab("Number of orgasms")