Run analysis.r