+ `infer_independent` was refactored: It is now implemented in `CalibrationModel` and calculates both ETI and HDI by default.
+ New method `CalibrationModel.likelihood` that can be parametrized with `scan_x=True` to evaluate the likelihood at all entries in `x`, instead of doing elementwise evaluation along `x` and `y`.
+ Imports work with `theano-pymc v1.1.2` (used by PyMC3 `v3.11.1`) and `aesara v2`
+ Added option to plot relative residuals in `plot_model` and `plot_t_band`
+ Guessing functions were removed because they are not very useful in practice
+ License document (and setup.py) were updated.