wrote function for heteroskedasticity test for general e,X as well as for a BART model... started examining heteroskedasticity in the BHD, wrote functions of significance for gammas, augmented kfoldcv to handle Z matrices