R script and results