Go Back   SEQanswers > General

Similar Threads
Thread Thread Starter Forum Replies Last Post
Converting VCF to PLINK PED format, need help? ketan_bnf Bioinformatics 2 03-21-2017 02:35 AM
creating binary ped files Rabu Bioinformatics 0 11-13-2014 10:02 AM
MISO output - mismatch between number of transcripts and number of assigned counts MikeStubbington Bioinformatics 1 01-29-2014 08:44 AM
Splitting a BAM file by every x number of reads, not lines. carmeyeii Bioinformatics 10 05-06-2013 12:02 PM
BWA aln _ option for number of mismatch ??? Sam64 Bioinformatics 2 10-25-2011 12:59 AM

Thread Tools
Old 04-10-2015, 03:51 AM   #1
Junior Member
Location: Oslo, Norway

Join Date: Apr 2015
Posts: 1
Default Mismatch between number of lines in plink .ped files and corresponding .map file

I have a problem with a mismatch between the number of columns in a plink .ped file and the corresponding map file; when I try to produce a .bed file I get the error message below.

I would be really really grateful for any input on this!

17028030 (of 17028030) markers to be included from [ ]

A problem with line 1 in [ raw_phased.ped ]
Expecting 6 + 2 * 17028030 = 34056066 columns, but found 20352684

The ped.file has the first 6 columns and they seem correct. The map file seems fine too. Code used below:

I have used an unfiltered vcf file to make sure that it is not due to filtering issues, and phased it using beagle:

module load beagle

java -Xmx4000m -jar /projects/cees/bin/beagle-4.0/b4.r1274.jar gt=../islands_parents_treeHM60_1.vcf ibd=true out=raw_phased nthreads=16 phase-its=10

Then I use vcftools to produce the plink .ped and map files

vcftools --gzvcf raw_phased.vcf.gz --plink --out raw_phased

Then the problem arises when I try to produce .bed and .raw files using plink
(version 1.07, the only one on the cluster, so that is how it has to be). In these steps no files are produced and the error message on top is displayed in the plink.log

#making bed file with plink
plink --file raw_phased --make-bed

#produce a .raw file that can be read into R using read.table
plink --recodeA --file raw_phased
AnnaRu is offline   Reply With Quote

ped-files, plink, raw-file

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 02:59 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO