add script to reproduce the entire analysis