This function calculates the edge-wise partial correlation difference for a single bootstrap.
single.boot(i, z, n, tY.org, P, levels.z, w.upper)
iteration number. This is not used within this function, but necessary for parSapply within scoring.boot.parallel function.
a length n vector representing a binary covariate
the number of rows in data
the transformed standardized data
the global correlation component
the levels of the covariates
the upper triangular of Omega
the difference for this bootstrap