from itertools import combinations
from scipy.stats import ttest_ind
for list1, list2 in combinations(list(data.keys()), 2):
    t, p = ttest_ind(data[list1], data[list2])
print(list1, list2, p)
How to calculate p-value for two lists of floats?