Attach trial wise parameters to idata¶

HSSM automatically cleans up the idata object we return from samplers, to avoid unnecessarily huge objects from being passed around (the decision was taken after having observed many inadvertent out of memory errors).

Instead, if so desired, we can attach the trial wise parameters to the idata object ourselves, using the add_likelihood_parameters_to_idata function.

This quick tutorial shows you how to do this.

Load Modules¶

In [1]:

import pandas as pd
import hssm

import pandas as pd import hssm

Simulate Data¶

In [2]:

# Condition 1
condition_1 = hssm.simulate_data(
    model="ddm", theta=dict(v=0.5, a=1.5, z=0.5, t=0.1), size=500
)

# Condition 2
condition_2 = hssm.simulate_data(
    model="ddm", theta=dict(v=1.0, a=1.5, z=0.5, t=0.1), size=500
)

condition_1["condition"] = "C1"
condition_2["condition"] = "C2"

data = pd.concat([condition_1, condition_2]).reset_index(drop=True)

# Condition 1 condition_1 = hssm.simulate_data( model="ddm", theta=dict(v=0.5, a=1.5, z=0.5, t=0.1), size=500 ) # Condition 2 condition_2 = hssm.simulate_data( model="ddm", theta=dict(v=1.0, a=1.5, z=0.5, t=0.1), size=500 ) condition_1["condition"] = "C1" condition_2["condition"] = "C2" data = pd.concat([condition_1, condition_2]).reset_index(drop=True)

Define Model and Sample¶

In [3]:

model = hssm.HSSM(
    model="ddm",
    data=data,
    include=[
        {
            "name": "v",
            "formula": "v ~ 1 + condition",
        }
    ],
)

idata = model.sample(sampler="mcmc", tune=500, draws=500)

model = hssm.HSSM( model="ddm", data=data, include=[ { "name": "v", "formula": "v ~ 1 + condition", } ], ) idata = model.sample(sampler="mcmc", tune=500, draws=500)

Model initialized successfully.
Using default initvals. 

Initializing NUTS using adapt_diag...
Multiprocess sampling (4 chains in 4 jobs)
NUTS: [a, t, z, v_Intercept, v_condition]

Output()

Sampling 4 chains for 500 tune and 500 draw iterations (2_000 + 2_000 draws total) took 12 seconds.
/Users/afengler/Library/CloudStorage/OneDrive-Personal/proj_hssm/HSSM/.venv/lib/python3.12/site-packages/pymc/pytensorf.py:958: FutureWarning: compile_pymc was renamed to compile. Old name will be removed in a future release of PyMC
  warnings.warn(
100%|██████████| 2000/2000 [00:00<00:00, 2468.61it/s]

Checking the idata object¶

In [4]:

idata

idata

Out[4]:

arviz.InferenceData

• posterior

  <xarray.Dataset> Size: 84kB
  Dimensions:          (chain: 4, draw: 500, v_condition_dim: 1)
  Coordinates:
    * chain            (chain) int64 32B 0 1 2 3
    * draw             (draw) int64 4kB 0 1 2 3 4 5 6 ... 494 495 496 497 498 499
    * v_condition_dim  (v_condition_dim) <U2 8B 'C2'
  Data variables:
      v_Intercept      (chain, draw) float64 16kB 0.5171 0.54 ... 0.6887 0.6311
      v_condition      (chain, draw, v_condition_dim) float64 16kB 0.54 ... 0.4947
      z                (chain, draw) float64 16kB 0.511 0.506 ... 0.4611 0.4504
      a                (chain, draw) float64 16kB 1.474 1.45 1.437 ... 1.499 1.481
      t                (chain, draw) float64 16kB 0.157 0.1533 ... 0.1124 0.09461
  Attributes:
      created_at:                  2025-07-13T12:29:09.083594+00:00
      arviz_version:               0.21.0
      inference_library:           pymc
      inference_library_version:   5.21.1
      sampling_time:               12.336488962173462
      tuning_steps:                500
      modeling_interface:          bambi
      modeling_interface_version:  0.15.0

  xarray.Dataset

  • Dimensions:
    • chain: 4
    • draw: 500
    • v_condition_dim: 1
  • Coordinates: (3)
    • chain
      (chain)
      int64
      0 1 2 3

      array([0, 1, 2, 3])

    • draw
      (draw)
      int64
      0 1 2 3 4 5 ... 495 496 497 498 499

      array([  0,   1,   2, ..., 497, 498, 499])

    • v_condition_dim
      (v_condition_dim)
      <U2
      'C2'

      array(['C2'], dtype='<U2')

  • Data variables: (5)
    • v_Intercept
      (chain, draw)
      float64
      0.5171 0.54 ... 0.6887 0.6311

      array([[0.51714871, 0.54004777, 0.62149724, ..., 0.61966065, 0.59243802,
              0.59243802],
             [0.61869938, 0.65723774, 0.63021498, ..., 0.6533621 , 0.6533621 ,
              0.63131509],
             [0.59528702, 0.59528702, 0.66226081, ..., 0.65958662, 0.66126487,
              0.62014973],
             [0.69280671, 0.60980159, 0.64506599, ..., 0.5786166 , 0.68870133,
              0.63114493]])

    • v_condition
      (chain, draw, v_condition_dim)
      float64
      0.54 0.4807 ... 0.4821 0.4947

      array([[[0.54001922],
              [0.48071279],
              [0.40194477],
              ...,
              [0.40263069],
              [0.574422  ],
              [0.574422  ]],
      
             [[0.41762479],
              [0.40570275],
              [0.49960274],
              ...,
              [0.36916889],
              [0.36916889],
              [0.52625447]],
      
             [[0.48220378],
              [0.48220378],
              [0.44222329],
              ...,
              [0.41889775],
              [0.41845161],
              [0.39375802]],
      
             [[0.51630056],
              [0.47777804],
              [0.46842409],
              ...,
              [0.52539226],
              [0.48205419],
              [0.49468106]]])

    • z
      (chain, draw)
      float64
      0.511 0.506 ... 0.4611 0.4504

      array([[0.51100413, 0.50596584, 0.47633051, ..., 0.47176446, 0.46930924,
              0.46930924],
             [0.50300194, 0.51073864, 0.44209232, ..., 0.4921373 , 0.4921373 ,
              0.46790296],
             [0.48919294, 0.48919294, 0.48041059, ..., 0.48003018, 0.46176087,
              0.48517901],
             [0.47304107, 0.46580817, 0.4675487 , ..., 0.47526074, 0.46114424,
              0.45035494]])

    • a
      (chain, draw)
      float64
      1.474 1.45 1.437 ... 1.499 1.481

      array([[1.47383317, 1.44986107, 1.43683483, ..., 1.47438526, 1.45438234,
              1.45438234],
             [1.46982248, 1.44986046, 1.53209259, ..., 1.44006039, 1.44006039,
              1.5302841 ],
             [1.44286499, 1.44286499, 1.51295776, ..., 1.51369058, 1.53867973,
              1.55879449],
             [1.4633992 , 1.54405247, 1.43534039, ..., 1.48265867, 1.49921362,
              1.4806146 ]])

    • t
      (chain, draw)
      float64
      0.157 0.1533 ... 0.1124 0.09461

      array([[0.15697393, 0.15326962, 0.13602127, ..., 0.12049825, 0.14573151,
              0.14573151],
             [0.14721884, 0.13517113, 0.05917854, ..., 0.13224046, 0.13224046,
              0.09130519],
             [0.15908614, 0.15908614, 0.10971873, ..., 0.07846036, 0.07129759,
              0.05361361],
             [0.13873064, 0.06621154, 0.12420361, ..., 0.10721722, 0.11242791,
              0.09461212]])

  • Indexes: (3)
    • chain
      PandasIndex

      PandasIndex(Index([0, 1, 2, 3], dtype='int64', name='chain'))

    • draw
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             490, 491, 492, 493, 494, 495, 496, 497, 498, 499],
            dtype='int64', name='draw', length=500))

    • v_condition_dim
      PandasIndex

      PandasIndex(Index(['C2'], dtype='object', name='v_condition_dim'))

  • Attributes: (8)

    created_at :

    2025-07-13T12:29:09.083594+00:00

    arviz_version :

    0.21.0

    inference_library :

    pymc

    inference_library_version :

    5.21.1

    sampling_time :

    12.336488962173462

    tuning_steps :

    500

    modeling_interface :

    bambi

    modeling_interface_version :

    0.15.0

  
  
• log_likelihood

  <xarray.Dataset> Size: 16MB
  Dimensions:      (chain: 4, draw: 500, __obs__: 1000)
  Coordinates:
    * chain        (chain) int64 32B 0 1 2 3
    * draw         (draw) int64 4kB 0 1 2 3 4 5 6 ... 493 494 495 496 497 498 499
    * __obs__      (__obs__) int64 8kB 0 1 2 3 4 5 6 ... 994 995 996 997 998 999
  Data variables:
      rt,response  (chain, draw, __obs__) float64 16MB -0.8921 -1.164 ... -1.179
  Attributes:
      modeling_interface:          bambi
      modeling_interface_version:  0.15.0

  xarray.Dataset

  • Dimensions:
    • chain: 4
    • draw: 500
    • __obs__: 1000
  • Coordinates: (3)
    • chain
      (chain)
      int64
      0 1 2 3

      array([0, 1, 2, 3])

    • draw
      (draw)
      int64
      0 1 2 3 4 5 ... 495 496 497 498 499

      array([  0,   1,   2, ..., 497, 498, 499])

    • __obs__
      (__obs__)
      int64
      0 1 2 3 4 5 ... 995 996 997 998 999

      array([  0,   1,   2, ..., 997, 998, 999])

  • Data variables: (1)
    • rt,response
      (chain, draw, __obs__)
      float64
      -0.8921 -1.164 ... -2.164 -1.179

      array([[[-0.89211002, -1.16416096, -1.74157262, ..., -1.29116593,
               -2.19774012, -1.23217246],
              [-0.85939804, -1.10405186, -1.7483828 , ..., -1.2982562 ,
               -2.18674077, -1.24056155],
              [-0.83830892, -1.12891851, -1.71366241, ..., -1.26750297,
               -2.15073771, -1.21118905],
              ...,
              [-0.88827243, -1.21153371, -1.68716591, ..., -1.2422521 ,
               -2.09630508, -1.18816096],
              [-0.90950666, -1.29416587, -1.69070797, ..., -1.25426613,
               -2.24539089, -1.19126137],
              [-0.90950666, -1.29416587, -1.69070797, ..., -1.25426613,
               -2.24539089, -1.19126137]],
      
             [[-0.80220639, -1.06013377, -1.71239699, ..., -1.28337643,
               -2.1702024 , -1.22596392],
              [-0.73239504, -0.89093473, -1.74686954, ..., -1.33861634,
               -2.26417379, -1.27827023],
              [-1.01148242, -1.37209535, -1.65049286, ..., -1.21337361,
               -2.10477268, -1.15806577],
      ...
              [-0.85592306, -1.08739794, -1.68707219, ..., -1.27703432,
               -2.15854744, -1.22085344],
              [-0.92577891, -1.24395907, -1.64858187, ..., -1.22610385,
               -2.08160689, -1.17246637],
              [-0.92436085, -1.14660565, -1.6854833 , ..., -1.26721903,
               -2.07229792, -1.21598761]],
      
             [[-0.79648859, -1.13895111, -1.67101372, ..., -1.27627423,
               -2.30303456, -1.2107856 ],
              [-0.97835803, -1.28437205, -1.66772511, ..., -1.23601603,
               -2.09786856, -1.1818053 ],
              [-0.83382596, -1.12131307, -1.70876056, ..., -1.28053028,
               -2.23665188, -1.21970116],
              ...,
              [-0.93128737, -1.22248361, -1.70623697, ..., -1.26856145,
               -2.18840754, -1.21000046],
              [-0.86543486, -1.23657433, -1.6442152 , ..., -1.23265566,
               -2.192903  , -1.17212604],
              [-0.9384307 , -1.28507594, -1.67427846, ..., -1.23724561,
               -2.16361154, -1.17916162]]])

  • Indexes: (3)
    • chain
      PandasIndex

      PandasIndex(Index([0, 1, 2, 3], dtype='int64', name='chain'))

    • draw
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             490, 491, 492, 493, 494, 495, 496, 497, 498, 499],
            dtype='int64', name='draw', length=500))

    • __obs__
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             990, 991, 992, 993, 994, 995, 996, 997, 998, 999],
            dtype='int64', name='__obs__', length=1000))

  • Attributes: (2)

    modeling_interface :

    bambi

    modeling_interface_version :

    0.15.0

  
  
• sample_stats

  <xarray.Dataset> Size: 248kB
  Dimensions:                (chain: 4, draw: 500)
  Coordinates:
    * chain                  (chain) int64 32B 0 1 2 3
    * draw                   (draw) int64 4kB 0 1 2 3 4 5 ... 495 496 497 498 499
  Data variables: (12/17)
      smallest_eigval        (chain, draw) float64 16kB nan nan nan ... nan nan
      acceptance_rate        (chain, draw) float64 16kB 0.999 0.9843 ... 0.9941
      perf_counter_diff      (chain, draw) float64 16kB 0.007624 ... 0.008489
      step_size              (chain, draw) float64 16kB 0.5975 0.5975 ... 0.6306
      energy                 (chain, draw) float64 16kB 1.631e+03 ... 1.626e+03
      step_size_bar          (chain, draw) float64 16kB 0.5078 0.5078 ... 0.4455
      ...                     ...
      tree_depth             (chain, draw) int64 16kB 3 3 4 4 3 3 ... 4 3 3 4 3 3
      reached_max_treedepth  (chain, draw) bool 2kB False False ... False False
      perf_counter_start     (chain, draw) float64 16kB 1.427e+06 ... 1.427e+06
      largest_eigval         (chain, draw) float64 16kB nan nan nan ... nan nan
      diverging              (chain, draw) bool 2kB False False ... False False
      max_energy_error       (chain, draw) float64 16kB -0.97 -0.2514 ... -0.0425
  Attributes:
      created_at:                  2025-07-13T12:29:09.093660+00:00
      arviz_version:               0.21.0
      inference_library:           pymc
      inference_library_version:   5.21.1
      sampling_time:               12.336488962173462
      tuning_steps:                500
      modeling_interface:          bambi
      modeling_interface_version:  0.15.0

  xarray.Dataset

  • Dimensions:
    • chain: 4
    • draw: 500
  • Coordinates: (2)
    • chain
      (chain)
      int64
      0 1 2 3

      array([0, 1, 2, 3])

    • draw
      (draw)
      int64
      0 1 2 3 4 5 ... 495 496 497 498 499

      array([  0,   1,   2, ..., 497, 498, 499])

  • Data variables: (17)
    • smallest_eigval
      (chain, draw)
      float64
      nan nan nan nan ... nan nan nan nan

      array([[nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan]])

    • acceptance_rate
      (chain, draw)
      float64
      0.999 0.9843 ... 0.8754 0.9941

      array([[0.99899574, 0.98427376, 0.94007693, ..., 0.88913011, 0.38067985,
              0.17382386],
             [0.96963385, 0.35093074, 0.87165549, ..., 0.95839737, 0.56263805,
              0.98164271],
             [0.82402644, 0.35254835, 0.912601  , ..., 0.96508482, 0.91995735,
              0.86907564],
             [0.96060646, 1.        , 0.89623832, ..., 1.        , 0.87539908,
              0.99405597]])

    • perf_counter_diff
      (chain, draw)
      float64
      0.007624 0.007353 ... 0.008489

      array([[0.007624  , 0.00735321, 0.01615529, ..., 0.00753525, 0.00735346,
              0.00777696],
             [0.00738925, 0.00377908, 0.0084235 , ..., 0.00896825, 0.00190383,
              0.01592554],
             [0.00743321, 0.00786996, 0.00793646, ..., 0.00747246, 0.00717438,
              0.00358404],
             [0.00732617, 0.00734325, 0.00763408, ..., 0.01799921, 0.00955742,
              0.00848917]])

    • step_size
      (chain, draw)
      float64
      0.5975 0.5975 ... 0.6306 0.6306

      array([[0.59752358, 0.59752358, 0.59752358, ..., 0.59752358, 0.59752358,
              0.59752358],
             [0.38784754, 0.38784754, 0.38784754, ..., 0.38784754, 0.38784754,
              0.38784754],
             [0.49690078, 0.49690078, 0.49690078, ..., 0.49690078, 0.49690078,
              0.49690078],
             [0.63061791, 0.63061791, 0.63061791, ..., 0.63061791, 0.63061791,
              0.63061791]])

    • energy
      (chain, draw)
      float64
      1.631e+03 1.63e+03 ... 1.626e+03

      array([[1630.69631792, 1629.62734758, 1627.05643916, ..., 1625.80227866,
              1629.66579447, 1629.7756048 ],
             [1626.32647922, 1632.48370731, 1634.83004223, ..., 1629.90944826,
              1628.46312271, 1628.05114814],
             [1628.45969076, 1631.92495563, 1629.71887201, ..., 1627.93737279,
              1628.60177843, 1632.05979479],
             [1630.82339503, 1631.66538563, 1628.52456439, ..., 1627.95922118,
              1627.68811079, 1626.27266518]])

    • step_size_bar
      (chain, draw)
      float64
      0.5078 0.5078 ... 0.4455 0.4455

      array([[0.50781584, 0.50781584, 0.50781584, ..., 0.50781584, 0.50781584,
              0.50781584],
             [0.49008555, 0.49008555, 0.49008555, ..., 0.49008555, 0.49008555,
              0.49008555],
             [0.50400304, 0.50400304, 0.50400304, ..., 0.50400304, 0.50400304,
              0.50400304],
             [0.4454689 , 0.4454689 , 0.4454689 , ..., 0.4454689 , 0.4454689 ,
              0.4454689 ]])

    • index_in_trajectory
      (chain, draw)
      int64
      3 -1 5 3 -1 -4 ... 5 -1 5 -10 -5 4

      array([[  3,  -1,   5, ...,   4,  -4,   0],
             [  4,   1,  -6, ...,   6,   0,   7],
             [ -5,   0,   2, ...,  -4,   2,  -3],
             [  4,  -4,   7, ..., -10,  -5,   4]])

    • n_steps
      (chain, draw)
      float64
      7.0 7.0 15.0 15.0 ... 15.0 7.0 7.0

      array([[ 7.,  7., 15., ...,  7.,  7.,  7.],
             [ 7.,  3.,  7., ...,  7.,  1., 15.],
             [ 7.,  7.,  7., ...,  7.,  7.,  3.],
             [ 7.,  7.,  7., ..., 15.,  7.,  7.]])

    • process_time_diff
      (chain, draw)
      float64
      0.007551 0.007321 ... 0.008173

      array([[0.007551, 0.007321, 0.015356, ..., 0.00746 , 0.00733 , 0.007642],
             [0.007381, 0.003748, 0.007703, ..., 0.007894, 0.001904, 0.01537 ],
             [0.007412, 0.007561, 0.00774 , ..., 0.007458, 0.007158, 0.003584],
             [0.007253, 0.007335, 0.007509, ..., 0.016259, 0.008841, 0.008173]])

    • energy_error
      (chain, draw)
      float64
      -0.9617 -0.2514 ... 0.043 0.01939

      array([[-0.96167665, -0.25139885,  0.01219149, ...,  0.21640905,
              -0.20854996,  0.        ],
             [ 0.02708104,  2.59162301, -1.47941526, ..., -0.01433663,
               0.        , -0.33364469],
             [-0.20843983,  0.        , -0.01997814, ...,  0.01511287,
              -0.50103452,  0.16767406],
             [ 0.12858918, -0.36713549, -0.113043  , ..., -0.66845625,
               0.04299687,  0.01938668]])

    • lp
      (chain, draw)
      float64
      -1.627e+03 ... -1.625e+03

      array([[-1627.47078996, -1625.58003051, -1624.4655962 , ...,
              -1624.72595789, -1627.00291233, -1627.00291233],
             [-1625.22588946, -1632.01756713, -1629.07734271, ...,
              -1626.44270871, -1626.44270871, -1625.37626727],
             [-1625.12895309, -1625.12895309, -1624.4901249 , ...,
              -1627.14196366, -1626.631181  , -1631.73134296],
             [-1628.98405346, -1627.41478378, -1624.9096195 , ...,
              -1624.39036309, -1625.23857872, -1625.28751373]])

    • tree_depth
      (chain, draw)
      int64
      3 3 4 4 3 3 3 3 ... 3 3 4 3 3 4 3 3

      array([[3, 3, 4, ..., 3, 3, 3],
             [3, 2, 3, ..., 3, 1, 4],
             [3, 3, 3, ..., 3, 3, 2],
             [3, 3, 3, ..., 4, 3, 3]])

    • reached_max_treedepth
      (chain, draw)
      bool
      False False False ... False False

      array([[False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False]])

    • perf_counter_start
      (chain, draw)
      float64
      1.427e+06 1.427e+06 ... 1.427e+06

      array([[1427401.97826979, 1427401.98597425, 1427401.99338738, ...,
              1427406.12062838, 1427406.12823167, 1427406.13566033],
             [1427402.13910204, 1427402.14656562, 1427402.15040271, ...,
              1427406.61256283, 1427406.62165175, 1427406.624541  ],
             [1427402.14002608, 1427402.14753121, 1427402.15554638, ...,
              1427406.53236171, 1427406.53989321, 1427406.54712004],
             [1427402.21389933, 1427402.22129496, 1427402.22872138, ...,
              1427406.63153483, 1427406.64973854, 1427406.65941208]])

    • largest_eigval
      (chain, draw)
      float64
      nan nan nan nan ... nan nan nan nan

      array([[nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan]])

    • diverging
      (chain, draw)
      bool
      False False False ... False False

      array([[False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False]])

    • max_energy_error
      (chain, draw)
      float64
      -0.97 -0.2514 ... 0.2286 -0.0425

      array([[-0.96996264, -0.25139885,  0.19062238, ...,  0.27480448,
               2.12027768,  2.6128521 ],
             [ 0.08181599,  2.7727651 ,  2.22317105, ..., -0.81901335,
               0.57511875, -0.66149511],
             [ 0.47257644,  2.50204722,  0.33240975, ..., -0.91055443,
              -1.09468036,  0.21075921],
             [-0.30839104, -0.59334298,  0.33479679, ..., -0.66845625,
               0.22861677, -0.04249524]])

  • Indexes: (2)
    • chain
      PandasIndex

      PandasIndex(Index([0, 1, 2, 3], dtype='int64', name='chain'))

    • draw
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             490, 491, 492, 493, 494, 495, 496, 497, 498, 499],
            dtype='int64', name='draw', length=500))

  • Attributes: (8)

    created_at :

    2025-07-13T12:29:09.093660+00:00

    arviz_version :

    0.21.0

    inference_library :

    pymc

    inference_library_version :

    5.21.1

    sampling_time :

    12.336488962173462

    tuning_steps :

    500

    modeling_interface :

    bambi

    modeling_interface_version :

    0.15.0

  
  
• observed_data

  <xarray.Dataset> Size: 24kB
  Dimensions:                  (__obs__: 1000, rt,response_extra_dim_0: 2)
  Coordinates:
    * __obs__                  (__obs__) int64 8kB 0 1 2 3 4 ... 996 997 998 999
    * rt,response_extra_dim_0  (rt,response_extra_dim_0) int64 16B 0 1
  Data variables:
      rt,response              (__obs__, rt,response_extra_dim_0) float64 16kB ...
  Attributes:
      created_at:                  2025-07-13T12:29:09.096399+00:00
      arviz_version:               0.21.0
      inference_library:           pymc
      inference_library_version:   5.21.1
      modeling_interface:          bambi
      modeling_interface_version:  0.15.0

  xarray.Dataset

  • Dimensions:
    • __obs__: 1000
    • rt,response_extra_dim_0: 2
  • Coordinates: (2)
    • __obs__
      (__obs__)
      int64
      0 1 2 3 4 5 ... 995 996 997 998 999

      array([  0,   1,   2, ..., 997, 998, 999])

    • rt,response_extra_dim_0
      (rt,response_extra_dim_0)
      int64
      0 1

      array([0, 1])

  • Data variables: (1)
    • rt,response
      (__obs__, rt,response_extra_dim_0)
      float64
      0.8199 1.0 0.5436 ... 1.0 1.668 1.0

      array([[0.81994694, 1.        ],
             [0.54364127, 1.        ],
             [2.20367026, 1.        ],
             ...,
             [1.72080851, 1.        ],
             [2.53323507, 1.        ],
             [1.66819727, 1.        ]])

  • Indexes: (2)
    • __obs__
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             990, 991, 992, 993, 994, 995, 996, 997, 998, 999],
            dtype='int64', name='__obs__', length=1000))

    • rt,response_extra_dim_0
      PandasIndex

      PandasIndex(Index([0, 1], dtype='int64', name='rt,response_extra_dim_0'))

  • Attributes: (6)

    created_at :

    2025-07-13T12:29:09.096399+00:00

    arviz_version :

    0.21.0

    inference_library :

    pymc

    inference_library_version :

    5.21.1

    modeling_interface :

    bambi

    modeling_interface_version :

    0.15.0

  
  

Note that none of the parameters attached in the posterior have the observation coordinate (__obs__). However, looking at the model graph, we clearly see that for variables that are represented as regression targets, we in fact compute intermediate deterministics which do have dimensionality __obs__.

In [5]:

model.graph()

model.graph()

Out[5]:

Computing Trial Wise Parameters¶

We can use the compute_likelihood_parameters_to_idata function to recompute and include in our idata the trial wise deterministics, which are part of our model graph.

In [6]:

### Attach Trial Wise Parameters
idata_trialwise = model.add_likelihood_parameters_to_idata(idata)

### Checking the `idata` object
idata_trialwise

### Attach Trial Wise Parameters idata_trialwise = model.add_likelihood_parameters_to_idata(idata) ### Checking the `idata` object idata_trialwise

Out[6]:

arviz.InferenceData

• posterior

  <xarray.Dataset> Size: 16MB
  Dimensions:          (chain: 4, draw: 500, v_condition_dim: 1, __obs__: 1000)
  Coordinates:
    * chain            (chain) int64 32B 0 1 2 3
    * draw             (draw) int64 4kB 0 1 2 3 4 5 6 ... 494 495 496 497 498 499
    * v_condition_dim  (v_condition_dim) <U2 8B 'C2'
    * __obs__          (__obs__) int64 8kB 0 1 2 3 4 5 ... 994 995 996 997 998 999
  Data variables:
      v_Intercept      (chain, draw) float64 16kB 0.5171 0.54 ... 0.6887 0.6311
      v_condition      (chain, draw, v_condition_dim) float64 16kB 0.54 ... 0.4947
      z                (chain, draw) float64 16kB 0.511 0.506 ... 0.4611 0.4504
      a                (chain, draw) float64 16kB 1.474 1.45 1.437 ... 1.499 1.481
      t                (chain, draw) float64 16kB 0.157 0.1533 ... 0.1124 0.09461
      v                (chain, draw, __obs__) float64 16MB 0.5171 0.5171 ... 1.126
  Attributes:
      created_at:                  2025-07-13T12:29:09.083594+00:00
      arviz_version:               0.21.0
      inference_library:           pymc
      inference_library_version:   5.21.1
      sampling_time:               12.336488962173462
      tuning_steps:                500
      modeling_interface:          bambi
      modeling_interface_version:  0.15.0

  xarray.Dataset

  • Dimensions:
    • chain: 4
    • draw: 500
    • v_condition_dim: 1
    • __obs__: 1000
  • Coordinates: (4)
    • chain
      (chain)
      int64
      0 1 2 3

      array([0, 1, 2, 3])

    • draw
      (draw)
      int64
      0 1 2 3 4 5 ... 495 496 497 498 499

      array([  0,   1,   2, ..., 497, 498, 499])

    • v_condition_dim
      (v_condition_dim)
      <U2
      'C2'

      array(['C2'], dtype='<U2')

    • __obs__
      (__obs__)
      int64
      0 1 2 3 4 5 ... 995 996 997 998 999

      array([  0,   1,   2, ..., 997, 998, 999])

  • Data variables: (6)
    • v_Intercept
      (chain, draw)
      float64
      0.5171 0.54 ... 0.6887 0.6311

      array([[0.51714871, 0.54004777, 0.62149724, ..., 0.61966065, 0.59243802,
              0.59243802],
             [0.61869938, 0.65723774, 0.63021498, ..., 0.6533621 , 0.6533621 ,
              0.63131509],
             [0.59528702, 0.59528702, 0.66226081, ..., 0.65958662, 0.66126487,
              0.62014973],
             [0.69280671, 0.60980159, 0.64506599, ..., 0.5786166 , 0.68870133,
              0.63114493]])

    • v_condition
      (chain, draw, v_condition_dim)
      float64
      0.54 0.4807 ... 0.4821 0.4947

      array([[[0.54001922],
              [0.48071279],
              [0.40194477],
              ...,
              [0.40263069],
              [0.574422  ],
              [0.574422  ]],
      
             [[0.41762479],
              [0.40570275],
              [0.49960274],
              ...,
              [0.36916889],
              [0.36916889],
              [0.52625447]],
      
             [[0.48220378],
              [0.48220378],
              [0.44222329],
              ...,
              [0.41889775],
              [0.41845161],
              [0.39375802]],
      
             [[0.51630056],
              [0.47777804],
              [0.46842409],
              ...,
              [0.52539226],
              [0.48205419],
              [0.49468106]]])

    • z
      (chain, draw)
      float64
      0.511 0.506 ... 0.4611 0.4504

      array([[0.51100413, 0.50596584, 0.47633051, ..., 0.47176446, 0.46930924,
              0.46930924],
             [0.50300194, 0.51073864, 0.44209232, ..., 0.4921373 , 0.4921373 ,
              0.46790296],
             [0.48919294, 0.48919294, 0.48041059, ..., 0.48003018, 0.46176087,
              0.48517901],
             [0.47304107, 0.46580817, 0.4675487 , ..., 0.47526074, 0.46114424,
              0.45035494]])

    • a
      (chain, draw)
      float64
      1.474 1.45 1.437 ... 1.499 1.481

      array([[1.47383317, 1.44986107, 1.43683483, ..., 1.47438526, 1.45438234,
              1.45438234],
             [1.46982248, 1.44986046, 1.53209259, ..., 1.44006039, 1.44006039,
              1.5302841 ],
             [1.44286499, 1.44286499, 1.51295776, ..., 1.51369058, 1.53867973,
              1.55879449],
             [1.4633992 , 1.54405247, 1.43534039, ..., 1.48265867, 1.49921362,
              1.4806146 ]])

    • t
      (chain, draw)
      float64
      0.157 0.1533 ... 0.1124 0.09461

      array([[0.15697393, 0.15326962, 0.13602127, ..., 0.12049825, 0.14573151,
              0.14573151],
             [0.14721884, 0.13517113, 0.05917854, ..., 0.13224046, 0.13224046,
              0.09130519],
             [0.15908614, 0.15908614, 0.10971873, ..., 0.07846036, 0.07129759,
              0.05361361],
             [0.13873064, 0.06621154, 0.12420361, ..., 0.10721722, 0.11242791,
              0.09461212]])

    • v
      (chain, draw, __obs__)
      float64
      0.5171 0.5171 ... 1.126 1.126

      array([[[0.51714871, 0.51714871, 0.51714871, ..., 1.05716793,
               1.05716793, 1.05716793],
              [0.54004777, 0.54004777, 0.54004777, ..., 1.02076057,
               1.02076057, 1.02076057],
              [0.62149724, 0.62149724, 0.62149724, ..., 1.02344201,
               1.02344201, 1.02344201],
              ...,
              [0.61966065, 0.61966065, 0.61966065, ..., 1.02229134,
               1.02229134, 1.02229134],
              [0.59243802, 0.59243802, 0.59243802, ..., 1.16686002,
               1.16686002, 1.16686002],
              [0.59243802, 0.59243802, 0.59243802, ..., 1.16686002,
               1.16686002, 1.16686002]],
      
             [[0.61869938, 0.61869938, 0.61869938, ..., 1.03632417,
               1.03632417, 1.03632417],
              [0.65723774, 0.65723774, 0.65723774, ..., 1.06294049,
               1.06294049, 1.06294049],
              [0.63021498, 0.63021498, 0.63021498, ..., 1.12981772,
               1.12981772, 1.12981772],
      ...
              [0.65958662, 0.65958662, 0.65958662, ..., 1.07848436,
               1.07848436, 1.07848436],
              [0.66126487, 0.66126487, 0.66126487, ..., 1.07971648,
               1.07971648, 1.07971648],
              [0.62014973, 0.62014973, 0.62014973, ..., 1.01390775,
               1.01390775, 1.01390775]],
      
             [[0.69280671, 0.69280671, 0.69280671, ..., 1.20910726,
               1.20910726, 1.20910726],
              [0.60980159, 0.60980159, 0.60980159, ..., 1.08757964,
               1.08757964, 1.08757964],
              [0.64506599, 0.64506599, 0.64506599, ..., 1.11349008,
               1.11349008, 1.11349008],
              ...,
              [0.5786166 , 0.5786166 , 0.5786166 , ..., 1.10400886,
               1.10400886, 1.10400886],
              [0.68870133, 0.68870133, 0.68870133, ..., 1.17075552,
               1.17075552, 1.17075552],
              [0.63114493, 0.63114493, 0.63114493, ..., 1.12582599,
               1.12582599, 1.12582599]]])

  • Indexes: (4)
    • chain
      PandasIndex

      PandasIndex(Index([0, 1, 2, 3], dtype='int64', name='chain'))

    • draw
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             490, 491, 492, 493, 494, 495, 496, 497, 498, 499],
            dtype='int64', name='draw', length=500))

    • v_condition_dim
      PandasIndex

      PandasIndex(Index(['C2'], dtype='object', name='v_condition_dim'))

    • __obs__
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             990, 991, 992, 993, 994, 995, 996, 997, 998, 999],
            dtype='int64', name='__obs__', length=1000))

  • Attributes: (8)

    created_at :

    2025-07-13T12:29:09.083594+00:00

    arviz_version :

    0.21.0

    inference_library :

    pymc

    inference_library_version :

    5.21.1

    sampling_time :

    12.336488962173462

    tuning_steps :

    500

    modeling_interface :

    bambi

    modeling_interface_version :

    0.15.0

  
  
• log_likelihood

  <xarray.Dataset> Size: 16MB
  Dimensions:      (chain: 4, draw: 500, __obs__: 1000)
  Coordinates:
    * chain        (chain) int64 32B 0 1 2 3
    * draw         (draw) int64 4kB 0 1 2 3 4 5 6 ... 493 494 495 496 497 498 499
    * __obs__      (__obs__) int64 8kB 0 1 2 3 4 5 6 ... 994 995 996 997 998 999
  Data variables:
      rt,response  (chain, draw, __obs__) float64 16MB -0.8921 -1.164 ... -1.179
  Attributes:
      modeling_interface:          bambi
      modeling_interface_version:  0.15.0

  xarray.Dataset

  • Dimensions:
    • chain: 4
    • draw: 500
    • __obs__: 1000
  • Coordinates: (3)
    • chain
      (chain)
      int64
      0 1 2 3

      array([0, 1, 2, 3])

    • draw
      (draw)
      int64
      0 1 2 3 4 5 ... 495 496 497 498 499

      array([  0,   1,   2, ..., 497, 498, 499])

    • __obs__
      (__obs__)
      int64
      0 1 2 3 4 5 ... 995 996 997 998 999

      array([  0,   1,   2, ..., 997, 998, 999])

  • Data variables: (1)
    • rt,response
      (chain, draw, __obs__)
      float64
      -0.8921 -1.164 ... -2.164 -1.179

      array([[[-0.89211002, -1.16416096, -1.74157262, ..., -1.29116593,
               -2.19774012, -1.23217246],
              [-0.85939804, -1.10405186, -1.7483828 , ..., -1.2982562 ,
               -2.18674077, -1.24056155],
              [-0.83830892, -1.12891851, -1.71366241, ..., -1.26750297,
               -2.15073771, -1.21118905],
              ...,
              [-0.88827243, -1.21153371, -1.68716591, ..., -1.2422521 ,
               -2.09630508, -1.18816096],
              [-0.90950666, -1.29416587, -1.69070797, ..., -1.25426613,
               -2.24539089, -1.19126137],
              [-0.90950666, -1.29416587, -1.69070797, ..., -1.25426613,
               -2.24539089, -1.19126137]],
      
             [[-0.80220639, -1.06013377, -1.71239699, ..., -1.28337643,
               -2.1702024 , -1.22596392],
              [-0.73239504, -0.89093473, -1.74686954, ..., -1.33861634,
               -2.26417379, -1.27827023],
              [-1.01148242, -1.37209535, -1.65049286, ..., -1.21337361,
               -2.10477268, -1.15806577],
      ...
              [-0.85592306, -1.08739794, -1.68707219, ..., -1.27703432,
               -2.15854744, -1.22085344],
              [-0.92577891, -1.24395907, -1.64858187, ..., -1.22610385,
               -2.08160689, -1.17246637],
              [-0.92436085, -1.14660565, -1.6854833 , ..., -1.26721903,
               -2.07229792, -1.21598761]],
      
             [[-0.79648859, -1.13895111, -1.67101372, ..., -1.27627423,
               -2.30303456, -1.2107856 ],
              [-0.97835803, -1.28437205, -1.66772511, ..., -1.23601603,
               -2.09786856, -1.1818053 ],
              [-0.83382596, -1.12131307, -1.70876056, ..., -1.28053028,
               -2.23665188, -1.21970116],
              ...,
              [-0.93128737, -1.22248361, -1.70623697, ..., -1.26856145,
               -2.18840754, -1.21000046],
              [-0.86543486, -1.23657433, -1.6442152 , ..., -1.23265566,
               -2.192903  , -1.17212604],
              [-0.9384307 , -1.28507594, -1.67427846, ..., -1.23724561,
               -2.16361154, -1.17916162]]])

  • Indexes: (3)
    • chain
      PandasIndex

      PandasIndex(Index([0, 1, 2, 3], dtype='int64', name='chain'))

    • draw
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             490, 491, 492, 493, 494, 495, 496, 497, 498, 499],
            dtype='int64', name='draw', length=500))

    • __obs__
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             990, 991, 992, 993, 994, 995, 996, 997, 998, 999],
            dtype='int64', name='__obs__', length=1000))

  • Attributes: (2)

    modeling_interface :

    bambi

    modeling_interface_version :

    0.15.0

  
  
• sample_stats

  <xarray.Dataset> Size: 248kB
  Dimensions:                (chain: 4, draw: 500)
  Coordinates:
    * chain                  (chain) int64 32B 0 1 2 3
    * draw                   (draw) int64 4kB 0 1 2 3 4 5 ... 495 496 497 498 499
  Data variables: (12/17)
      smallest_eigval        (chain, draw) float64 16kB nan nan nan ... nan nan
      acceptance_rate        (chain, draw) float64 16kB 0.999 0.9843 ... 0.9941
      perf_counter_diff      (chain, draw) float64 16kB 0.007624 ... 0.008489
      step_size              (chain, draw) float64 16kB 0.5975 0.5975 ... 0.6306
      energy                 (chain, draw) float64 16kB 1.631e+03 ... 1.626e+03
      step_size_bar          (chain, draw) float64 16kB 0.5078 0.5078 ... 0.4455
      ...                     ...
      tree_depth             (chain, draw) int64 16kB 3 3 4 4 3 3 ... 4 3 3 4 3 3
      reached_max_treedepth  (chain, draw) bool 2kB False False ... False False
      perf_counter_start     (chain, draw) float64 16kB 1.427e+06 ... 1.427e+06
      largest_eigval         (chain, draw) float64 16kB nan nan nan ... nan nan
      diverging              (chain, draw) bool 2kB False False ... False False
      max_energy_error       (chain, draw) float64 16kB -0.97 -0.2514 ... -0.0425
  Attributes:
      created_at:                  2025-07-13T12:29:09.093660+00:00
      arviz_version:               0.21.0
      inference_library:           pymc
      inference_library_version:   5.21.1
      sampling_time:               12.336488962173462
      tuning_steps:                500
      modeling_interface:          bambi
      modeling_interface_version:  0.15.0

  xarray.Dataset

  • Dimensions:
    • chain: 4
    • draw: 500
  • Coordinates: (2)
    • chain
      (chain)
      int64
      0 1 2 3

      array([0, 1, 2, 3])

    • draw
      (draw)
      int64
      0 1 2 3 4 5 ... 495 496 497 498 499

      array([  0,   1,   2, ..., 497, 498, 499])

  • Data variables: (17)
    • smallest_eigval
      (chain, draw)
      float64
      nan nan nan nan ... nan nan nan nan

      array([[nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan]])

    • acceptance_rate
      (chain, draw)
      float64
      0.999 0.9843 ... 0.8754 0.9941

      array([[0.99899574, 0.98427376, 0.94007693, ..., 0.88913011, 0.38067985,
              0.17382386],
             [0.96963385, 0.35093074, 0.87165549, ..., 0.95839737, 0.56263805,
              0.98164271],
             [0.82402644, 0.35254835, 0.912601  , ..., 0.96508482, 0.91995735,
              0.86907564],
             [0.96060646, 1.        , 0.89623832, ..., 1.        , 0.87539908,
              0.99405597]])

    • perf_counter_diff
      (chain, draw)
      float64
      0.007624 0.007353 ... 0.008489

      array([[0.007624  , 0.00735321, 0.01615529, ..., 0.00753525, 0.00735346,
              0.00777696],
             [0.00738925, 0.00377908, 0.0084235 , ..., 0.00896825, 0.00190383,
              0.01592554],
             [0.00743321, 0.00786996, 0.00793646, ..., 0.00747246, 0.00717438,
              0.00358404],
             [0.00732617, 0.00734325, 0.00763408, ..., 0.01799921, 0.00955742,
              0.00848917]])

    • step_size
      (chain, draw)
      float64
      0.5975 0.5975 ... 0.6306 0.6306

      array([[0.59752358, 0.59752358, 0.59752358, ..., 0.59752358, 0.59752358,
              0.59752358],
             [0.38784754, 0.38784754, 0.38784754, ..., 0.38784754, 0.38784754,
              0.38784754],
             [0.49690078, 0.49690078, 0.49690078, ..., 0.49690078, 0.49690078,
              0.49690078],
             [0.63061791, 0.63061791, 0.63061791, ..., 0.63061791, 0.63061791,
              0.63061791]])

    • energy
      (chain, draw)
      float64
      1.631e+03 1.63e+03 ... 1.626e+03

      array([[1630.69631792, 1629.62734758, 1627.05643916, ..., 1625.80227866,
              1629.66579447, 1629.7756048 ],
             [1626.32647922, 1632.48370731, 1634.83004223, ..., 1629.90944826,
              1628.46312271, 1628.05114814],
             [1628.45969076, 1631.92495563, 1629.71887201, ..., 1627.93737279,
              1628.60177843, 1632.05979479],
             [1630.82339503, 1631.66538563, 1628.52456439, ..., 1627.95922118,
              1627.68811079, 1626.27266518]])

    • step_size_bar
      (chain, draw)
      float64
      0.5078 0.5078 ... 0.4455 0.4455

      array([[0.50781584, 0.50781584, 0.50781584, ..., 0.50781584, 0.50781584,
              0.50781584],
             [0.49008555, 0.49008555, 0.49008555, ..., 0.49008555, 0.49008555,
              0.49008555],
             [0.50400304, 0.50400304, 0.50400304, ..., 0.50400304, 0.50400304,
              0.50400304],
             [0.4454689 , 0.4454689 , 0.4454689 , ..., 0.4454689 , 0.4454689 ,
              0.4454689 ]])

    • index_in_trajectory
      (chain, draw)
      int64
      3 -1 5 3 -1 -4 ... 5 -1 5 -10 -5 4

      array([[  3,  -1,   5, ...,   4,  -4,   0],
             [  4,   1,  -6, ...,   6,   0,   7],
             [ -5,   0,   2, ...,  -4,   2,  -3],
             [  4,  -4,   7, ..., -10,  -5,   4]])

    • n_steps
      (chain, draw)
      float64
      7.0 7.0 15.0 15.0 ... 15.0 7.0 7.0

      array([[ 7.,  7., 15., ...,  7.,  7.,  7.],
             [ 7.,  3.,  7., ...,  7.,  1., 15.],
             [ 7.,  7.,  7., ...,  7.,  7.,  3.],
             [ 7.,  7.,  7., ..., 15.,  7.,  7.]])

    • process_time_diff
      (chain, draw)
      float64
      0.007551 0.007321 ... 0.008173

      array([[0.007551, 0.007321, 0.015356, ..., 0.00746 , 0.00733 , 0.007642],
             [0.007381, 0.003748, 0.007703, ..., 0.007894, 0.001904, 0.01537 ],
             [0.007412, 0.007561, 0.00774 , ..., 0.007458, 0.007158, 0.003584],
             [0.007253, 0.007335, 0.007509, ..., 0.016259, 0.008841, 0.008173]])

    • energy_error
      (chain, draw)
      float64
      -0.9617 -0.2514 ... 0.043 0.01939

      array([[-0.96167665, -0.25139885,  0.01219149, ...,  0.21640905,
              -0.20854996,  0.        ],
             [ 0.02708104,  2.59162301, -1.47941526, ..., -0.01433663,
               0.        , -0.33364469],
             [-0.20843983,  0.        , -0.01997814, ...,  0.01511287,
              -0.50103452,  0.16767406],
             [ 0.12858918, -0.36713549, -0.113043  , ..., -0.66845625,
               0.04299687,  0.01938668]])

    • lp
      (chain, draw)
      float64
      -1.627e+03 ... -1.625e+03

      array([[-1627.47078996, -1625.58003051, -1624.4655962 , ...,
              -1624.72595789, -1627.00291233, -1627.00291233],
             [-1625.22588946, -1632.01756713, -1629.07734271, ...,
              -1626.44270871, -1626.44270871, -1625.37626727],
             [-1625.12895309, -1625.12895309, -1624.4901249 , ...,
              -1627.14196366, -1626.631181  , -1631.73134296],
             [-1628.98405346, -1627.41478378, -1624.9096195 , ...,
              -1624.39036309, -1625.23857872, -1625.28751373]])

    • tree_depth
      (chain, draw)
      int64
      3 3 4 4 3 3 3 3 ... 3 3 4 3 3 4 3 3

      array([[3, 3, 4, ..., 3, 3, 3],
             [3, 2, 3, ..., 3, 1, 4],
             [3, 3, 3, ..., 3, 3, 2],
             [3, 3, 3, ..., 4, 3, 3]])

    • reached_max_treedepth
      (chain, draw)
      bool
      False False False ... False False

      array([[False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False]])

    • perf_counter_start
      (chain, draw)
      float64
      1.427e+06 1.427e+06 ... 1.427e+06

      array([[1427401.97826979, 1427401.98597425, 1427401.99338738, ...,
              1427406.12062838, 1427406.12823167, 1427406.13566033],
             [1427402.13910204, 1427402.14656562, 1427402.15040271, ...,
              1427406.61256283, 1427406.62165175, 1427406.624541  ],
             [1427402.14002608, 1427402.14753121, 1427402.15554638, ...,
              1427406.53236171, 1427406.53989321, 1427406.54712004],
             [1427402.21389933, 1427402.22129496, 1427402.22872138, ...,
              1427406.63153483, 1427406.64973854, 1427406.65941208]])

    • largest_eigval
      (chain, draw)
      float64
      nan nan nan nan ... nan nan nan nan

      array([[nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan],
             [nan, nan, nan, ..., nan, nan, nan]])

    • diverging
      (chain, draw)
      bool
      False False False ... False False

      array([[False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False],
             [False, False, False, ..., False, False, False]])

    • max_energy_error
      (chain, draw)
      float64
      -0.97 -0.2514 ... 0.2286 -0.0425

      array([[-0.96996264, -0.25139885,  0.19062238, ...,  0.27480448,
               2.12027768,  2.6128521 ],
             [ 0.08181599,  2.7727651 ,  2.22317105, ..., -0.81901335,
               0.57511875, -0.66149511],
             [ 0.47257644,  2.50204722,  0.33240975, ..., -0.91055443,
              -1.09468036,  0.21075921],
             [-0.30839104, -0.59334298,  0.33479679, ..., -0.66845625,
               0.22861677, -0.04249524]])

  • Indexes: (2)
    • chain
      PandasIndex

      PandasIndex(Index([0, 1, 2, 3], dtype='int64', name='chain'))

    • draw
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             490, 491, 492, 493, 494, 495, 496, 497, 498, 499],
            dtype='int64', name='draw', length=500))

  • Attributes: (8)

    created_at :

    2025-07-13T12:29:09.093660+00:00

    arviz_version :

    0.21.0

    inference_library :

    pymc

    inference_library_version :

    5.21.1

    sampling_time :

    12.336488962173462

    tuning_steps :

    500

    modeling_interface :

    bambi

    modeling_interface_version :

    0.15.0

  
  
• observed_data

  <xarray.Dataset> Size: 24kB
  Dimensions:                  (__obs__: 1000, rt,response_extra_dim_0: 2)
  Coordinates:
    * __obs__                  (__obs__) int64 8kB 0 1 2 3 4 ... 996 997 998 999
    * rt,response_extra_dim_0  (rt,response_extra_dim_0) int64 16B 0 1
  Data variables:
      rt,response              (__obs__, rt,response_extra_dim_0) float64 16kB ...
  Attributes:
      created_at:                  2025-07-13T12:29:09.096399+00:00
      arviz_version:               0.21.0
      inference_library:           pymc
      inference_library_version:   5.21.1
      modeling_interface:          bambi
      modeling_interface_version:  0.15.0

  xarray.Dataset

  • Dimensions:
    • __obs__: 1000
    • rt,response_extra_dim_0: 2
  • Coordinates: (2)
    • __obs__
      (__obs__)
      int64
      0 1 2 3 4 5 ... 995 996 997 998 999

      array([  0,   1,   2, ..., 997, 998, 999])

    • rt,response_extra_dim_0
      (rt,response_extra_dim_0)
      int64
      0 1

      array([0, 1])

  • Data variables: (1)
    • rt,response
      (__obs__, rt,response_extra_dim_0)
      float64
      0.8199 1.0 0.5436 ... 1.0 1.668 1.0

      array([[0.81994694, 1.        ],
             [0.54364127, 1.        ],
             [2.20367026, 1.        ],
             ...,
             [1.72080851, 1.        ],
             [2.53323507, 1.        ],
             [1.66819727, 1.        ]])

  • Indexes: (2)
    • __obs__
      PandasIndex

      PandasIndex(Index([  0,   1,   2,   3,   4,   5,   6,   7,   8,   9,
             ...
             990, 991, 992, 993, 994, 995, 996, 997, 998, 999],
            dtype='int64', name='__obs__', length=1000))

    • rt,response_extra_dim_0
      PandasIndex

      PandasIndex(Index([0, 1], dtype='int64', name='rt,response_extra_dim_0'))

  • Attributes: (6)

    created_at :

    2025-07-13T12:29:09.096399+00:00

    arviz_version :

    0.21.0

    inference_library :

    pymc

    inference_library_version :

    5.21.1

    modeling_interface :

    bambi

    modeling_interface_version :

    0.15.0

  
  

Note how the idata include the v parameter which has the __obs__ coordinate. Computing these trial-wise parameters, can serve us well for plotting and other post-hoc analyses, here we would otherwise struggle to be effective without incurring some pre-computation pain (which in addition would be error prone). See e.g. the tutorial on hierarchical Variational Inference for an example.