Hi - hail newbie here. I am trying to run linear regression on 8 phenotypes. They can be arranged into 4 sets of 2 each, where each set of 2 has the same missingness structure. My code is:
lin_reg = hl.linear_regression_rows(
y=[[filtered_mt.pheno.A1,filtered_mt.pheno.A2],
[filtered_mt.pheno.B1,filtered_mt.pheno.B2],
[filtered_mt.pheno.C1,filtered_mt.pheno.C2],
[filtered_mt.pheno.D1,filtered_mt.pheno.D2]],
x=filtered_mt.GT.n_alt_alleles(),
covariates=covariates
)
where A1 & A2 have the same missingness pattern, as do B1/B2, etc.
The output for beta, t-stat, p-value, etc is an array, however, it seems to duplicate information:[[1.94e-01,1.93e-01],[1.94e-01,1.93e-01],[-9.65e-02,1.41e-01],[-9.65e-02,1.41e-01], …
Sorry I can’t figure out how to show the entire string for the first locus, but it has 8 sets within , rather than just 4 (the first 4 are displayed above). You can see that the results for the phenotypes A1 and A2 are duplicated (two sets with 1.94e-01 and 1.93e-01) as are the results for B1 and B2 (two sets with -9.65e-02 and 1.41e-01).
What am I doing wrong to have the output duplicated? And how do I parse these to get at the individual results?
Apologies if my question is unclear!