geomean.m

% Copyright (C) 2010-2017, Raytheon BBN Technologies and contributors listed
% in the AUTHORS file in TASBE analytics package distribution's top directory.
%
% This file is part of the TASBE analytics package, and is distributed
% under the terms of the GNU General Public License, with a linking
% exception, as described in the file LICENSE in the TASBE analytics
% package distribution's top directory.

function m = geomean(data,weights)
    if (isempty(data))
        m = NaN;
    else
        if nargin < 2
            m = 10.^mean(log10(data));
        else
            valid = weights>0;
            m = 10.^wmean(log10(data(valid)),weights(valid));
        end
    end
end