add run_analysis.R script