Question: Writing PLINK PED and MAP files... Error: Could not open temporary file.
1
gravatar for angelaparody
2.1 years ago by
angelaparody50
angelaparody50 wrote:

Hi,

I have edited #CHROM of my .vcf file so it is purely numeric (from 1 to 13635, which correspond to scaffolds#). I need this to avoid a error message on a downstream step. I need to get .ped/.fam format and using plink2 it doesn't allow me to keep these CHROM notation (I need to use --allow-extra-chr, and then it changes all CHROM into 0). Then I found out that vcftools can handle this, but it only works when I convert it into transpose ped/fam (tped/tfam). See below when I try to convert to .ped/.map (unsuccessfully) and how it worked when using --plink-tped:

MacBook-Pro-de-Angela:UsingVCFtools angelaparodymerino$ vcftools --vcf mac3_minDP3_maxmeanDP289_maf005_minQ40_minGQ30_hwe005_265ind_IDs2_CHROM.vcf --out plinkvcf --plink

VCFtools - 0.1.15
(C) Adam Auton and Anthony Marcketta 2009

Parameters as interpreted:
    --vcf mac3_minDP3_maxmeanDP289_maf005_minQ40_minGQ30_hwe005_265ind_IDs2_CHROM.vcf
    --out plinkvcf
    --plink

After filtering, kept 265 out of 265 Individuals
Writing PLINK PED and MAP files ... 
Error:  Could not open temporary file.

MacBook-Pro-de-Angela:UsingVCFtools angelaparodymerino$ vcftools --vcf mac3_minDP3_maxmeanDP289_maf005_minQ40_minGQ30_hwe005_265ind_IDs2_CHROM.vcf --out plinkvcf --plink-tped

VCFtools - 0.1.15
(C) Adam Auton and Anthony Marcketta 2009

Parameters as interpreted:
    --vcf mac3_minDP3_maxmeanDP289_maf005_minQ40_minGQ30_hwe005_265ind_IDs2_CHROM.vcf
    --out plinkvcf
    --plink-tped

After filtering, kept 265 out of 265 Individuals
Writing PLINK TPED file ... Writing PLINK TFAM file ... Done.
After filtering, kept 3412 out of a possible 3412 Sites
Run Time = 1.00 seconds

However, I need ped/map (not tped/tmap, because it is not an accepted input file in Mega2). The error message says that: "Could not open temporary file". Anyone can help? In short, what I need is to convert my .vcf file into .ped/.map but keeping the CHROM as it is in my .vcf file (numeric and from 1 to 13635). In general I am finding problematic to work without real chromosome information (just with scaffolds) in most of the steps/conversions.

Thanks in advance,

Regards,

Ángela Parody-Merino

plink format • 2.3k views
ADD COMMENTlink modified 9 months ago by zx87548.4k • written 2.1 years ago by angelaparody50
0
gravatar for Zhilong Jia
2.1 years ago by
Zhilong Jia1.5k
London
Zhilong Jia1.5k wrote:

A related post: Vcftools Can'T Create Ped Map File Due To Ulimit Of Allowed Open Files.

ADD COMMENTlink written 2.1 years ago by Zhilong Jia1.5k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2098 users visited in the last hour