The idea of the regression code lifted up to Python sounds like it would be beneficial to others as well for writing their own custom regression function. My understanding is that we only need one element of the score vector, and one element of the inverse of the Fisher’s information matrix, similar to how the Wald test returns one number for beta, instead of the full beta vector. I’ve looked at LogisticRegressionModel.scala, and I think we’ll try modifying that to return the two extra values. If we run into questions, would the Hail dev forum be the right place to post them?
Thanks!