tinder_experiments_monte_carlo.py

import random

N_EXPERIMENTS = 10000
N_SWIPES = 250
A_REAL_RATE = 66 / 250
B_REAL_RATE = 70 / 250


def experiment():
    """ observable rate for A and B """
    a = sum(random.random() < A_REAL_RATE for i in range(N_SWIPES))
    b = sum(random.random() < B_REAL_RATE for i in range(N_SWIPES))
    return a, b


experiments_resuls = [experiment() for i in range(N_EXPERIMENTS)]
error_rate = sum(a > b for a, b in experiments_resuls) / N_EXPERIMENTS
print(f'error_rate = {error_rate}')