R, package: ldblock, function stack1kg error in opening BCFFile
0
0
Entering edit mode
ting835698 • 0
@534bda6c
Last seen 3.2 years ago

I am trying to get a VCF stack from the 1000 genome project using package "ldblock" on R (version 4.0.3), but when I used function

sta = stack1kg()

the following error was shown.

Error in open.BcfFile(BcfFile(file, character(0), ...)) : 'open' VCF/BCF failed filename: http://ftp.1000genomes.ebi.ac.uk/vol1/ftp/data_collections/1000_genomes_project/release/20190312_biallelic_SNV_and_INDEL/ALL.chr1.shapeit2_integrated_snvindels_v2a_27022019.GRCh38.phased.vcf.gz

I tried stack1kg(chrs = as.character(1:22)) and manually change the number from 1 to 22, all shown similar message as above (the only thing that changed was the filenames). So I assumed it is caused by missing package(s). I downloaded Rsamtools, VariantAnnotation, updated xfun, GenomicFiles and some others, none seems to work.

I dont know whether the material I am learning from is a bit too old (the last update is in 2017) , does everyone still using ldblock to obtain VCF stacks? If anyone had experienced something similar to this or if there is a solution pls let me know. Many Thanks.

ldblock • 734 views
ADD COMMENT
0
Entering edit mode

Hi, I'm facing the same problem right now. Did you find any solutions? Best, Mehdi

ADD REPLY

Login before adding your answer.

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