Я хочу извлечь положение пика из спектра (энергетический спектр рассеянных фотонов). Для этого я использую scipy.optimize.curve_fit, чтобы подогнать гауссиану к области спектра, которая напоминает гауссову.
Как найти неопределенность пикового значения? Само пиковое значение будет дано результатом для среднего параметра из регрессии Гаусса.
Есть две вещи, которые пришли мне на ум:
- Я получаю значения ковариации из процедуры минимизации, из которой я получаю ошибку среднего параметра.
- Кроме того, я мог бы подумать об использовании сигмы Гаусса, чтобы получить ошибку среднего.
Мои мысли по этому поводу заключались бы в том, что ошибка в среднем параметре не может быть неправильным путем. И я бы также поспорил, что стандартная ошибка на самом деле не говорит нам о неопределенности, с которой мы знаем пиковое значение. Он говорит нам о форме распределения, но не о неопределенности пикового значения (которое для простоты мы считаем истинным, четко определенным значением).
(Это репост вопроса, который я первоначально разместил на stats.stackoverflow, где я не получил ответов через 2 дня.)
