Plink 1.9 file-LD based SNP pruning- -problems with autoconvert option??
746 views
Skip to first unread message
Hsarin
unread,
Jul 7, 2015, 7:22:18 AM7/7/15
Reply to author
Sign in to reply to author
Forward
Sign in to forward
Delete
You do not have permission to delete messages in this group
Copy link
Report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to plink2...@googlegroups.com
Hi, I'm trying to perform LD based SNP pruning on my data so I could use the output SNP list on determinating the PCA's from the dataset. The dataset is imputated genotype data.The command I've been trying to use is
In result, I however get temporary .bed .bim and .fam files(and also the .log file) but not the result files : prune.in or prune.out. Could somebody help with this what might be the problem? Really appreciate the answers!
Hsarin
unread,
Jul 7, 2015, 8:19:46 AM7/7/15
Reply to author
Sign in to reply to author
Forward
Sign in to forward
Delete
You do not have permission to delete messages in this group
Copy link
Report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to plink2...@googlegroups.com
The log file says: Total genotyping rate is 0.982489.
489301 variants and 2660 people pass filters and QC.
Error: Too few valid variants for --indep-pairwise.
Christopher Chang
unread,
Jul 7, 2015, 8:29:57 AM7/7/15
Reply to author
Sign in to reply to author
Forward
Sign in to forward
Delete
You do not have permission to delete messages in this group
Copy link
Report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to plink2...@googlegroups.com, heikki...@hotmail.com
--indep-pairwise requires chromosome codes and bp positions to function, since it uses a sliding window. The "chromosome 0 variants" line in the log implies that chromosome codes are missing from the .bgen.
If they're actually present, can you send me the .bgen and .sample files so I can investigate further? Thanks!
Hsarin
unread,
Jul 8, 2015, 8:48:19 AM7/8/15
Reply to author
Sign in to reply to author
Forward
Sign in to forward
Delete
You do not have permission to delete messages in this group
Copy link
Report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to plink2...@googlegroups.com, heikki...@hotmail.com
Ok, so you were right. The chromosome codes are actually missing and also the two first columns are NA and ----- before rsXXXXX. What would be the easiest way to change file to the right form? So I have to delete the first to columns and then add the chromosome codes? Really appreaciate your help!
Daniel Poveda Martínez
unread,
Aug 12, 2021, 7:15:05 PM8/12/21
Reply to author
Sign in to reply to author
Forward
Sign in to forward
Delete
You do not have permission to delete messages in this group
Copy link
Report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message