Added new analysis script
brought back two last paragaraphs in ms.Rmd. #188da13