Sum of distributions

give the following PyMC3 code:

import pymc3 as pm

with pm.Model() as model:
    a = pm.Normal('a', 3, 10)
    b = pm.Normal('b', 10, 5)
    c = a+b

How the c distribution is computed?

1 Like

Technically speaking, c is a random variable, which follow a distribution that is the result of the convolution between a and b.

In practice, PyMC3 samples from a and b (a single realization of the random variable), and add the number together