added analysis code and figure