Folks,
I'm trying to get ethnic-specific MAF for a set of SNPs inside R (relatively novice R user).
So, I searched BioConductor's AnnotationHub for human, dbsnp, etc and found AH47013.
However, when I tried to load it
library(AnnotationHub)
ah = AnnotationHub()
# launch interactive dataset query interface
d <- display(ah)
AH47013 <- ah[["AH47013"]] # load dbSNP 141
warnings()
that showed that it was failing to load https://annotationhub.bioconductor.org/fetch/52464
which is actually a redirect to
ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606_b141_GRCh37p13/VCF/common_all.vcf.gz
which no longer exists and has been replaced with dbSNP144:
ftp://ftp.ncbi.nih.gov/snp/organisms/human_9606_b144_GRCh37p13/VCF/common_all_20150605.vcf.gz
So, the questions are :
- how does AnnotationHub get updated?
- how do I work around this in the mean time?
- was there a different resource or a better way to get this done in BioConductor?
- (yes, we've found SNPsnap server @ Broad, but want to do this locally in R)
Thanks!
Curtis
This is a good question and I also tried to use AnnotationHub without success. I believe this should be done through the SNPLocs data packages, however the most recent package has not been updated since dnSNP144 (we are now at dbSNP146).