print MeanOutliers.nummeandevs(None, 90, 2, 3) print MeanOutliers.nummeandevs(None, 20, 2, 3) print MeanOutliers.nummeandevs(None, 21, 2, 3) print MeanOutliers.nummeandevs(None, 23, 2, 3) print MeanOutliers.nummeandevs(None, 25, 2, 3) print MeanOutliers.nummeandevs(None, 15, 2, 3) print MeanOutliers.nummeandevs(None, 20, 2, 3) print MeanOutliers.nummeandevs(None, 25, 2, 3) print MeanOutliers.nummeandevs(None, 25, 2, 3) print MeanOutliers.nummeandevs(None, 25, 2, 3) print MeanOutliers.nummeandevs(None, 60, 2, 3) """ meansum = 0.0 meandevsum = 0.0 count = 0 std = 40 incount = 0 for i in range(1,10000): val = random.gauss(0,std) if abs(val) < 2*std: incount += 1 meansum += 1.0*val count += 1 meandevsum += abs(1.0*(val - meansum/count)) MeanOutliers.nummeandevs(None, val, 2, 3) print meansum/count, meandevsum/count print "incount", incount print MeanOutliers.nummeandevs(None, 80, 2, 3) pass