I was tasked to create this plot twinData %>% ggplot(mapping = aes(ht1, ht2)) + geom_point() + facet_wrap(cohort ~ zygosity) this is the resulting plot
twinData %>% ggplot(mapping = aes(ht1, ht2)) + geom_point() + facet_wrap(cohort ~ zygosity)