# cancer1 
# Fig. 33.5 

import random
from reader import *
from corrCoef import dot

NUM_TRIALS = 32000

rdr = readVals("cancers.txt")
xs = readFloat1D(rdr)
print("X:", xs)
ys = readFloat1D(rdr)
print("Y:", ys)
cor = dot(xs,ys)
print(f"dot corr = {cor:.4f}")
count = 0
perms = ys.copy()
for j in range(NUM_TRIALS): 
  random.shuffle(perms)
  if dot(xs, perms) >= cor:
    count += 1
print(f"count = {count}; p = {(count/NUM_TRIALS):.5f}")