added r script to plot area