MI Jul 2013 weighted_mean [ weighted_mean, weighted_mean error] = untitled( x, errors ) gives weighted mean and error on weighted mean where weights = 1/errors^2
0001 function [wmean error_wm] = weighted_mean( x, errors ) 0002 % MI Jul 2013 0003 %weighted_mean [ weighted_mean, weighted_mean error] = untitled( x, errors ) 0004 % gives weighted mean and error on weighted mean where weights = 1/errors^2 0005 0006 weights = 1./(errors).^2; 0007 0008 wmean = sum(weights.*x)/sum(weights); 0009 0010 error_wm = 1./sqrt(sum(weights)); 0011 0012 end 0013