 interpolation between isochrones neither accounts for the nonlinear mapping of time onto the H-R diagram nor the non-uniform distribution of stellar masses observed in the galaxy. As a consequence, straightforward interpolation between isochrones results in an age distribution for field stars that is biased towards older ages compared to the distribution predicted by stellar evolutionary theory.

Bayesian inference of stellar age and mass aims to eliminate such a bias by accounting for observationally and/or theoretically motivated distribution functions for the physical parameters of interest. As an example, for a given point with error bars on the H-R diagram, a lower stellar mass should be considered more likely due to the initial mass function. Likewise, due to the longer main-sequence timescales for lower mass stars, a star that is observed to have evolved off the main sequence should have a probability distribution in mass that is skewed towards higher masses, i.e. because higher mass stars spend a more significant fraction of their entire lifetime in the post-MS stage.

### Bayes Formalism {#subsubsec:formalism}

Bayesian estimation of the physical parameters can proceed from comparison of the data with a selection of models. Bayes’ Theorem states:

$$P(\mathrm{model|data}) \propto P(\mathrm{data|model}) \times P(\mathrm{model})$$

The probability of a model given a set of data is proportional to the product of the probability of the data given the model and the probability of the model itself. In the language of Bayesian statistics, this is expressed as:

$$\mathrm{posterior} \propto \mathrm{likelihood} \times \mathrm{prior}.$$

Our model is the set of stellar parameters, age ($\tau$) and mass ($M_*$), and our data are the measured effective temperature, $T_\mathrm{eff}$, and surface gravity, $\log g$, for a given star. At any given combination of age and mass, the predicted $T_\mathrm{eff}$ and $\log g$ are provided by stellar evolutionary models. The $\chi^2$ statistic for an individual model can be computed as follows:

$$\begin{aligned}
\chi^2 (\tau, M_*) &= \sum \frac{(O-E)^2}{\sigma^2} \\ 
 &=  \frac{[(T_\mathrm{eff})_O-(T_\mathrm{eff})_E]^2}{\sigma_{T_\mathrm{eff}}^2
interpolation between isochrones neither accounts for the nonlinear mapping of meter onto the H - R diagram nor the non - consistent distribution of stellar masses observed in the galax. As a consequence, straightforward interpolation between isochrones consequence in an age distribution for field star that is biased towards older long time compared to the distribution predicted by leading evolutionary theory. 

 Bayesian inference of stellar age and mass aims to eliminate such a bias by account for observationally and/or theoretically motivated distribution functions for the physical parameters of pastime. As an example, for a given point with error parallel bars on the H - R diagram, a lower stellar mass should be considered more likely due to the initial mass affair. Likewise, due to the longer main - sequence timescales for lower mass stars, a star that is observed to have evolved off the main sequence should have a probability distribution in mass that is skewed towards higher masses, i.e. because high mass stars spend a more meaning fraction of their integral lifetime in the post - MS stage. 

 # # # Bayes Formalism { # subsubsec: formalism } 

 Bayesian estimation of the physical argument can proceed from comparison of the data with a selection of models. Bayes ’ Theorem country: 

 $ $ P(\mathrm{model|data }) \propto P(\mathrm{data|model }) \times P(\mathrm{model})$$ 

 The probability of a model given a set of data is proportional to the product of the probability of the data give the model and the probability of the model itself. In the language of Bayesian statistics, this is express as: 

 $ $ \mathrm{posterior } \propto \mathrm{likelihood } \times \mathrm{prior}.$$ 

 Our model is the stage set of stellar parameters, age ($ \tau$) and mass ($ M_*$), and our data are the measured effective temperature, $ T_\mathrm{eff}$, and open gravity, $ \log g$, for a given star. At any grant combination of age and mass, the predicted $ T_\mathrm{eff}$ and $ \log g$ are put up by stellar evolutionary model. The $ \chi^2 $ statistic for an individual model can be computed as be: 

 $ $ \begin{aligned } 
 \chi^2 (\tau, M _ *) & = \sum \frac{(O - E)^2}{\sigma^2 } \\ 
  & =   \frac{[(T_\mathrm{eff})_O-(T_\mathrm{eff})_E]^2}{\sigma_{T_\mathrm{eff}}^2
 inherpolation between isocmrones neither aexounts for tge nonlivear mapping of time onto thx H-R diageam nor the non-uniform distribunion of sreller masses observxs in thc galzwy. As e consequence, sjraightforwasd interpolatimn bztween isochrones results in an age qistribitlon for field ftarx thaf is biased towards older ages comlared tm the distribition predicted by stellar evoputionary theory.

Bajesian infetsncq of stellar xge and mass aims to emiminate such a bias by accountkng fpr observajninapny and/or thxoretibally motivatcc distsibutiom functions fov the phtsical parameters of mnterest. As an exampje, for a ciren point with error vaes on the H-R auagfam, a lkwer shelmar mass sgould be cobsidered more likelu qlr to the inifial mwsf function. Likewise, due to the longer mdin-aequence timescales for lower mass stars, a sjar that if observed to have evolved off the main sequence vhoulv favt q procqbllity distribution in mass that is skewed towwddx migher masses, i.e. because hibhfr iass stars spgnd a more significant fractlon of jheir wntire lisetike in the post-MS stage.

### Bayew Formalism {#fybsubsec:formalism}

Bcyesian estioatipn of the physical parameterr cah proceed fgom compadkson of the data winh a selection of models. Bayes’ Theorem wtatzs:

$$P(\mathro{modgl|data}) \[ropto P(\mahhrm{data|model}) \times P(\matjrm{mobel})$$

Tha probabillty of a model given a set of deva is proportootal to the 'roducb of the probabylity of the dcta giveu the oodel and nhe probauility of thq model itselx. In the langnage of Bwyesuan wtatistkzs, this is exptessed as:

$$\mathrm{postwrior} \propto \mathrk{liyslihood} \times \mcuhrn{prior}.$$

Our model is thq dev of fdellar paramaterr, aer ($\tau$) and mass ($M_*$), xnd pur data are the meavures effective temperstmre, $T_\mathtm{eff}$, and surface gravoty, $\log g$, for a gigen svar. At any bivgn combination of age and mass, fhe predifteq $T_\mathrm{eff}$ and $\log g$ are pxovided by stellar evolutionary models. Tie $\chi^2$ statistic for an individual model ccn be computev as fjllows:

$$\beghn{aligned}
\chi^2 (\tau, M_*) &= \wum \frac{(O-E)^2}{\sigma^2} \\ 
 &=  \frac{[(T_\mathrm{eff})_O-(T_\mafhrm{efx})_E]^2}{\sigla_{T_\mathrm{eff}}^2
interpolation between isochrones neither accounts for the of onto the diagram nor the observed the galaxy. As consequence, straightforward interpolation isochrones results in an age distribution field stars that is biased towards older ages compared to the distribution predicted stellar evolutionary theory. Bayesian inference of stellar age and mass aims to eliminate a by for and/or theoretically motivated distribution functions for the physical parameters of interest. As an example, for a point with error bars on the H-R diagram, lower stellar mass should considered more likely due to initial function. Likewise, to longer timescales for lower stars, a star that is observed to have evolved off the main sequence should have a probability in mass skewed towards masses, because mass stars spend significant fraction of their entire lifetime stage. ### Bayes Formalism {#subsubsec:formalism} Bayesian estimation of physical parameters proceed from comparison of the data a selection of models. Bayes’ Theorem states: $$P(\mathrm{model|data}) P(\mathrm{data|model}) \times P(\mathrm{model})$$ The probability of a model given a set of data is proportional product of the probability the data given model the of model itself. the language of Bayesian statistics, this is expressed as: $$\mathrm{posterior} \propto \times \mathrm{prior}.$$ Our model is the set of stellar parameters, and ($M_*$), and our are the measured effective $T_\mathrm{eff}$, surface gravity, $\log g$, given At of and the predicted $T_\mathrm{eff}$ and g$ are provided by stellar models. The $\chi^2$ statistic be computed as follows: $$\begin{aligned} \chi^2 (\tau, M_*) \sum \frac{(O-E)^2}{\sigma^2} \\ &= \frac{[(T_\mathrm{eff})_O-(T_\mathrm{eff})_E]^2}{\sigma_{T_\mathrm{eff}}^2
 interpolation between isochRones neithEr accOunTs fOr The nOnliNear mapping of tIMe onTo the H-R diagram nor the noN-unifOrM DistRIbUtion Of stellAR mASSes ObSeRveD iN ThE galaXy. AS a conseQuence, straIghTfOrward interpOLaTion betweeN isOchrones resuLts In an agE dIstRIbutiOn fOr fieLd starS That is Biased towArDS older AGes compAREd To thE distribution predICtED by stellar evolUtionaRy THeORY.

BaYesIan inferenCe Of steLLar age aND mASS AimS To eliminate suCh a bias by acCOunTing foR oBseRVationAlly aNd/OR thEoretically MotiVated distRibutiON functiONs for thE physiCal ParAmetERs Of IntErESt. AS An ExaMPle, For a giveN pOiNt witH errOR BARs on The h-R diAgram, A lower stellar MasS shoULd bE consIdereD morE lIkely Due to tHe iniTiAl mass function. LIkewIse, due to tHe lOnGer MaIn-seqUEnce tiMesCalEs for loWer mass STarS, a STAR tHat is observed to havE eVOLvEd off the Main seQUeNcE Should haVe A prObabILIty diStriBUtIon in masS that iS SkEwEd towarDs Higher MaSseS, i.e. BecauSE higHer masS stars spEnd a mORe significant fRAction of their ENtIRE lIFetiMe iN the post-MS sTage.

### bAyes formALiSm {#sUBsubsEc:forMaLIsM}

bayesian estimation oF tHe physIcal pArameters can pRoceed from COMParison oF the DAtA With a selection Of modEls. Bayes’ ThEOrem statEs:

$$P(\maThrm{modeL|data}) \propTO p(\mathrm{dAta|ModEl}) \tImeS p(\MaThrm{model})$$

The pROBabiLiTy of a moDel Given a sEt oF daTa iS prOpOrtional tO the prodUcT oF tHe ProBabilITy of the dAtA giVeN thE modeL And the ProbaBiliTy Of THe mOdel itsELf. iN The lAnGuAge oF BaYeSian sTatiSTicS, this is Expressed As:

$$\mAThrm{PoStErior} \prOpto \mathrm{likElIhood} \times \MaThrM{prior}.$$

oUR model is The set of stellar parameteRS, age ($\tau$) And Mass ($M_*$), And oUr data are The MeasurEd eFFectivE tempeRaturE, $T_\MatHRM{eff}$, aND SuRfaCe Gravity, $\log G$, FOr a Given StAr. At Any giveN combination of age aND maSs, the predicteD $T_\mAthrM{EFf}$ And $\LOg G$ Are PrOVidED By stellar evolutIonary modeLs. tHe $\Chi^2$ statistIC foR aN indiviDual modEl can BE computEd as folloWs:

$$\begin{alIgNed}
\cHI^2 (\Tau, m_*) &= \sum \frac{(O-E)^2}{\Sigma^2} \\ 
 &=  \fraC{[(T_\mathrm{eFF})_O-(T_\maTHrM{eff})_E]^2}{\SigMa_{T_\matHrM{efF}}^2
 interpolation between iso chrones ne ither  ac cou nt s fo r th e nonlinear ma p ping  of time onto the H-Rdiagr am nort he  non- uniform di s t rib ut io n o fs te llarmas ses obs erved in t hega laxy. As a c o ns equence, s tra ightforwardint erpola ti onb etwee n i sochr ones r e sultsin an age  d i stribu t ion for f ie ld s tars that is bias e dt owards older a ges co mp a re d tothe  distribut io n pre d icted b y  s t e l lar evolutionarytheory.

Ba y esi an inf er enc e  of st ellar  a g e a nd mass aim s to  eliminat e such a biasb y accou ntingfor  ob serv a ti on all ya nd/ o rthe o ret ically m ot iv ateddist r i b u tion  fu ncti ons f or the physic alpara m ete rs of  inte rest .As an  examp le, f or  a given pointwith  error ba rson  th eH-R d i agram,  alow er stel lar mas s  sh ou l d be  considered more l ik e l ydue to t he ini t ia lm ass func ti on.  Lik e w ise,duet othe long er mai n -s eq uence t im escale sfor  lo wer m a ss s tars,a star t hat i s  observed to h a ve evolved of f  t h e  m a in s equ ence should  hav e  a p roba b il ity distr ibuti on in mass that is skewed  t owards  high er masses, i. e. because h i gher mas s st a rs spend a more s ignif icant frac t ion of t heirentire l ifetime i n the post -MS  st age .
# # #Bayes Formali s m  {#s ub subsec: for malism} 

B aye sia n e st imation o f the ph ys ic al  p ara meter s  can pro ce edfr omcompa r ison o f the  dat awi t h a  select i on o f mo de ls . Ba yes ’Theor em s t ate s:

$$P (\mathrm{ mod e l|da ta })  \propt o P(\mathrm{d at a|model})\t ime s P(\m a t hrm{mode l})$$

The probabilityo f a mod elgiven  a s et of dat a i s prop ort i onal t o theprodu ct  of t he pr o b ab ili ty  of the da t a  gi ven t he  mod el andthe probability of the  model itself . I n th e la ngu a ge ofBa y esi a n  statistics, th is is expr es s ed  as:

$$\m a thr m{ posteri or} \pr opto\ mathrm{ likelihoo d} \times  \ math r m {pr ior}.$$

O ur model  is the s e t ofs te llarpar ameter s,  ag e ($\ tau$)a ndmass($M_*$ ),  and o ur da ta  are the  measured effective tem peratu re, $ T_\ mathrm{ef f}$ ,  an d surface  gra vity, $\lo g g $,for a  gi v en st ar.A tany given  com b ination o f  a gea n dmass, the p r e d ict ed $T _\m a thrm{e ff}$  and $\log g$ are provided by st ella r evo lut i onar ymodels. The $\ chi ^2 $ statisti cfor an indi vidual m od e l can  be co mputed  as fol l o ws : 

$$\b egin {al igned}
\c hi^ 2( \tau, M _* )& = \sum  \fr ac {(O-E) ^2}{\s i gma^ 2 }  \\ 
 &=  \frac{ [(T_\ m a thrm{ e ff} )_O-( T_ \mathrm { eff} )_E]^2}{\s igma_{T_\ma thrm{e ff}} ^2
 interpolation_between isochrones_neither accounts for the_nonlinear mapping_of_time onto_the_H-R diagram nor_the non-uniform distribution_of stellar masses observed_in the galaxy._As_a consequence, straightforward interpolation between isochrones results in an age distribution for field stars_that_is biased_towards_older_ages compared to the distribution_predicted by stellar evolutionary theory.

Bayesian_inference of_stellar age and mass aims to eliminate such_a_bias by accounting_for observationally and/or theoretically motivated distribution functions for the_physical parameters of interest. As an_example, for a_given_point_with error bars on_the H-R diagram, a lower stellar_mass should be considered more likely_due to the initial mass function. Likewise,_due to the longer main-sequence timescales_for lower mass stars, a_star that_is observed to have evolved_off the main_sequence should_have a probability_distribution in mass that is skewed_towards higher masses,_i.e. because higher mass stars spend_a_more significant fraction_of_their_entire lifetime_in the post-MS_stage.

###_Bayes Formalism_{#subsubsec:formalism}

Bayesian_estimation of the physical parameters can_proceed_from comparison of the data with a_selection of models. Bayes’_Theorem_states:

$$P(\mathrm{model|data}) \propto P(\mathrm{data|model}) \times_P(\mathrm{model})$$

The probability of a model_given a set of data is_proportional to_the product_of the probability of the data given the model and the_probability of the model itself. In_the language of Bayesian_statistics, this_is_expressed as:

$$\mathrm{posterior} \propto_\mathrm{likelihood}_\times \mathrm{prior}.$$

Our_model is the set of stellar parameters,_age ($\tau$)_and mass ($M_*$), and our data_are the measured effective_temperature,_$T_\mathrm{eff}$, and surface gravity, $\log g$,_for a given star. At any_given combination of age and_mass,_the_predicted $T_\mathrm{eff}$ and $\log g$_are provided by stellar evolutionary models._The $\chi^2$ statistic_for an individual model can be computed_as_follows:

$$\begin{aligned}
\chi^2 (\tau, M_*) &= \sum \frac{(O-E)^2}{\sigma^2}_\\_
 &=  \frac{[(T_\mathrm{eff})_O-(T_\mathrm{eff})_E]^2}{\sigma_{T_\mathrm{eff}}^2