Add R script and output for 4th plot (four panels)