added first R script and results