McCullagh_compute_regression_weights function

Computes regression weights w; R_dot_j * (N - R_dot_j[j]) * (n_do_j[j] a= na_dot_j[j+ 1] )