differential_abundance

Now written with Rscript instead of rpy2.  tests also added.
replace script name with ``%prog``