Неопределенность в повторяющихся измерениях

Предположим, вы получили некоторые данные измерения и у вас есть следующие числа:

  • измерение 1: 1,510 ± 0,085
  • измерение 2: 1,608 ± 0,089
  • измерение 3: 1,566 ± 0,059
  • измерение 4: 1,638 ± 0,066
  • измерение 5: 1,660 ± 0,071

Итак, у меня есть 5 значений с соответствующими неопределенностями. Теперь я хотел бы взять среднее значение этих 5 и сообщить о стандартном отклонении. Перед глазами три сценария:

  1. Я беру среднее значение и вычисляю sd, используя только значения без ошибок, и получаю среднее значение 1,596 ± 0,060 ;

  2. Я вычисляю sd, используя формулу распространения ошибок, где Вопрос "=" А + Б , затем д Вопрос "=" ( д А 2 + д Б 2 ) 1 / 2 , и я суммирую в квадратуре все 5 ошибок, указанных выше, и получаю значение 1,596 ± 0,168 , значительно больше по сравнению со случаем 1;

  3. Я беру среднее значение ошибок 5 измерений, которые у меня есть, со значением 1,596 ± 0,074 .

Я действительно запутался, потому что не могу понять, как справиться со случаем, когда вам нужно усреднить некоторые данные с ошибками. Если у кого-то есть какие-то подсказки и намеки, это было бы очень признательно.

Ответы (4)

В общем, всегда следует использовать все доступные данные, чтобы делать выводы. Если вам так повезло, что у вас есть границы неопределенности для ваших повторных измерений, это означает, что эти неопределенности также являются своего рода данными, и в идеале их следует использовать для уточнения заключения.

Что бы я сделал, так это использовал взвешенный метод: если мы предполагаем, что шум и неопределенность являются гауссовыми, то точкам с низкой неопределенностью следует присвоить больший вес, чем точкам с большей неопределенностью. См. этот вопрос и ответ на него . В этом случае мы хотим максимизировать логарифмическую вероятность бревно ( л ) "=" с ( 1 / 2 ) я ( у я мю ) 2 / о я 2 где а искомое среднее и о я неопределенность для каждой из точек данных. Производная я ( у я мю ) / о я 2 что равно нулю для взвешенного среднего

мю "=" я у я / о я 2 я ( 1 / о я 2 ) .
Точно так же можно рассчитать взвешенное стандартное отклонение :
о "=" Н я ( у я мю ) 2 / о я 2 ( Н 1 ) я ( 1 / о я 2 ) .

Но учтите, что многое зависит от модели, которую вы используете. Неопределенности не обязательно должны быть гауссовыми.

Большое спасибо за разъясняющий ответ! не только вам, Андерс, но и другим пользователям. Это мне очень помогло.
Вышеупомянутое является средним значением и стандартным отклонением для выборок и является наилучшей оценкой неизвестного среднего значения и стандартного отклонения совокупности.

Здесь нужно быть осторожным. Предположим, что данные, которые вы сообщаете в своем вопросе, представляют собой средние значения и стандартные отклонения средних значений для серии случайных выборок. То есть, 1,510 ± 0,085 это среднее ± стандартное отклонение среднего значения для одного образца, состоящего из ряда отдельных измерений, 1,608±0,089 — среднее значение ± стандартное отклонение среднего значения для другого образца, состоящего из ряда отдельных измерений, и т. д.

Рассмотрим ряд я "=" 1 , 2 , . . . , к случайные выборки, т. я т час образец, состоящий из н я конкретные значения для интересующей случайной переменной. Каждая выборка может иметь разное количество конкретных значений. Для каждого образца вы оцениваете и сообщаете среднее значение и стандартное отклонение среднего значения. Среднее значение я т час образец м я "=" 1 н я Дж н я у Дж я где у Дж я это Дж т час значение в я т час образец. Стандартное отклонение среднего для я т час образец С я "=" с я 2 н я где с я "=" Дж н я ( у Дж я м я ) 2 н я 1 является стандартным отклонением для я т час образец.

Наилучшей оценкой среднего является м б е с т "=" я "=" 1 к м я / С я 2 я "=" 1 к 1 С я 2 и наилучшая оценка стандартного отклонения среднего С б е с т "=" ( я "=" 1 к 1 С я 2 ) 1 / 2 . Вы сообщаете м б е с т ± С б е с т для вашего конечного результата.

Примечание. Примерные значения означают м б е с т и стандартное отклонение С б е с т являются наилучшей оценкой неизвестных средних значений генеральной совокупности мю и стандартное отклонение среднего о мю .

(Например, подробности см. в тексте «Анализ данных для ученых и инженеров» Мейера.)

Андерс С уже дал вам очень хороший ответ, но я хочу показать вам, где вы ошиблись в своих рассуждениях.

Я вычисляю sd, используя формулу распространения ошибок, где Вопрос "=" А + Б , затем д Вопрос "=" ( д А 2 + д Б 2 ) 1 / 2 , и я суммирую в квадратуре все 5 ошибок, указанных выше, и получаю значение 1,596 ± 0,168 , значительно больше по сравнению со случаем 1;

Это вычисляет ошибку в сумме ваших измерений. С помощью этого метода вы получаете, что сумма ваших измерений равна 7,982 ± 0,168 .

Но чтобы получить среднее значение, вы взяли эту сумму, а затем разделили на количество измерений . Поскольку количество измерений является точным числом без неопределенности, вы можете получить ошибку в среднем, разделив ошибку в сумме на то же число.

Так что у вас должно быть 1,596 ± 0,036 этим методом, а не ± 0,168 . И вы бы улучшили свою ошибку по сравнению с одним измерением.

Однако, используя метод Андерса, вы получите еще меньшую ошибку.

Если ваши измерения нормально распределены (и систематическими погрешностями можно пренебречь), вы можете использовать среднее значение значений и стандартное отклонение среднего (SDOM) в качестве окончательной неопределенности. Это просто

С Д О М : о Икс ¯ "=" о Икс Н
где о Икс - это просто стандартное отклонение ваших измерений, а N - ваше количество измерений.

В книге Джона Тейлора «Введение в анализ ошибок» это описывается как «неопределенность в нашей наилучшей оценке x (а именно Икс ¯ ) является... SDOM" с учетом вышеуказанных условий.