AbstractMCMC.jl icon indicating copy to clipboard operation
AbstractMCMC.jl copied to clipboard

Correct syntax for chainscat and chainsstack

Open astro-kevin opened this issue 2 years ago • 3 comments

I'm not sure what's the correct syntax for the functions chainscat and chainsstack. I tried calling chainscat(a,b) and chainsstack([a,b]) and both return an error. (As an aside, is the only difference between the two one takes chains as arguments and one takes a vector of chains?) In my test code, both a and b are MCMCChains objects.

astro-kevin avatar May 23 '23 15:05 astro-kevin

Could you paste the error you're experiencing?

On Tue, May 23, 2023, at 8:22 AM, astro-kevin wrote:

I'm not sure what's the correct syntax for the functions chainscat and chainsstack. I tried calling chainscat(a,b) and chainsstack([a,b]) and both return an error. (As an aside, is the only difference between the two one takes chains as arguments and one takes a vector of chains?) In my test code, both a and b are MCMCChains objects.

— Reply to this email directly, view it on GitHub https://github.com/TuringLang/AbstractMCMC.jl/issues/122, or unsubscribe https://github.com/notifications/unsubscribe-auth/AADHITX4W724DJTDPPHUTL3XHTI5FANCNFSM6AAAAAAYMCAPTA. You are receiving this because you are subscribed to this thread.Message ID: @.***>

cpfiffer avatar May 24 '23 22:05 cpfiffer

Here is the error: ERROR: ArgumentError: chain ranges differ. I think that's because of the burn in which is discarded. However, I'm just not really understanding how this works because the object which it returns when it does work (when I'm stacking any except the first chain), it returns an object with multiple chains, rather than one long chain as I expect. Are the stats presented on the multiple chains equivalent to the one with very many samples?

astro-kevin avatar May 25 '23 16:05 astro-kevin

For this it might be more helpful if you posted the whole script you're using, or, even better, a minimum working example -- is that something you're able to do?

cpfiffer avatar May 25 '23 19:05 cpfiffer