added analysis.R file