search

s-goldman / Dusty-Evolved-Star-Kit

SED-fitting python package for fitting evolved stars
https://dusty-evolved-star-kit.readthedocs.io
Other
15 stars 2 forks source link

Bayesian fitting #88

Open s-goldman opened 4 years ago

s-goldman commented 4 years ago

Use a prior, maybe the luminosity functions from the Magellanic Clouds, to constrain the luminosities. Then use observations to generate a posterior distribution. This would be a problem for highly-evolved stars which may need a flag in the user interface.

s-goldman commented 4 years ago

This would solve the problems of fitting both spectra and photometry, the unrealistic weighting of data in the far-IR and optical, unrealistic DUSTY models, and finally yield some version of an errorbar

s-goldman commented 4 years ago

I would need an option in the user interface to choose the luminosity function for either LMC, SMC, (or M31)

s-goldman commented 4 years ago

Functions have been created to calculate the grid weights, and next is priors. The way that the code now works is that it creates a grid of luminosities between 1,000 and 100,000 Lsun with a specified sampling in logarithmic spacing. I need to take the specified grid and duplicate and append it for each of the luminosities.

s-goldman commented 3 years ago

Also add option to optimize with emcee.