new script

1 files changed, 5 insertions, 2 deletions

diff --git a/data/combined/mean.py b/data/combined/mean.py
index 2499beb..dd1f56a 100755
--- a/data/combined/mean.py
+++ b/data/combined/mean.py
@@ -23,15 +23,18 @@ for line in sk_file:
     if "-1" not in line and 2 < float(line[labels.index('z-value')]) < 3:
         id = line[labels.index('name')]
         for x in labels[beg:]: data[id][x] += [float(line[labels.index(x)])]
-
+var = []
 for i in data:
     out = [i]
     for l in labels[beg:]:
         if len(data[i][l]) > 1:
+            var += [np.var(data[i][l])]
             data[i][l] = np.average(data[i][l])
             out += [str(np.average(data[i][l]))]
     #if len(out) > 1:
     #    print ','.join(out)
+sys.stderr.write(str(div**0.5)+'\n')
+sys.stderr.write(str(np.average(var))+'\n')
 
 sk_file.seek(0)
 print sk_file.readline().strip()
@@ -42,6 +45,6 @@ for line in sk_file:
         obs = float(line[labels.index(l)])
         if obs > 0:
             diff = obs-data[id][l]
-            line[labels.index(l)] = str(obs - diff/div)
+            line[labels.index(l)] = str(data[id][l] + diff/(div**0.5))
     print ','.join(line)