Thank you very much zx8754 :) it works.
• 0 views
•
link
written
by
aarvaBioinfo •
Sorting of SNPs which has less than 0.05 threshold P-Value
How can we get total number of SNPs in plink association output file as "qassoc" and "adjusted file" which has less than 0.05 P values?
• 2,607 views
•
link
updated
by
zx8754
•
written
by
aarvaBioinfo •
•
written
by
aarvaBioinfo •
1 answer
From the plink manuals the pvalue column for qassoc file is the 9th. So we can use awk to keep the 1st row header, then filter on 9th column value, see example:
awk '(NR==1) || ($9 < 0.05) ' myfile.qassoc > myfile_subset.qassoc
To get the number of rows:
awk '$9 < 0.05' myfile.qassoc | wc -l
Hi aarvaBioinfo,
If an answer was helpful you should upvote it, if the answer resolved your question you should mark it as accepted.
Log in to answer this question.
More posts like this
-
p-adjusted threshold in pathfindR
written by mrashad
-
PLINK meta-analysis - multiple parameters after '+'
written by tarumuranen •
-
How to reduce number of entries in QQ plot?
written by anamaria
-
Get SNP information with association test with PLINK
written by landscape95
-
what to do when adjusted p-values are not acceptable in RNASeq?
written by nazaninhoseinkhan
-
Extracting row and column names from matrix in R.
written by kandoigaurav
-
Quantitative trait analysis
written by toyan.j.p •
-
Logistic regression with PLINK, p-values
written by monsegd14 •
-
Plink And Missing Values
written by RawRaw •
-
Clarification On P-Value Threshold For Gwas
written by mlscmahe
please post some example lines from that files.