Make returned object hierarchical to give names to all objects

[billdenney]

Fix #8

The return value of the main plot() routine now returns a named list.

Also this PR speeds up testing by estimating fewer model iterations.