-
Notifications
You must be signed in to change notification settings - Fork 34
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
Priors and Hyperpriors #28
Comments
@arwelHughes Thanks for the feedback. Good idea. It may take some time to add to the documentation as people are currently focused on getting the next release out, which has been long overdue. |
May I intrude in this ticket to ask a somewhat related question: what kind of distribution domains can Paramonte handle? It doesn't seem clear from the help & manual pages. For example, can it deal with a distribution defined on the Cartesian product of the real line, the integers, and a discrete (nominal/categorical) set? Cheers and thank you for the great project! |
@pglpm Thank you for your inquiry, and I apologize for the delayed response, as some of us have been traveling. Cartesian is the only one supported in the current release. We have had internal discussions of adding support for discrete domains that have so far not gained traction as we did not need it, and no one has asked for it. If this interests you, I'd appreciate it if you create an issue dedicated to this request to materialize the need and justify the effort to add it. We have been working on a new major release that has been delayed for too long, but hopefully not much longer. |
@fagheri Thank you for the info! Then I'll create an issue ticket. Looking forward to the new release. I've been working on "nonparametric density inference" for inference problems in neuroscience and medicine. In medicine it often happens that the variates are a combination of nominal, ordinal, interval types (psychometric tests, family attributes, clinical tests, and similar). Ordinal variates can be dealt with with a continuous distribution, by using latent variables. But nominal variates really require a sampler that works on discrete spaces. I'm trying to build a package for nonparametric density inference with this general kind of data, especially for use in medicine. The goal is to make it as user-friendly and transparent as possible (the pre-package is available here, an example application is here). I'm working in R and the Monte Carlo sampling is taken care of by the package Nimble. But eventually I'd like to write something directly in Fortran. This is why I'm interested in Paramonte. |
I would also be very interested in this. |
Thanks again for your feedback @Peku995 and @pglpm ! |
A concrete use case for me would be a scattered light simulation where a sphere size distribution scatters the light. |
I can imagine a similar use-case, where we have sets of Molecular Dynamics simulations of biomembranes at say, different pressures or salt concentrations. From these we build the model to fit our data, but they are also a discreet set (which could be in a database), but the final model is constructed of continuous functions. But, I guess the way of approaching this classically would be to use Bayesian model selection (i.e. nested sampling - hint! ;) ) to choose between the possibilities using a set of fits.... |
Could we have a general discussion / examples of how one might implement priors and hyperpriors in the objective function somewhere in the documentation?
The text was updated successfully, but these errors were encountered: