remove cobra.stats

The code was not referenced anywhere else in cobrapy. It provided a set
of convenient wrappers to statistical functions like p value adjustment
in rpy. However, that functionality can now be replaced by native python
implementations like statsmodels
doc update: phenotype phase plane

Added an example notebook for phenotype phase plane

Also fixed a few docstrings in core.