Enable rhat and ess for tidy.brmsfit() #149
Merged
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
|
Hmm, this is a nice idea and I really appreciate PRs, but ... could you see if you can easily resolve the checks that are failing ... ? |
alyst
force-pushed
the
brmsfit_rhat_ess
branch
from
67a2241
to
4b38575
Compare
added 2 commits
February 8, 2024 14:29
somehow applying rhat/ess_bulk on the draws_array directly resulted in wrong values, it needs to be passed to summarize_draws()
|
Tests fixed |
|
Since it turned out that using posterior::summarise_draws() was required to correctly calculate both rhat and ess, I can change the code to also use summarise_draws() for mean/median, std/mad and quantile calculation -- so that all the metrics (except HDPinterval) are calculated by one call. |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
Add support for calculating Rhat and ESS convergence metrics to tidy.brmsfit2.
For Rhat, it uses brms::rhat() method, for ESS, it uses posterior::ess_basic().
Actually, the posterior package has support for vectorized calculation of many metrics, including mean, std, quantiles etc (using MatrixStats under the hood), so it could be potentially used in the other places instead of
apply(samples, 2, fun)to improve the performance.