ggplot(na.omit(d), aes(x=log10(ncarehomes),y=covid,col=RegionNation)) +geom_point() +geom_smooth(method="lm") + facet_wrap(~RegionNation)
## `geom_smooth()` using formula 'y ~ x'