grouped_ggscatterstats function

Scatterplot with marginal distributions for all levels of a grouping variable