def paired_difference_test(a_counts, a_nums, b_counts, b_nums):
""" Run a paired difference t-test on the event counts
represented by the two input parameters.
Args:
Returns:
Raises:
"""
assert len(set(a_counts.keys()) ^ set(b_counts.keys())) == 0
assert len(set(a_nums.keys()) ^ set(b_nums.keys())) == 0
a_means = []
b_means = []
for src_id in sorted(a_counts.keys()):
a_cnt = a_counts[src_id]
a_num = a_nums[src_id]
a_mean = float(a_cnt) / float(a_num)
b_cnt = b_counts[src_id]
b_num = b_nums[src_id]
b_mean = float(b_cnt) / float(b_num)
print '%.2f\t%.2f\t%d\t%d' % (a_mean, b_mean, a_num, b_num)
a_means.append(a_mean)
b_means.append(b_mean)
paired_diff_t_test(a_means, b_means)
def paired_difference_test(a_counts, a_nums, b_counts, b_nums):
""" Run a paired difference t-test on the event counts
represented by the two input parameters.
Args:
Returns:
Raises:
"""
assert len(set(a_counts.keys()) ^ set(b_counts.keys())) == 0
assert len(set(a_nums.keys()) ^ set(b_nums.keys())) == 0
a_means = []
b_means = []
for src_id in sorted(a_counts.keys()):
a_cnt = a_counts[src_id]
a_num = a_nums[src_id]
a_mean = float(a_cnt) / float(a_num)
b_cnt = b_counts[src_id]
b_num = b_nums[src_id]
b_mean = float(b_cnt) / float(b_num)
print '%.2f\t%.2f\t%d\t%d' % (a_mean,b_mean,a_num,b_num)
a_means.append(a_mean)
b_means.append(b_mean)
paired_diff_t_test(a_means, b_means)
def diff_test(counts_a, nums_a, counts_b, nums_b):
""" Paired difference t-test.
"""
# Sanity check
keys_a = counts_a.keys()
keys_b = counts_b.keys()
assert len(set(keys_a) ^ set(keys_b)) == 0
means_a = []
means_b = []
print 'A\tB\t#A\t#B'
for i in sorted(keys_a):
means_a.append(float(counts_a[i]) / float(nums_a[i]))
means_b.append(float(counts_b[i]) / float(nums_b[i]))
print '%.2f\t%.2f\t%d\t%d' % (means_a[i],means_b[i],nums_a[i],nums_b[i])
paired_diff_t_test(means_a,means_b)
def diff_test(counts_a, nums_a, counts_b, nums_b):
""" Paired difference t-test.
"""
# Sanity check
keys_a = counts_a.keys()
keys_b = counts_b.keys()
assert len(set(keys_a) ^ set(keys_b)) == 0
means_a = []
means_b = []
print 'A\tB\t#A\t#B'
for i in sorted(keys_a):
means_a.append(float(counts_a[i]) / float(nums_a[i]))
means_b.append(float(counts_b[i]) / float(nums_b[i]))
print '%.2f\t%.2f\t%d\t%d' % (means_a[i], means_b[i], nums_a[i],
nums_b[i])
paired_diff_t_test(means_a, means_b)
