The results of Polygenic Risk Scores(PRS)
1
0
Entering edit mode
5.0 years ago
ysx0226 • 0

PRSice was used to calculate the PRS in ADNI dataset, but the results are quite abnormal ,see belows:

When I use the ADNI data in PRSice dir,the results R2 were quite large (eg:0.749622);

When I use the example data in PRSice dir,the results R2 were quite small (eg:0.0133696);

PRSice.prsice on ADNI

Set Threshold   R2  P   Coefficient Standard.Error  Num_SNP

Base    0.0001  0.749622    8.43404e-23 56.3465 5.7326  33

Base    0.00015 0.829959    4.56913e-21 80.9499 8.59464 54

PRSice.prsice on example

Set Threshold   R2  P   Coefficient Standard.Error  Num_SNP

Base    0.00025 0.0133696   8.43169e-06 -0.197266   0.0442903   2

Base    0.0003  0.00824473  0.000456434 -0.225204   0.0642503   3
SNP prsice • 2.0k views
ADD COMMENT
0
Entering edit mode

I added prsice tag, as I believe the developer of the program watches this tag.

ADD REPLY
0
Entering edit mode

what do you mean by "add prsice tag", is this a parameter of the program?

ADD REPLY
1
Entering edit mode

Nope, he mean adding a tag to your post so that I can be notified.

ADD REPLY
0
Entering edit mode

Thank you for your advice~

ADD REPLY
1
Entering edit mode
4.9 years ago
Sam ★ 4.7k

You are working on two different dataset, it is therefore normal that they are going to have a different result (different R2).

However, here, we have a special case where you've obtain an abnormally high R2 for your ADNI dataset. When it comes to PRS analysis, we usually expect a small R2, with a high R2, it is usually because you have use the same dataset to generate the GWAS summary statistic AND use the same set of sample to generate the PRS. You should make sure that an independent dataset is used to generate the GWAS summary statistics.

ADD COMMENT
0
Entering edit mode

Thank you, according to your comment, my results seem to be over-fitting

I think I am using a dataset to predict itself, that is a mistake

ADD REPLY

Login before adding your answer.

Traffic: 2457 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6