add run_analysis.r script