Example code from MIT microbiota time-series