Finished new analysis script

Run from command line with Rscript; work on knitr report…