Second R script