annayqho / TheCannon

a data-driven method for determining stellar parameters and abundances from stellar spectra
MIT License
39 stars 16 forks source link

Label estimate uncertainties #43

Closed tingyuansen closed 8 years ago

tingyuansen commented 9 years ago

The code formally derives uncertainties for the label estimates. If I understand correctly, the calculation does not include the coefficient uncertainties. Only the flux uncertainties are included. I am guessing the uncertainty of these "transferred-labels" might be too optimistic.

davidwhogg commented 9 years ago

To use these (to first order), we need to compute, at the best-fit label position, the full derivative d(label) / d(coefficient) and then multiply them around the covariance matrix. Or something like that...

tingyuansen commented 9 years ago

I agree. But I think it should be something like d(label)/d(coefficient,flux), since flux also contains uncertainties in this case. We can probably also assume cov(coefficient,flux) = 0. The question is however whether we can estimate cov(coefficient) in the first place -- the topic in another post. We should continue our fight there :) .

davidwhogg commented 9 years ago

Yes, we already have cov(coefficient).

tingyuansen commented 9 years ago

I would argue that we might not have that (see the other post on coefficient uncertainties).