Is there an efficient way to get all of my linear models parameters distributions mean? I.e. analogue to sklearn's .coef_?

Questions
2

The summary function got moved to Arviz. You can still use it from there - see the code below.

import pymc3 as pm
import arviz as az

with pm.Model() as model:
    x = pm.Normal('x')
    y = pm.Deterministic('y', x**2)
    trace = pm.sample()
    
az.summary(trace)
1 Like