Entering edit mode
4 months ago
Hello,
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:
Rscript PRSice.R \
--dir . \
--prsice ./PRSice_linux \
--base target \
--target summary_stats.txt \
--thread 1 \
--stat BETA \
--binary-target F
plink \
--bfile target \
--score summary_stats.txt 1 4 6 header \
--out final
In both versions, it says no phenotypes found. 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!