Add run_analysis.R script