R code for plog2
picture for plot2