added analysis.R