powered by
bias corrected double centering # included for speed comparison
doubleCenterBiasCorrected(x)