-
Notifications
You must be signed in to change notification settings - Fork 0
/
ttest.py
69 lines (51 loc) · 2.23 KB
/
ttest.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
"""Approximately simulates trec_eval using pytrec_eval."""
import argparse
from audioop import avg
import os
import sys
from scipy import stats
import pytrec_eval
def get_result(run, qrel, args):
relevance_threshold = 2 if args.tag == 'cast20' else 1
evaluator = pytrec_eval.RelevanceEvaluator(
qrel, {'ndcg_cut.3', 'recip_rank'}, relevance_level=relevance_threshold)
results = evaluator.evaluate(run)
ndcg_cut_3 = []
recip_rank = []
for qid in results:
recip_rank.append(results[qid]['recip_rank'])
ndcg_cut_3.append(results[qid]['ndcg_cut_3'])
return recip_rank, ndcg_cut_3
def main():
parser = argparse.ArgumentParser()
parser.add_argument('--qrel')
parser.add_argument('--run1')
parser.add_argument('--run2')
parser.add_argument('--tag')
args = parser.parse_args()
assert os.path.exists(args.run1)
assert os.path.exists(args.run2)
with open(args.qrel, 'r') as f_qrel:
qrel = pytrec_eval.parse_qrel(f_qrel)
with open(args.run1, 'r') as f_run1:
run1 = pytrec_eval.parse_run(f_run1)
with open(args.run2, 'r') as f_run2:
run2 = pytrec_eval.parse_run(f_run2)
run1_recip_rank, run1_ndcg_cut_3 = get_result(run1, qrel, args)
run2_recip_rank, run2_ndcg_cut_3 = get_result(run2, qrel, args)
# for i in range(len(run1_recip_rank)):
# print(run1_recip_rank[i] - run2_recip_rank[i])
# print(run1_recip_rank - run2_recip_rank)
# run1_recip_rank = [float('%.3f' % item) for item in run1_recip_rank]
# run1_ndcg_cut_3 = [float('%.3f' % item) for item in run1_ndcg_cut_3]
# run2_recip_rank = [float('%.3f' % item) for item in run2_recip_rank]
# run2_ndcg_cut_3 = [float('%.3f' % item) for item in run2_ndcg_cut_3]
print(sum(run2_recip_rank)/len(run2_recip_rank))
print(sum(run1_recip_rank)/len(run1_recip_rank))
print(sum(run2_ndcg_cut_3)/len(run1_ndcg_cut_3))
print(sum(run1_ndcg_cut_3)/len(run1_ndcg_cut_3))
# print(stats.levene(run1_recip_rank, run2_recip_rank))
print(stats.ttest_rel(run1_ndcg_cut_3,run2_ndcg_cut_3,alternative='greater'))
print(stats.ttest_rel(run1_recip_rank,run2_recip_rank,alternative='greater'))
if __name__ == "__main__":
sys.exit(main())