readSnpMatrix in FACETS gives error
0
0
Entering edit mode
6.5 years ago
hsbinf ▴ 30

Hi! I'm trying to run FACETS (allele-specific copy number and clonal heterogeneity analysis tool for high-throughput DNA sequencing).

I created a SNP matrix using the perl script snp-pileup.pl. But when I run it in R it gives an error.

rcmat <- readSnpMatrix("sample1_facets_input", skip=0L, err.thresh=Inf, del.thresh=Inf, perl.pileup=TRUE) 
Error in scan(filename, what = list(Chromosome = "", Position = 0, NOR.DP = 0,  : 
scan() expected 'a real', got 'chr1'

Any ideas why? I can run most of the samples, its just 1 sample that has this error.

Thanks!

CNV Purity Ploidy FACETS • 1.5k views
ADD COMMENT
0
Entering edit mode

Just a guess: in that particular sample, your chromosomes have the 'chr' prefix, when they should be just encoded numerically,. Strip the 'chr' prefix from the input file and then retry.

Kevin

ADD REPLY

Login before adding your answer.

Traffic: 2066 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