Hello I use following command to get 012 matrix of vcf.gz
plink --vcf test.vcf.gz --remove ID_of_miss_indv_for_plink --recode A --out TEST
But I only want genotypes of SNPs (I want to keep the columns rs149437228_C, rs179827_A). I do not want to plink print the columns with names "._C" "._G" ,"._AC" and so on.
Can some one help me?
plink --vcf test.vcf.gz --remove ID_of_miss_indv_for_plink --recode A --out TEST
But I only want genotypes of SNPs (I want to keep the columns rs149437228_C, rs179827_A). I do not want to plink print the columns with names "._C" "._G" ,"._AC" and so on.
Can some one help me?
Comment