add R scripts and plot result