Question: SNP IDs changed when making .ped files using PLINK
gravatar for kianalee
6 months ago by
kianalee0 wrote:

I have data on genetic variants for multiple individuals for a family in the form of a .vcf file. I am successfully able to convert this .vcf file into a .ped file using plink:

plink --vcf  myvcf file --recode --out myplink

But when I read the .ped file into R and look at $genotype or $map, the SNPs are not named by either chromosomal location or an rs number, but are rather labeled "locus.1", "loucs.2", etc.

For example, in R:

>head(sample$map) allele.1 allele.2
locus.1  locus.1        G        A
locus.2  locus.2        G        A
locus.3  locus.3        C        T
locus.4  locus.4        A     <NA>
locus.5  locus.5        C    CCCCT
locus.6  locus.6        T        G

How can I either trace back the actual genomic position of these loci, or make .ped files with either the chromosomal location or rs number?

snp plink • 222 views
ADD COMMENTlink modified 6 months ago by chrchang5235.7k • written 6 months ago by kianalee0
gravatar for chrchang523
6 months ago by
United States
chrchang5235.7k wrote:

If the chromosomal locations and rsIDs are in your VCF file, they'll also be in the .ped file generated by your plink command. I'd double-check how the R package works.

ADD COMMENTlink written 6 months ago by chrchang5235.7k
Please log in to add an answer.


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