All four plots + R code