Comments (1)
This is a bit late, but I think you raised an interesting question.
What you call linearity (Dac <= (Dab + Dbc)) is sometimes called triangle inequality, and you are right in that both f2 and (∑|aj-bj| )/M have that property (regardless of how a, b, and c are related to each other), whereas Fst doesn't have it.
f2 doesn't have the proportionality property as you define it, but it has another property which is perhaps more important, and which (∑|aj-bj| )/M doesn't satisfy: If a, b, and c are related like a -> b -> c or like a <- b -> c then f2(a, c) = f2(a, b) + f2(b, c).
This is because f2 is the variance of allele frequency differences between two populations, and the sum of the variances of two independent random variables (drift between a and b and drift between b and c) is equal to the variance of the sum of both (drift between a and c).
See Ancient Admixture in Human History, section Additivity of F2 along a tree branch for more details.
from admixtools.
Related Issues (20)
- 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
- running many replicates of find_graphs? HOT 1
- Issue with extract_f2 HOT 2
- Problem making f2 file from plink .bed .fam and .bim files HOT 6
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.