Hello All, I want to convert TCGA MAF to VCF analysis with PLINK I worked on the conversion of MAF2VCF using maf2vcf.pl, I got converted VCF file, but while conversion it has missing phenotype data and some gene id. due to this I could not able to process with PLINK analysis. please suggest me to exact conversation or adding phenotype data. Thank you.
To better help, you should show all commands that you have used.
Some general advice: when converting a VCF to PLINK format, PLINK will not know what are the phenotypes. You need to create a custom FAM file that matches exactly to the order of samples in your PLINK object - you can then use this FAM file in all downstream analysis via the --fam command line parameter.
Please note: PLINK will NOT store your samples in the same order as they appear in the VCF. You can control this, however, via the --indiv-sort file parameter when converting. Please read my previous answer, here: linkage disequilibrium analysis
Kevin
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
version 2.3.6
Hello Sir, Thanks for the reply and suggestion. the below Perl command I used to convert maf to vcf
it has converted as
TCGA.HNSC.varscan.5296cf00-4d8c-4db3-80d7-930a4b44f90d.DR-10.0.somatic.vcfand tsv smaple file likeNow i want to custom .fam file how i can use this .tsv sample file to generate.fam file indivally please let me know best way. If i generate .fam custom is it have any problem internal connection with .bim .bam to process further work.
You will generate the FAM file on your own, by studying how it is structured: https://www.cog-genomics.org/plink/1.9/formats#fam
Still, I am not getting to generate .fam file because I don't have patient family details, I have only gender phenotype values. I am new to this work not getting idea could you explain how I can process this work. Thanks.
Still, I am not getting to generate .fam file because I don't have patient family details, I have only gender phenotype values. I am new to this work not getting idea could you explain how I can process this work. Thanks.
You do not have to have FID (family ID) in order to use PLINK. Also, this process is simply not something that is easy to explain across a web forum, like Biostars. I have provided enough details already in my previous answer: linkage disequilibrium analysis
Hi,
I am new to this field and I have a problem with converting MAF file to VCF format with perl command. Is the maf2vcf.pl still valid and how did you download and use this file to convert your MAF file to VCF format? Any help is really appreciated!