Comments (1)
The scores are approximate log-likelihoods (with the minus sign omitted). They are difficult to interpret on their own and are more useful when comparing models. In theory, an absolute log-likelihood score difference of one suggests that the probability of the data under one model is about 2.7 times as high as under the other model. However, the likelihood (probability of the data given the model) is very different from the posterior probability (probability of the model given the data).
Another way of looking at the score is to consider that it is calculated as a weighted sum of f-statistic z-scores residuals (difference in z-scores of fitted minus observed f-statistics). A score of 694 suggests that one or more f-statistic z-scores are much greater than 4, so there is strong evidence that at least one observed f-statistic is not compatible with this model. But that is assuming that the goal is to find a model which fits the data perfectly. The problem with that is that in the limit of infinite data, any minor deviations between the real history and the model will lead to large f-statistic z-score residuals and large absolute scores. At the same time, noisy data and large models with too many degrees of freedom will result in "better" scores, closer to zero, even if these models have nothing to do with the actual history.
So I would focus more on the difference in scores of competing models, not on the absolute score of a single model.
from admixtools.
Related Issues (20)
- Optimization unsuccessful (constraints inconsistent). Try increasing 'diag', when running qpgraph with a graph and upper and lower bounds HOT 2
- Error in parse_qpadm_output HOT 1
- The order of pops in outgroup f3 test HOT 2
- Explanation of the functions qpdstat and qpf4ratio HOT 1
- Protocol Advice: Establishing a t hresholds for the number of allowed admixture events HOT 3
- qp3pop unique_only does not work as expected HOT 2
- Do qpWave, qpAdm, and f5 have requirements on the minimum number of individuals for each population? HOT 1
- est_to_boo does not preserve the SNP block names (number of SNPs), but est_to_loo() does. Is that intended? HOT 2
- large difference of f3 result between admixtools and admixtools2 HOT 1
- Zero drift edges HOT 2
- Behavior affected by other R packages? HOT 2
- Error in if (gimp > 0 && gimp%%plusminus_generations == 0) { : missing value where TRUE/FALSE needed HOT 2
- please add tags/releases HOT 2
- Inconsistent bootstrap significance testing HOT 5
- keyword argument typo lack of error HOT 3
- Comparing graphs with compare_fits HOT 2
- auto_only: change default value to FALSE HOT 1
- qpWave returning one fewer ranks than expected HOT 1
- qpAdm computation stucks when the num. of letf pops below 3
- Issue with f3 and f4 HOT 2
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 admixtools.