def std(y):
    """ Returns unbiased standard deviation. """
    from math import sqrt
    import numpy as N
    l = len(y)
    s = N.std(y)
    if l == 1:
        return s
    else:
        return s * sqrt(float(l) / (l - 1))
