Thank you @jessegrabowski , I used the .sel(chain=0) and below is the output:
This is an xarray with 291 rows and each row containing 500 samples of the likelihood , I think.
How do I get 291 rows with a single value , something like np.average on each row might work.
But not sure how to do that on an xarray.
Can you please let me know?