solution

Cho hai mẫu A (n_A phần tử) và B (n_B phần tử). Thống kê T = |mean(A) - mean(B)|. Hãy thực hiện B_perm = 5000 hoán vị: với rng = np.random.default_rng(seed), mỗi lần dùng rng.permutation(n_A + n_B) để hoán vị mảng kết hợp [A,B], sau đó tách thành n_A phần tử đầu (gán A') và n_B phần tử sau (gán B'), tính |mean(A') - mean(B')|. P-value = (số lần |mean(A')-mean(B')| >= T) / B_perm. In T và p_value cách nhau dấu cách, làm tròn 4 chữ số thập phân.