function [error] = stdEstimator(vec)
    error = sqrt( sum((mean(vec) - vec).^2) / (size(vec, 1)*(size(vec, 1)-1)) );
end