Exposed p-values and other summary stats.

Added an option to return the df in summary() instead of printing. Default behavior is to print a summary (unchanged).
Should this be added to other fitters as well? Or is there a better way to do this? Perhaps a "_pvalue" property set by fit()?