diff options
Diffstat (limited to 'data/combined/mean.py')
| -rwxr-xr-x | data/combined/mean.py | 7 |
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) |