def std(self, axis):
        if axis not in (0,1):
            raise ValueError('Axis must be 0 or 1.')

        other_axis = axis
        length = self.shape[other_axis]
        sub_mean = self - self.mean(other_axis)
        return (sub_mean.multiply(sub_mean).sum(other_axis)/length).sqrt()
Returns the uncorrected sample standard deviation along the given axis.
dataset/ETHPy150Open jimmycallin/pydsm/pydsm/indexmatrix.py/IndexMatrix.std