Create analysis.R