I have a normally distributed parameter and pm.summary() reports a smaller sd then when I use trace.posterior[‘’].std(axis=(0, 1)) (2 chains). Does pm.summary() not pool across chains when calculating SD?

Hello, thanks for you question. Could you post a code snippet that reproduces your problem? For example, I tried: with pm.Model(): pm.Normal('m') trace = pm.sample() pm.summary(trace) and trace.posterior['m'].std(axis=(0, 1)) But I found identical results.

Different SD between pm.summary() and trace.posterior.std()

Questions version agnostic

daniel-saunders-phil September 6, 2024, 1:10pm 4

ah no worries - glad you sorted it out.

Topic		Replies	Views
Bayesian regression in PyMC3 Questions	1	399	October 16, 2019
How to obtain the variance estimations in Linear Regression models Questions	5	889	June 26, 2017
Compute reduced summary statistics Questions	3	1612	April 30, 2020
Prevent roundoff on pm.summary Questions	1	442	May 5, 2020
Getting statistics from a single chain Questions	1	387	January 3, 2020

Different SD between pm.summary() and trace.posterior.std()

Related topics