modified:   analysis.R