Entering edit mode

3 months ago

Alejandro Rojo
▴
20

Hi all,

I have an annotated vcf file for one individual which I want to estimate his polygenic risk score (PRS) for a certain trait, using the scoring files from the PGSCatalog. The scoring file contains the SNP ID, reference and alternative allele, and weights. How can I estimate the PRS using the scoring file without using the classical approach of having GWAS summary data, target data, etc...?

Thank you very much for your consideration.

Thank you very much for your answer. One last question. Since the allele represent the single point mutation, should I use dummy coding to transform the nucleotides and then perform the multiplication with the betas or is there another approach?

Yes, if effect allele is

`"A"`

and genotype is`"A A"`

, then`2 * beta`

Does anyone have any software or script that performs these calculations?

It is a one-liner, multiply genotypes with coefficients and sum them

Maybe try PRSice-2: https://www.prsice.info/