Our group would like to perform logistic regression of rare variants across 10 or so moderate sized data sets. In order to perform a score test combining summary statistics over these data sets, we would need the components of the score test for each data set: the score function and information matrix.

This is implemented in, for example, the RVTESTS and RAREMETAL software.

I’m not a statistician myself, but it looks like for single variant analysis, that software generated an information matrix which was a 1x1 matrix.

Can you add that feature to Hail? Or is there another way to output those components?

Thanks!