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
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 

/gpfs/homes/hsarin/PLINK_1.9/plink1.9 --bgen ../dilgom_1000g_chr1.bgen –sample ../dilgom_1000g_chr1.sample --indep-pairwise 100 20 0.1 –out SNPs_pruned

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
to plink2...@googlegroups.com
The log file says: Total genotyping rate is 0.982489.
489301 variants and 2660 people pass filters and QC.
Note: No phenotypes present.
--indep-pairwise: Ignoring 489301 chromosome 0 variants.
Error: Too few valid variants for --indep-pairwise.

Christopher Chang

unread,
Jul 7, 2015, 8:29:57 AM7/7/15
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
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
to plink2-users
did you fix this? I also have the same problem

--indep: Ignoring 92460 chromosome 0 variants.
Error: Too few valid variants for --indep.


Reply all
Reply to author
Forward
0 new messages