Comments (4)
Yes, I can observe the same behaviour, though the order can be the opposite as well, and vary from sample to sample.
The number of evidence points for BOLFI affects the shape of the fitted GP, which in turn affects the shape of the posterior. The NUTS sampler is a self-adjusting, recursive and random algorithm, and possibly chooses different internal settings in these cases. Due to the recursion the differences in runtime can be significant. Also, the default initialization points are taken from the evidence with minimal discrepancy, and these can have a big effect as well.
from elfi.
oh great, makes sense! I forgot about the NUTS peculiarities....
Thank you!
from elfi.
Is NUTS the only available sampler within bolfi.sample? How do I find the list of parameters/options for bolfi.sample in the documentation?
from elfi.
There's an open PR #296 to add the existing Metropolis sampler to the BOLFI.sample
interface, but for the time being, NUTS is your only built-in option. Basically you could of course give the BOLFI posterior logpdf to any external sampler (such as elfi.methods.mcmc.metropolis).
The standard Python docstring shows the available arguments for BOLFI.sample
, available also here:
https://elfi.readthedocs.io/en/latest/api.html#elfi.BOLFI.sample
from elfi.
Related Issues (20)
- threshold value returned when specifying 'n_sim' in Rejection sampler HOT 1
- Use Dask for parallel computing? HOT 5
- Can rejection sampling be done with no explicit passing of observed data to a simulator? HOT 2
- Regression adjustment returns point estimate in specific case HOT 3
- networkx 2 compatibility HOT 3
- stochastic_optimization call to differential_evolution broken in newer scipy HOT 1
- Implement initial standard BSL HOT 1
- Implement initial semiBSL method HOT 1
- Tightly pinned requirements are becoming out of date HOT 2
- elfi.examples docstrings are missing details HOT 1
- BOLFIRE target nodes
- Extend sample result summaries HOT 1
- Issue in trying to setup RandMaxVar in BOLFI HOT 4
- Require `plot_gp` inputs `const` and `bounds` to be given as dictionaries
- Random variation in default kernel parameters
- BOLFI posterior.py logging
- Processing failed simulations / samples HOT 2
- BOLFI fails if numpy version is >= 1.20 HOT 1
- Basic rejection method of inferring parameters of an SIR disease model HOT 4
- import elfi fails HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from elfi.