R Script for above charts