This function gives the gradient of the log likelihood of the Box-Cox model. Main purpose is to be an input to the maxLik function.
ScoreFun(bb, ixx, iyy, iw, iZZ)the gradient of the log likelihood estimate for the coefficients in `bb`
initial values for the intercept and slope coefficients
continuous predictor
binary outcome
sample weight
covariates to be incorporated in the model