Question: R, PopGenome Package. Cannot open file
0
gravatar for pifferdavide
4.2 years ago by
pifferdavide90
Italy
pifferdavide90 wrote:

I downloaded the Chromosome in vcf format along with the tbi file from 1000 Genomes ftp server. I installed the packaged "PopGenome". However, I cannot open the files. I type this command and get an error message. Am I missing something? 

readVCF("C:/Users/Davide/Documents/chr22/Chr22.vcf", 1000,"22",1,100000)
open: No such file or directory
Caught exception inside whop_tabix::open('C:/Users/Davide/Documents/chr22/Chr22.vcf.gz"):
    'whop_tabix::open : Failed to open tabix index file'
return FALSE from whoptabix_open
vcff::open : could not open tabix-index!
VCF_open : Could not open file 'C:/Users/Davide/Documents/chr22/Chr22.vcf' as tabix-indexed!
[1] FALSE

 

snp popgenome R • 2.4k views
ADD COMMENTlink modified 3.9 years ago by bastianxpfeifer50 • written 4.2 years ago by pifferdavide90

Did you try this: GENOME.class=readVCF() instead of just readVCF()?

By the way, the example .VCF file didn't have a tabix-index file with it, why it can load in successfully?

 

ADD REPLYlink written 4.1 years ago by huan192010
0
gravatar for RamRS
4.2 years ago by
RamRS20k
Houston, TX
RamRS20k wrote:

You should have a Chr22.vcf.tbi tabix index in the same folder as the VCF file. If that doesn't exist, you can always create one as shown here: http://samtools.sourceforge.net/tabix.shtml

ADD COMMENTlink written 4.2 years ago by RamRS20k

I do have a Chr22.vcf.tbi tabix index in the same folder.

ADD REPLYlink written 4.2 years ago by pifferdavide90

Can you check the permissions on the file please? Also Windows might set security flags on downloaded files that you might have to "unblock".

ADD REPLYlink written 4.2 years ago by RamRS20k

You were indeed right but it's still not working, even after I unblocked both files...still getting same error message!

ADD REPLYlink written 4.2 years ago by pifferdavide90

Can you try, on R:

Sys.chmod("C:/Users/Davide/Documents/chr22/Chr22.vcf", mode = "0777",showWarnings=TRUE)
Sys.chmod("C:/Users/Davide/Documents/chr22/Chr22.vcf.tbi", mode = "0777",showWarnings=TRUE)
ADD REPLYlink modified 4.2 years ago • written 4.2 years ago by RamRS20k

I tired it but I get no output

ADD REPLYlink written 4.2 years ago by pifferdavide90
0
gravatar for pifferdavide
4.2 years ago by
pifferdavide90
Italy
pifferdavide90 wrote:

I have both files in the same folder and I downloaded from here: ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/

ADD COMMENTlink written 4.2 years ago by pifferdavide90
0
gravatar for bastianxpfeifer
3.9 years ago by
Germany
bastianxpfeifer50 wrote:

Hi ALL, see also http://popgenome.weebly.com/

THX, all the best.

Bastian

ADD COMMENTlink written 3.9 years ago by bastianxpfeifer50
Please log in to add an answer.

Help
Access

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