Better summary functions

We should have a summary function that mirrors the arviz summary object
https://python.arviz.org/en/latest/api/generated/arviz.summary.html

I'm not sure if this is directly applicable out of the box (since our objects are a bit higher dimensional than the examples).  