powered by
PQL function for binary data
binomial_PQL(Y, X_sig1 = NULL, Beta, Z, Alpha)
PQL estimate for binary data