function ss = stderr(data)
    ss = std(data)/sqrt(length(data));
end