oceanic_computational_asr's People
oceanic_computational_asr's Issues
dialect match
dialect - dialect match for akla1240. What's up?
0 branch lengths in Gray et al 2009 MCCT-tree
I have a problem when running phytools::make.simmap() with the Gray et al 2009 MCCT-tree. make.simmap()
throws this error:
Error in 1:nrow(L) : argument of length 0
For full clarity, I'm running:
phytools::make.simmap(
model = "ARD",
pi = "estimated",
method = "optim")
I am fairly confident I've diagnosed the issue, there are branch lengths of length zero in the tree. The 0 branch lengths are there already in the MCCT-tree, they're not the result of pruning or some other manipulation of the tree.
In this project, I'm doing ASR with parsimony (castor::asr_max_parsimony()
, script here), ML (corHMM::corHMM()
, script here) and SCM (phytools::make.simmap()
, script here). This problem with the 0 branch lengths only stops the analysis with make.simmap()
, meaning the other two functions are eating up the tree just fine - but that I should probably still be concerned since there are still 0 branch lengths in all analysis with this tree.
Now, I know of three ways of solving this:
compute.brlen()
- replacing the 0 branch lengths with something tiny, but not zero. ( for example:
tree$edge.length[tree$edge.length==0] <- max(nodeHeights(tree ))*1e-6
) - sampling the posterior
I don't want to do (1) because I don't want Grafen branch lengths when I have actual branch lengths to work with. That option is scrapped right away.
The choice is between (2) and (3). I checked all the posterior trees, and as far as I can tell none of them have 0 branch lengths. But, I'm not 100% confident about the way I was checking them so if someone else has a manner of checking them and finds a number greater than 0, holler.
The argument against doing (3) is mainly that this is a project that's aimed at being understood by traditional linguists, and I think they would struggle a bit with the idea that I'm randomly choosing a tree out of a posterior each time I'm doing the ASR and comparing to their findings.
The argument against (2) is that even teeny-tiny branch lengths will do weird things with SCM, since trying to model change along such very very short branches makes for weird results.
I'm leaning towards doing (3). The question then is when to randomly sample:
a) every time the ASR is run on each feature (once per feature + methods combo, i.e. 201 *3)
b) once for all ASR on all features per method (once for parsimony, once for ML and once for SCM, i.e 3 times)
c) once for all ASR on all features on all methods (once in total)
(a) seems the methodologically soundest to me.
Advice @king-ben @SimonGreenhill?
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.