``` r d <- diamonds[sample(nrow(diamonds), 1000), ] p <- qplot(carat, price, colour=cut, data=d) + facet_wrap(~cut) + stat_smooth(method=lm) p ```  Here is the [plotly result](https://plot.ly/~cpsievert/246/fair-good-very-good-premium-ideal-fair-good-very-good-premium-ideal/) 