function coeffOfVariation = calcCoeffOfVariation(x, varargin)

if nargin > 1
    unbiased = varargin{1};
else
    unbiased = 0;
end

if unbiased
    coeffOfVariation = nanstd(x) / nanmean(x) * (1 + 1 / 4 / numel(x(~isnan(x))));
else
    coeffOfVariation = nanstd(x) / nanmean(x);
end

end