Comments (4)
Hi YangChao,
The tree you get out of ClonalFrameML is not the same as the tree you would get by removing the recombinant regions.
The lengths of the branches in the ClonalFrameML output tree take into account both mutation and recombination events that have been found on that branch. If a lot of recombination events have been detected on a given branch, excluding all these regions would result in a very short branch, but not so in ClonalFrameML, and indeed if we see a lot of recombination events it must mean that a significant amount of evolution has taken place and so the branch can't be that short.
Also, excluding recombinant regions before building the tree would mean that recombination regions are excluded for all branches. So if you have that all sites are recombinant on at least one branch, you would be left with nothing to build your tree, whereas ClonalFrameML accounts for recombination and the branches where they were inferred.
If you really want to build a tree based on only the non-recombinant regions then you can excluded the regions in importation_status and build a tree from the remainder, but for the reasons described above the ClonalFrameML output tree should be closer to the truth.
Best wishes,
Xavier
from clonalframeml.
Hi Xavier,
If ClonalFrameML output tree take into account both mutation and recombination events, how to calculate the length of each branch? For example, we see 10 SNPs in a branch, 5 SNPs was introduced by 5 mutations, the other 5 SNPs were introduced by 1 recombination. In the ClonalFrameML output tree, how many SNPs would present? How to correct branch length in ClonalFrameML?
Best wishes,
Yangchao
from clonalframeml.
The length of each branch is calculated as a function of the number of mutation and recombination events found on that branch, and the rates at which mutation and recombination events happen. Cf the ClonalFrameML paper for a more exact description.
from clonalframeml.
Thank you so much.
from clonalframeml.
Related Issues (20)
- cfml_results.R error HOT 8
- xmfa format question; error: Segmentation fault (core dumped) HOT 1
- Question: xmfa format? HOT 7
- Recombination events total HOT 1
- convert_NewickTree_to_marginal_tree() HOT 2
- ERROR: incompatible with a bifurcating unrooted tree HOT 1
- node labels HOT 2
- ClonalFrameML for analyses of gene families HOT 2
- How to solve problems about "makefile", Thanks a lot! HOT 2
- input alignment for ClonalFrameML HOT 1
- segmentation fault (core dumped) HOT 1
- how to calculate the r/m value for subclade of a tree๏ผ HOT 4
- R values for calculation of R*delta HOT 2
- Recombination filtered output vs importation_status.txt HOT 1
- why the r/m I calculated is so small using a subclone of subclone of klebsiella pneumoniae HOT 1
- Is it better to order genes by the reference genome for a xmfa file ??? HOT 2
- std::bad_alloc problem HOT 1
- Are there any detailed description of the output pdf from cfml_results.R? HOT 1
- ERROR: Token: ((( Found left bracket but no right bracket HOT 1
- High recombination sequence detection HOT 3
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 clonalframeml.