<p>I'm a little confused about what you mean by there being no phenotype. To build the score you need a phenotype, although to apply it you might not. (Usually PRS is build with a base dataset and then applied using a target dataset). Here's a good tutorial about doing PRS in plink that might help you:</p>
<p><a href="https://choishingwan.github.io/PRS-Tutorial/plink/" rel="nofollow noreferrer">https://choishingwan.github.io/PRS-Tutorial/plink/</a></p>
<p>Also this paper is a really good paper that helped me better understand how building PRS works:</p>
<p><a href="https://www.nature.com/articles/s41596-020-0353-1" rel="nofollow noreferrer">https://www.nature.com/articles/s41596-020-0353-1</a></p>
<p>Hope this helps!</p>

CC BY-SA 4.0
2024-02-29T17:30:56.033Z
2024-02-29T17:30:56.033Z
Sarah
<p>I've been struggling with constructing a PRS score. In the pics you can see how my summary stats and my target data look like.</p>
<p>Target Data Summary Stats</p>
<p>I've tried the following code which should be pretty straightforward:</p>
<pre><code>plink \
  --bfile target \
  --score summary_stats.txt 1 4 6 header \
  --out final
</code></pre>
<p>When I run it, I get this message:</p>
<pre><code>10499994 variants and 1198 people pass filters and QC. Note: No phenotypes present.
</code></pre>
<p>The data is binary and the fam file really doesn't have phenotypes (it's -9). Is there any way to still run a PRS?</p>
<p>I'm unsure what I'm doing wrong and would really appreciate some insight! Also, the data has been harmonized already so it should technically work!</p>
<p><a href="https://i.sstatic.net/WIFMj.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/WIFMj.png" alt="Summary Statistics" /></a>
<a href="https://i.sstatic.net/hZrhp.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/hZrhp.png" alt="My data" /></a></p>

CC BY-SA 4.0
2023-12-21T09:19:02.343Z
Please don't put images of text into posts. Copy-paste the actual text and add it as a code block. If formatting is difficult, someone else can edit the post to fix that up.
gringer
<plink><polygenic-risk-score>
Problem with Calculating Polygenic Risk Score (PRS) in Plink
# Problem with Calculating Polygenic Risk Score (PRS) in Plink

(Asked by: username_0 on 2023-12-20)

I've been struggling with constructing a PRS score. In the pics you can see how my summary stats and my target data look like.
Target Data Summary Stats
I've tried the following code which should be pretty straightforward:
plink \
  --bfile target \
  --score summary_stats.txt 1 4 6 header \
  --out final

When I run it, I get this message:
10499994 variants and 1198 people pass filters and QC. Note: No phenotypes present.

The data is binary and the fam file really doesn't have phenotypes (it's -9). Is there any way to still run a PRS?
I'm unsure what I'm doing wrong and would really appreciate some insight! Also, the data has been harmonized already so it should technically work!

**Tags:** plink, polygenic-risk-score
> **Comments:**
> **username_1**: Please don't put images of text into posts. Copy-paste the actual text and add it as a code block. If formatting is difficult, someone else can edit the post to fix that up.  

---

## 1 answers

### Answer by username_2 on 2024-02-29. Score: 0

I'm a little confused about what you mean by there being no phenotype. To build the score you need a phenotype, although to apply it you might not. (Usually PRS is build with a base dataset and then applied using a target dataset). Here's a good tutorial about doing PRS in plink that might help you:

<https://choishingwan.github.io/PRS-Tutorial/plink/>

Also this paper is a really good paper that helped me better understand how building PRS works:

<https://www.nature.com/articles/s41596-020-0353-1>

Hope this helps!


---

