analysis.R