Removed geometric dist from function.  Will add back in later after troubleshooting. Made output abundances and lm on log10 scale.