Hello,

I’m new to Bayesian modelling so forgive if this question is obvious. I initialize a model, `m1`

, and an identical model `m2`

, as

```
samples = scipy.stats.norm.rvs(loc=1,sd=1,size=10000)
with pm.Model() as m1:
mu=pm.Normal("mu",mu=0, sd=1)
std = pm.Gamma("std",mu=0.5,sd=3)
output=pm.Normal("output",mu=mu,sd=std, observed = samples)
```

However, when I do

```
logp = m1.logp
logps = [logp(trace[i]) for i in range(len(trace))]
```

after training `m1`

, and comparing against:

```
logp_2 = m2.logp
logps_2 = [logp_2(trace[i]) for i in range(len(trace))]
```

I get the same answers, even though `m1`

has been trained but `m2`

has not. Can someone please explain why this is?

Thank you!