update analyze.R