# cancer2.py 

import random
from reader import *
from corrCoef import *

NUM_TRIALS = 32000

rdr = readVals("cancers.txt")
xs = readFloat1D(rdr)
print("X:", xs)
ys = readFloat1D(rdr)
print("Y:", ys)
cor = corrCoef(xs,ys)
print(f"pearson corr = {cor:.4f}")

count = 0
perms = ys.copy()
for j in range(NUM_TRIALS): 
  random.shuffle(perms)
  if corrCoef(xs, perms) >= cor:
    count += 1
print(f"count = {count}; p = {(count/NUM_TRIALS):.5f}")