snpStats:::read.pedfile won't read in file.map for snps =
0
0
Entering edit mode
5.0 years ago
kianalee ▴ 10

I am having problems with reading in .ped files using snpsStats::read.pedfile.

I try to read in the pedfile into R:

>sampleA <- snpStats::read.pedfile(file = "sampleA.ped", snps = "sampleA.map")

I used this before on a different pedfile and was able to identify the SNPs (snp.names) by RS#, however, now I receive the following error:

Error in snpStats::read.pedfile(file = "sampleA.ped", snps = "sampleA.map") :    
No unambiguous snp names found on file

If I just read it in with

>sampleA <- snpStats::read.pedfile(file = "sampleA.ped")
>head(sampleA$map)

I get:

snp.name allele.1 allele.2
locus.1  locus.1        G     <NA>
locus.2  locus.2        G     <NA>
locus.3  locus.3        C        T
locus.4  locus.4        A     <NA>
locus.5  locus.5    CCCCT        C
locus.6  locus.6        G        T

Which doesn't identify the snp.names as RS #'s.

Does anyone have any suggestions on how to fix this?

read.pedfile genome R • 1.6k views
ADD COMMENT

Login before adding your answer.

Traffic: 2771 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6