Create analysis.r