Question on computing PGS from beta_auto and G

[scienception]

Hi Florian!
I had a question regarding this matrix, vector multiplication in file code/example-with-provided-LD.R. Why is beta_auto multiplied by map_pgs2$beta? I thought it would be beta_auto only.

pred_auto <- big_prodVec(G, beta_auto * map_pgs2$beta,
                         ind.col = map_pgs2[["_NUM_ID_"]],
                         ncores = NCORES)

[scienception]

map_pgs2$beta seems to be a vector of 1s. I guess my question is why bother multiplying by 1?

[privefl]

It's in case there is some allele switching needed to match between the two datasets.