After reading other discussions, I have found out that the regression functionality is suitable for any general regression regardless of the actual outcome/exposure (not only genetic variants).
My current analysis plan is to run a linear regression with the same set of exposure and covariates but with (many many) different outcomes. I’m actually trying to run a regression with all common variants as outcomes.
I tried to supply ‘mt.GT.n_alt_alleles()’ in the y variable but it didn’t work. After looking at the documentation, it seems that the regression function currently gets only 1d arrays as an input (Hail | Statistics).
Is there a solution? I think people who wants to run a PheWAS also require having multiple outcomes option when running a regression.