bug fixed in generating Y, using for loop now. Moved the process to compute CI into pset2_util.R