Gnuplot средневзвешенная ошибка

Question

У меня есть набор данных с y баров ошибок. Я могу найти среднее значение данных без ошибок, используя следующее:

f(x)=mean_y
fit f(x) "data" via mean_y

Тем не менее, я хочу найти взвешенное среднее с учетом ошибок. Возможно ли что-то подобное в gnuplot или мне нужно это кодировать?

Заранее спасибо.

Answer 1

Вы можете использовать using параметр для указания столбца для ошибок. С тремя using спецификаторы, третий интерпретируется как стандартное отклонение s и используется для вычисления веса 1/s**2 для соответствующего значения:

f(x) = mean_y
fit f(x) "data" using 1:2:3 via mean_y

Это предполагает, что ваш data файл имеет три столбца, x, z а также stdev, Если у вас есть только два столбца, z а также stdev, вы должны использовать

fit f(x) "data" using 0:1:2 via mean_y