Uhh… good question! It seems to be working as expected with the order from `DynamicPPL.syms(DynamicPPL.VarInfo(m))`

. It’s also the order of the data in the chain’s underlying AxisArrays (see `ch.value`

).

I wonder what the `ch.name_map.parameters`

order is about… Just displaying the chain (`display(ch)`

) shows this order in the header, and the other order in the summary stats:

```
Chains MCMC chain (1000×14×1 Array{Float64, 3}):
Iterations = 1:1:1000
Number of chains = 1
Samples per chain = 1000
Wall duration = 10.55 seconds
Compute duration = 10.55 seconds
parameters = μ, σ
internals = lp, n_steps, is_accept, acceptance_rate, log_density, hamiltonian_energy, hamiltonian_energy_error, max_hamiltonian_energy_error, tree_depth, numerical_error, step_size, nom_step_size
Summary Statistics
parameters mean std naive_se mcse ess rhat e ⋯
Symbol Float64 Float64 Float64 Float64 Float64 Float64 ⋯
σ 1.1765 0.5388 0.0170 0.0246 411.3738 1.0006 ⋯
μ 1.9466 0.6001 0.0190 0.0182 295.0682 0.9990 ⋯
```

I had filed a related issue here but the most relevant is probably this one.