Calculates the first two moments of the bivariate distribution of NB_model and NB_all
calc_NB_moments(Y, pi, z, weights = NULL)
Two means, two SDs, and one correlation coefficient. First element is for the model and second is for treating all
Vector of the binary response variable
Vector of predicted risks
Decision threshold at which the NBs are calculated
Optinal - observation weights