def weighted_std(points, weights): mean = weighted_mean(points, weights) std = sp.sqrt(((points - mean)**2 * weights).sum()) return std