Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

WAIC #230

Open
wants to merge 3 commits into
base: develop
Choose a base branch
from
Open

WAIC #230

wants to merge 3 commits into from

Conversation

KSkwarczynski
Copy link
Member

@KSkwarczynski KSkwarczynski commented Nov 29, 2024

Pull request description

Implement WAIC based on Gellman paper
image
https://arxiv.org/pdf/1307.5928

Changes or fixes

-Extend usage of vectors instead of raw pointers

Examples

Using some dummy values I obtained WAIC very similiar to DIC. It make sense both methods give similair results
image
[Edit]
WAIC after Fix
image

Copy link

Hi @KSkwarczynski, thank you for contributing to MaCh3! Please wait for MaCh3 developers to review your PR. If no one answers within a week, please message people from this list: https://github.com/orgs/mach3-software/teams/mach3admin . While waiting, please enjoy this

meme

Use this action on your projects. Use jokes on issues instead.

@github-actions github-actions bot added Documentation Improvements or additions to documentation MCMC MCMC related Samples labels Nov 29, 2024
// ****************
MACH3LOG_INFO("******************************");
switch(Criterion) {
// TSpline3 (third order spline in ROOT)
Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Does this comment really belong here?

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

no crap...

@@ -2615,7 +2541,61 @@ void SampleSummary::StudyDIC() {
const double DIC_stat = Dhat + 2 * p_D;
MACH3LOG_INFO("Effective number of parameters following DIC formalism is equal to: {:.2f}", p_D);
MACH3LOG_INFO("DIC test statistic = {:.2f}", DIC_stat);
MACH3LOG_INFO("******************************");
}

Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Since you refer to "equation 3" later on in the code, please put the reference in the header before the function declaration, or here in the implementation before the function declaration

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

They are in header file this way it should be readable in doxygen

/// @cite Gelman2014

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Documentation Improvements or additions to documentation MCMC MCMC related Samples
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants