first commit

Michele's R script from Nov 2014 meetup