jmackereth / asteroestimate

python-based detection probabilities and scaling relations for asteroseismology (Work In Progress!)
MIT License
3 stars 4 forks source link

Write a wrapper function to Monte Carlo over input uncertainties. #5

Open grd349 opened 4 years ago

grd349 commented 4 years ago

So this will be the main way that people can interact with the software. It will be a function (inside a class?) with inputs and out the prior as the summary statistics of a Gaussian distribution and the samples from the Monte Carlo.

grd349 commented 4 years ago

So the function will look like ...

def NameOfFunc(J, H, K, parallax, ...) ''' J, H, K, parallax all tuples where first value is estimate and second value is uncertainty '''

Monte Carlo over inputs 

return (numax_mean, numax_std), samples

Where samples is a list of the values from the MC.

jmackereth commented 4 years ago

Ok, this is implemented in the branch mcmc_inputs here: #7