Answers.Body: <p>The coefficient is the third column from the right, you can tell because some are positive and other negative. When the value is negative the association is a negative correlation. When it is positive (last two) that is more like the result you are seeking. The further away from zero the stronger the correlation (positive or negative)</p>
  <p>These are very large values for the coefficients. The SNP is the last column on the right. I would disgard the top two results because the probabilities are too borderline, e.g. if you used Bonferoni's correction.</p>
  <p>There isn't enough information on this calculation nor the data set to make any further insights other than you have 7 SNPs showing PRS</p>
  
Answers.ContentLicense: CC BY-SA 4.0
Answers.CreationDate: 2020-10-21T15:49:36.963Z
Answers.LastActivityDate: 2020-10-21T15:49:36.963Z
Answers.UserDisplayName: M__
Body: <p>I am using PRSice to compute the PRS over a train set and want to use the coefficient used on the train set to apply it on another set which I will call the test set.</p>
  <p>Once I compute the PRS I get a file called analysis.prsice</p>
  <p>This is the few first lines</p>
  <hr />
  <pre><code>Pheno Set Threshold R2 P Coefficient Standard.Error Num_SNP
  
  - Base 5e-08 0.00550157 0.0659631 -21.0513 11.4492 43
  
  - Base 5.005e-05 0.00396637 0.113759 -61.8545 39.1108 184
  
  - Base 0.00010005 0.00236122 0.220463 -68.9527 56.2742 286
  
  - Base 0.00015005 0.00113035 0.394503 -59.4058 69.7677 382
  
  - Base 0.00020005 0.00072253 0.495129 -55.2605 81.0064 471
  
  - Base 0.00025005 0.000334179 0.642121 -41.5674 89.4432 538
  
  - Base 0.00030005 7.5306e-05 0.825111 -21.391 96.8021 601
  
  - Base 0.00035005 3.95773e-06 0.959527 5.25804 103.612 668
  
  - Base 0.00040005 0.000175506 0.735018 37.723 111.455 743
  </code></pre>
  <hr />
  <p>I thought that I could apply the PRS on another dataset by multiplying the coefficient in the result by the allelic dosage of the corresponding SNPs. However, I am not sure how to know which are the corresponding SNPs of the coefficients.</p>
  <p>Also, if there was a command that does that and it would be easier, could you please tell me how to do it?</p>
  <p>Thanks</p>
  
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2020-10-17T21:06:37.847Z
Comments.Text: How does PRSice work, i.e. what ML is it using? What is the accurarcy of your predictions? Applying a trained data set to predict another data set should be straight forward. You don't need to specify individual coefficients
Comments.UserDisplayName: M__
Comments.ContentLicense: CC BY-SA 4.0
Comments.CreationDate: 2020-10-19T12:37:50.587Z
Comments.Text: I am not sure what you mean? actually, I am not even sure how you measure PRS accuracy. I know in the C+T method that you add the allelic dosage weighted by the effect size in the GWAS
Comments.UserDisplayName: lalaland
Tags: <snp><gwas><phylogenetics>
Title: Apply trained PRS on another dataset
ThreadText: # Apply trained PRS on another dataset
  
  (Asked by: username_0 on 2020-10-17)
  
  I am using PRSice to compute the PRS over a train set and want to use the coefficient used on the train set to apply it on another set which I will call the test set.
  
  Once I compute the PRS I get a file called analysis.prsice
  
  This is the few first lines
  
  ---
  
  ```
  Pheno Set Threshold R2 P Coefficient Standard.Error Num_SNP
  
  - Base 5e-08 0.00550157 0.0659631 -21.0513 11.4492 43
  
  - Base 5.005e-05 0.00396637 0.113759 -61.8545 39.1108 184
  
  - Base 0.00010005 0.00236122 0.220463 -68.9527 56.2742 286
  
  - Base 0.00015005 0.00113035 0.394503 -59.4058 69.7677 382
  
  - Base 0.00020005 0.00072253 0.495129 -55.2605 81.0064 471
  
  - Base 0.00025005 0.000334179 0.642121 -41.5674 89.4432 538
  
  - Base 0.00030005 7.5306e-05 0.825111 -21.391 96.8021 601
  
  - Base 0.00035005 3.95773e-06 0.959527 5.25804 103.612 668
  
  - Base 0.00040005 0.000175506 0.735018 37.723 111.455 743
  
  ```
  
  ---
  
  I thought that I could apply the PRS on another dataset by multiplying the coefficient in the result by the allelic dosage of the corresponding SNPs. However, I am not sure how to know which are the corresponding SNPs of the coefficients.
  
  Also, if there was a command that does that and it would be easier, could you please tell me how to do it?
  
  Thanks
  
  **Tags:** snp, gwas, phylogenetics
  
  > **Comments:**
  > **username_1**: How does PRSice work, i.e. what ML is it using? What is the accurarcy of your predictions? Applying a trained data set to predict another data set should be straight forward. You don't need to specify individual coefficients  
  > **username_0**: I am not sure what you mean? actually, I am not even sure how you measure PRS accuracy. I know in the C+T method that you add the allelic dosage weighted by the effect size in the GWAS  
  
  ---
  
  ## 1 answers
  
  ### Answer by username_1 on 2020-10-21. Score: 0 | Accepted answer
  
  The coefficient is the third column from the right, you can tell because some are positive and other negative. When the value is negative the association is a negative correlation. When it is positive (last two) that is more like the result you are seeking. The further away from zero the stronger the correlation (positive or negative)
  
  These are very large values for the coefficients. The SNP is the last column on the right. I would disgard the top two results because the probabilities are too borderline, e.g. if you used Bonferoni's correction.
  
  There isn't enough information on this calculation nor the data set to make any further insights other than you have 7 SNPs showing PRS
  
  
  ---
  
  