Comments (5)
Dear Vinay,
could you please type the following command in the R Console on the machine you are using
to perform divergence_stratigraphy() and post the output?
parallel::detectCores()
Depending on the output I will know what the problem might be.
Cheers,
Hajk
from orthologr.
Dear HajkD,
The output I am getting for the above command is [1] 64.
from orthologr.
Dear Vinay,
thanks for the fast response!
I will have a look at it and come back to you.
Cheers,
Hajk
from orthologr.
Dear Hajk,
For unknown reasons the calculation of DI completes in time. Following are my observations which needs your help in interpreting it.
Query sequences number = 49000
Subject sequences number = 33164
Result are obtained for 21474 sequences only.
htop still shows lesser number of threads running, however the job completes.
Does this means that rest of the genes do not have any different?
from orthologr.
Dear Vinay,
Yes, your 20K sequences make absolute sense.
For the divergence stratum assignment first orthologous genes
are inferred using BLAST best reciprocal hit.
This means that if you have Species A and Species B,
first BLAST(A,B) is performed and only the best hit for each CDS is retrieved,
second BLAST(B,A) is performed and again only the best hit for each CDS is retrieved.
Only the CDS that have the best hit in both directions BLAST(A,B) = BLAST(B,A) are
considered as orthologs and are used to compute dNdS values.
So since your smaller organism has 34K seqs the maximum possible number of best hits
would already be 34K, but since of course not all CDS are orthologs, this possibility space
will be reduced to much less than 34K. This is why I think that 20K make absolute sense and
reflect my experience when I am computing dNdS values between species.
I hope this helps you to interpret the results.
Thanks for clarifying the threads issue, I also couldn't find any problems with multicore processing.
All cores are being used.
Best wishes,
Hajk
from orthologr.
Related Issues (20)
- How to set PATH to OrthoFinder2 correctly? HOT 5
- Control output value of orthologs function HOT 1
- Please check the correct path to Comeron... the interface call did not work properly HOT 4
- File blastresult_AllColdShockProteins.fasta.csv could not be read correctly. HOT 3
- Installation error HOT 9
- Ensembl gff does not work HOT 1
- "path" option for Blast calls is non-functional HOT 1
- Clarification on pre-processing data and data types for dNdS estimation HOT 8
- From orthologs to multiple align HOT 8
- ERROR: number of input seqs differ (aa: 1; nuc: 2)!! HOT 1
- Missing headings blast_rec HOT 7
- Missing sequences even with delete_corrupt_cds = FALSE HOT 8
- unable to install on R4.2 HOT 1
- Integrate blastx and tblastx options HOT 2
- caught segfault, address cause 'invalid permissions' HOT 1
- Issue running blastp in dnds HOT 2
- h(simpleErrpr(msg, call)) error HOT 4
- Docker image request HOT 1
- Package release request HOT 4
- Dependency question HOT 4
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 orthologr.