changed gemmFit so that it only calculates pearson's R on the last rep of each datagen. This involved adding an argument to gemmFit, (along with if checks for each time "r" was referenced), and a logical variable to gemmFit that is set to FALSE at the beginning of each datagen loop and TRUE at the beginning of the reps loop whenever reps == n.reps