ALT family_peds %>% ggplot(aes(x = birth_year_scaled, y = y)) + geom_point(aes( fill = post_1776, color = fam, shape = sex), alpha = .3) + geom_mark_hull( aes( x = birth_year_scaled, y = y, group = interaction(fam, gen), fill = NULL, color = fam ), concavity = 5, expand = unit(3, "mm"), radius = unit(1, "mm"), alpha = .3, fill = NA, linewidth = 0.6 ) + theme_bw() + labs(title = "Simulated Phenotypes by Family", x = "Scaled Birth Year", y = "Phenotype (y)") + # add viridis color scale_color_viridis_d(option = "plasma") + scale_fill_viridis_c()+ scale_shape_manual(values = c("F" = 21, "M" = 24)) + # remove color legend theme(legend.position = "bottom") + guides(color = "none", shape = guide_legend(override.aes = list(size = 3, alpha = 1, fill = "black")))