def adc_var ( ip_addr, ip_port, timeout, channel, packetNumber, buffer_size, verbose): s_adc = socket.socket( socket.AF_INET, socket.SOCK_STREAM ) s_adc.settimeout( int(timeout) ) # try connection try: s_adc.connect( (ip_addr, int( ip_port ) ) ) except: print "adc_var(): Error, no ADC connection, port: ", ip_port, s_adc.close() sys.exit() # open connection message = sendReceive ( s_adc, channel, packetNumber, buffer_size, verbose ) # close connection s_adc.close() # unpacking message msg_len = len(message) if msg_len > 0: num = unpack(str(msg_len)+'B', message) # add offset, [0 255] -> [-128 128] num = map(lambda x: x-128, num) if verbose: print "Min=", min(num), "Max=", max(num), "Mean=", stats.mean(num), "Var=", stats.var(num), "Std=", stats.stdev(num) var = stats.var(num) else: var = "N/A" return var
def main(): (opts, args) = parser.parse_args() if len(args) != 1: print "\nPlease, select channel: 0/1\n" parser.print_help() else: if opts.verbose: verbose = True print "connecting to ", opts.ip_addr print "port number", opts.ip_port message = args[0] message = sendReceive(opts.ip_addr, opts.ip_port, opts.timeout, message, opts.verbose) #print "String:" #print message #print "Hex:" hexMessage = b2a_hex(message) #print hexMessage print "No. bytes: ", len(message) print "Done!" #Unpacking message msg_len = len(message) num = unpack(str(msg_len) + 'B', message) #print num num = map(lambda x: x - 128, num) #print num print "Maximum: ", max(num), "Minimum: ", min( num), "Mean: ", stats.mean( num), "Standard deviation: ", stats.stdev( num), "Variance: ", stats.var(num)
print(stats.relfreq(af)) print('\nVARIATION') print('obrientransform:') l = list(range(1, 21)) a = N.array(l) ll = [l] * 5 aa = N.array(ll) print(stats.obrientransform(l, l, l, l, l)) print(stats.obrientransform(a, a, a, a, a)) print('samplevar:', stats.samplevar(l), stats.samplevar(a)) print('samplestdev:', stats.samplestdev(l), stats.samplestdev(a)) print('var:', stats.var(l), stats.var(a)) print('stdev:', stats.stdev(l), stats.stdev(a)) print('sterr:', stats.sterr(l), stats.sterr(a)) print('sem:', stats.sem(l), stats.sem(a)) print('z:', stats.z(l, 4), stats.z(a, 4)) print('zs:') print(stats.zs(l)) print(stats.zs(a)) print('\nTRIMMING') print('trimboth:') print(stats.trimboth(l, .2)) print(stats.trimboth(lf, .2)) print(stats.trimboth(a, .2)) print(stats.trimboth(af, .2)) print('trim1:')
print stats.relfreq(af) print '\nVARIATION' print 'obrientransform:' l = range(1,21) a = N.array(l) ll = [l]*5 aa = N.array(ll) print stats.obrientransform(l,l,l,l,l) print stats.obrientransform(a,a,a,a,a) print 'samplevar:',stats.samplevar(l),stats.samplevar(a) print 'samplestdev:',stats.samplestdev(l),stats.samplestdev(a) print 'var:',stats.var(l),stats.var(a) print 'stdev:',stats.stdev(l),stats.stdev(a) print 'sterr:',stats.sterr(l),stats.sterr(a) print 'sem:',stats.sem(l),stats.sem(a) print 'z:',stats.z(l,4),stats.z(a,4) print 'zs:' print stats.zs(l) print stats.zs(a) print '\nTRIMMING' print 'trimboth:' print stats.trimboth(l,.2) print stats.trimboth(lf,.2) print stats.trimboth(a,.2) print stats.trimboth(af,.2) print 'trim1:'
def test_var(self): "Testing var" data = [self.L, self.A] for d in data: self.EQ(stats.var(d), 35.0)
def evaluate(self, *args, **params): return _stats.var(*args, **params)
def test_var(self): "Testing var" data = [ self.L, self.A ] for d in data: self.EQ( stats.var( d ), 35.0 )
print stats.relfreq(af) print '\nVARIATION' print 'obrientransform:' l = range(1, 21) a = N.array(l) ll = [l] * 5 aa = N.array(ll) print stats.obrientransform(l, l, l, l, l) print stats.obrientransform(a, a, a, a, a) print 'samplevar:', stats.samplevar(l), stats.samplevar(a) print 'samplestdev:', stats.samplestdev(l), stats.samplestdev(a) print 'var:', stats.var(l), stats.var(a) print 'stdev:', stats.stdev(l), stats.stdev(a) print 'sterr:', stats.sterr(l), stats.sterr(a) print 'sem:', stats.sem(l), stats.sem(a) print 'z:', stats.z(l, 4), stats.z(a, 4) print 'zs:' print stats.zs(l) print stats.zs(a) print '\nTRIMMING' print 'trimboth:' print stats.trimboth(l, .2) print stats.trimboth(lf, .2) print stats.trimboth(a, .2) print stats.trimboth(af, .2) print 'trim1:'