What is the mean of $S_i$?

In [15]:

    

import urllib2
import csv
import numpy

data = open('../data/data.csv', 'r').readlines()
fieldnames = ['x', 'y', 'z', 'unmasked', 'synapses']
reader = csv.reader(data)
reader.next()

rows = [[int(col) for col in row] for row in reader]


    

In [16]:

    

rows[0]


    

    
Out[16]:





[19, 1369, 55, 5063, 0]

In [20]:

    

mean = numpy.mean([r[-1] for r in rows])
no_um_mean = numpy.mean([r[-1] for r in rows if r[-2] == 0])
only_um_mean = numpy.mean([r[-1] for r in rows if r[-2] != 0])


    

In [22]:

    

print(mean, no_um_mean, only_um_mean)


    

    




(124.71150608650609, 0.0, 139.61649843243146)

In [23]:

    

var = numpy.var([r[-1] for r in rows])
no_um_var = numpy.var([r[-1] for r in rows if r[-2] == 0])
only_um_var = numpy.var([r[-1] for r in rows if r[-2] != 0])


    

In [24]:

    

print(var, no_um_var, only_um_var)


    

    




(8463.553588148694, 0.0, 7394.0989170154053)