add R script for paper