This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Writing PLINK PED and MAP files... Error: Could not open temporary file.

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

format plink

1 answer

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

Log in to answer this question.