updated analysis.R