This is the first R script in this project