bd-j / prospector

Python code for Stellar Population Inference from Spectra and SEDs
http://prospect.readthedocs.io
MIT License
153 stars 71 forks source link

Dirichlet SFH prior #332

Closed MohammadRih closed 1 month ago

MohammadRih commented 2 months ago

Hello,

I'm trying to set up a model to use dirichlet prior for sfh but have trouble doing that. The problem is that I cannot find concentration parameter, $\alpha_{D}$, in the TemplateLibrary["dirichlet_sfh"]. How can I set up this parameter in my work?!

bd-j commented 2 months ago

Hi @MohammadRih the concentration parameter vector alpha can be adjusting by changing the parameters alpha and beta of the prior distributions for the elements of z_fraction at https://github.com/bd-j/prospector/blob/main/prospect/models/templates.py#L676

For symmetric dirichlet all alpha and beta should be changed to the same value.

See https://arxiv.org/abs/1010.3436 and https://arxiv.org/abs/1811.03637 for more info.

MohammadRih commented 2 months ago

Thanks for your kind and useful comment.