## Load the script for the principal component (PC) and the correlation-based statistics

source("ldtstk.r")

## Now call the function. 
## Parameters below: data_file; 3 numbers of PC; number of permutations

DeltaTest("testdat-recoded.txt", 1, 2, 3, 1000)