Statistics

See in #22547

MatplotlibDeprecationWarning: Axes3D(fig) adding itself to the figure is deprecated since 3.4. Pass the keyword argument auto_add_to_figure=False and use fig.add_axes(ax) to suppress this warning. The default value of auto_add_to_figure will change to False in mpl3.5 and True values will no longer work in 3.6.  This is consistent with other Axes classes

We need to rep

Collection of follow-ups to #5827. These can/should be broken out into individual PRs. Many are relatively straightforward and would make a good first PR.

General

• Documentation (none was added in original PR).
• Release notes.
• Example notebook.
• Double-check how sm.tsa.arima.ARIMA works with fix_params (it should fail except when the fit method is statespace

See:
gonum/gonum#845 (comment)

Are there any plans to add a Zero-Inflated Poisson (ZIP) and Zero-Inflated Negative Binomial (ZINB) to TFP? Those are usually very common distributions in other packages, and it shouldn't be hard to implement.

Since the default output is meant to be human-readable, would it make sense to add thousands separators to make the output more easily readable?