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

Add section about Pathfinder diagnostic and using for inits #833

Merged
merged 3 commits into from
Nov 18, 2024

Conversation

avehtari
Copy link
Contributor

Summary

Related to #751 and recent discussion about using Pathfinder for inits in https://discourse.mc-stan.org/t/using-pathfinder-to-initalize-sampling/37221

I now just added brief sections on Pareto-k diagnostic and using Pathfinder for inits. It would be good to add much description about the Pathfinder overall.

Copyright and Licensing

Please list the copyright holder for the work you are submitting (this will be you or your assignee, such as a university or company):

By submitting this pull request, the copyright holder is agreeing to license the submitted work under the following licenses:

Copy link
Contributor

@bob-carpenter bob-carpenter left a comment

Choose a reason for hiding this comment

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

Thanks for writing this up. Just some small clarifying comments and typos and references.

## Diagnosing Pathfinder

Pathfinder diagnoses the accuracy of the approximation by computing the density ratio of the true posterior and
the approximation and using Pareto-$\hat{k}$ diagnostic (Vehtari et al., 2024) to assess whether these ratios can
Copy link
Contributor

Choose a reason for hiding this comment

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

We should use actual references here like "@zhang_pathfinder:2022"


Pathfinder diagnoses the accuracy of the approximation by computing the density ratio of the true posterior and
the approximation and using Pareto-$\hat{k}$ diagnostic (Vehtari et al., 2024) to assess whether these ratios can
be used to improve the approximation via resmapling. /, the
Copy link
Contributor

Choose a reason for hiding this comment

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

resmapling is misspelled.

I don't know what the / is doing here---should this be a different character?

be used to improve the approximation via resmapling. /, the
normalization for the posterior can be estimated reliably (Section 3, Vehtari et al., 2024), which is the
first requirement for reliable resampling. If estimated Pareto-$\hat{k}$ for the ratios is smaller than 0.7,
there is still need to further diagnose importance sampling estimates by taking into account also the expetant
Copy link
Contributor

Choose a reason for hiding this comment

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

I don't know the word "expetant"---looks like it's Latin. Can we just write this in English?

function (Section 2.2, Vehtari et al., 2024). If estimated Pareto-$\hat{k}$ is larger than 0.7, then the
estimate for the normalization is unreliable and any Mote Carlo estimate may have a big error. The resampled draws
can still contain some useful information about the location and shape of the posterior which can be used in early
parts of Bayesian workflow (Gelman et al, 2020).
Copy link
Contributor

Choose a reason for hiding this comment

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

All of these should get turned into citations.

@avehtari
Copy link
Contributor Author

I fixed references and several typos

Copy link
Contributor

@bob-carpenter bob-carpenter left a comment

Choose a reason for hiding this comment

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

Looks good---thanks.

@bob-carpenter bob-carpenter merged commit 82c0051 into master Nov 18, 2024
@bob-carpenter
Copy link
Contributor

P.S. I fixed the merge conflict in the .bib file. Hope I didn't break anything in the process.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

2 participants