Comments (4)
Hi!
For taxonomic groups with large amounts of genetic data the download phase can be very slow and sometimes errors can arise due to records getting malformed (in your case it looks like a record has lost its description). Sometimes simply trying to download on a different day helps.
I have a few suggestions:
- Try restarting the
phylotaR
pipeline usingrestart()
- Try updating to the latest version of
phylotaR
(see below code snippet) and re-running. - Use
restez
. This package downloads large chunks of GenBank to your local machine andphylotaR
can query that instead of the remote database. It should be much faster and reduces the risk of malformed records. Downloading withrestez
is much faster because it fetches highly compressed sequence files. One downside is you may need to download large chunks of sequence data, you may need a lot of free disk space.
Currently restez
is not available via CRAN due to a key dependency being archived. To use restez
you will need to install packages via GitHub.
# devtools required for installtion via GitHub
install.packages("devtools")
# install key dependency no longer available on CRAN
devtools::install_github("hannesmuehleisen/MonetDBLite-R")
# install restez
devtools::install_github("ropensci/restez")
# install latest version of phylotaR (the one on CRAN does not make use of restez)
devtools::install_github("ropensci/phylotaR")
from phylotar.
Thank you very much for your answer. However, now when using PhylotaR (after restez download), I have the error: Error in stages_run(wd = wd, frm = 1, to = nstages, stgs_msg = stgs_msg) :
Unexpected Error in typ_location[[2]] : subscript out of bounds.
I attached the code used
library(phylotaR)
library(restez)
local_genbank<-getwd()
restez_path_set(local_genbank)
db_download(db = 'nucleotide')
db_create(db = 'nucleotide')
###phylotaR
ncbi_dr <- "/home/evolecolab/Software/ncbi-blast-2.9.0+-x64-linux/ncbi-blast-2.9.0+/bin"
restez_path_set(filepath = local_genbank)
txid <- 9443
setup(wd = getwd(), txid = txid, ncbi_dr = ncbi_dr)
run(wd = getwd())
from phylotar.
For info, we're looking for a new maintainer / a new maintainer team for this package, see #57 and feel free to volunteer, we'd be happy to help.
from phylotar.
Thank you very much for your answer. However, now when using PhylotaR (after restez download), I have the error: Error in stages_run(wd = wd, frm = 1, to = nstages, stgs_msg = stgs_msg) : Unexpected Error in typ_location[[2]] : subscript out of bounds. I attached the code used
library(phylotaR)
library(restez)
local_genbank<-getwd()
restez_path_set(local_genbank)
db_download(db = 'nucleotide')
db_create(db = 'nucleotide')
###phylotaR
ncbi_dr <- "/home/evolecolab/Software/ncbi-blast-2.9.0+-x64-linux/ncbi-blast-2.9.0+/bin"
restez_path_set(filepath = local_genbank)
txid <- 9443
setup(wd = getwd(), txid = txid, ncbi_dr = ncbi_dr)
run(wd = getwd())
This is similar to #39, it should work in the latest version. If not, please file a new issue.
from phylotar.
Related Issues (20)
- problem with the restez/phylota integration
- Possibility to supply NCBI API key when running pipeline to access NCBI data faster? HOT 1
- Use phylotaR for retrieving orthologous protein sequences HOT 1
- No periods or tildes in file paths HOT 3
- Multiple markers/loci per cluster issue HOT 4
- Can not install the package in windows HOT 2
- Clusters with different parent number but the same MAD, seed, length, n_taxa, n_seqs
- Description - not enough information
- New Maintainer Wanted HOT 5
- Next release HOT 6
- CRAN archive of package outsider
- Fix pkgdown configuration for reference HOT 3
- No need for restez_connect() HOT 1
- Remove treeman as a dependency? HOT 2
- Could get something new from https://github.com/pachterlab/gget
- test failure HOT 1
- Generating taxonomic dictionary error HOT 1
- Fetch and error seem to be stuck HOT 2
- drop taxa functionality HOT 1
- Identical clusters making it into alignments HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from phylotar.