Applications of Analytic Combinatorics in Mathematical Biology (joint with H. Chang, M. Drmota, E. Y. Jin, and Y.-W. Lee)
|
|
- Wilfrid Cox
- 5 years ago
- Views:
Transcription
1 Applications of Analytic Combinatorics in Mathematical Biology (joint with H. Chang, M. Drmota, E. Y. Jin, and Y.-W. Lee) Michael Fuchs Department of Applied Mathematics National Chiao Tung University June 30th, 2015 Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
2 Analytic Combinatorics (i) Combinatorialists use recurrence, generating functions, and such transformations as the Vandermonde convolution; others, to my horror, use contour integrals, differential equations, and other resources of mathematical analysis. - John Riordan (1968). Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
3 Analytic Combinatorics (ii) Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
4 Outline of the Talk 1 Introduction Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
5 Outline of the Talk 1 Introduction 2 Patterns in Phylogenetic Trees H. Chang and M. Fuchs (2010). Limit theorems for patterns in phylogenetic trees, J. Math. Biol., 60:4, Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
6 Outline of the Talk 1 Introduction 2 Patterns in Phylogenetic Trees H. Chang and M. Fuchs (2010). Limit theorems for patterns in phylogenetic trees, J. Math. Biol., 60:4, Shapley Value and Fair Proportion Index M. Fuchs and E. Y. Jin (2015). Equality of Shapley value and fair proportion index in phylogenetic trees, J. Math. Biol., in press. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
7 Outline of the Talk 1 Introduction 2 Patterns in Phylogenetic Trees H. Chang and M. Fuchs (2010). Limit theorems for patterns in phylogenetic trees, J. Math. Biol., 60:4, Shapley Value and Fair Proportion Index M. Fuchs and E. Y. Jin (2015). Equality of Shapley value and fair proportion index in phylogenetic trees, J. Math. Biol., in press. 4 Number of Groups formed by Social Animals M. Drmota, M. Fuchs, Y.-W. Lee (2016+). Stochastic analysis of the extra clustering model for animal grouping, in revision. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
8 Evolutionary Biology Charles Darwin ( ) Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
9 Evolutionary Biology First notebook on Transmutation of Species (1837) Charles Darwin ( ) Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
10 Phylogenetic trees (=PTs) Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
11 Phylogenetic trees (=PTs) Phylogenetic tree of size n: rooted, plane, unlabelled binary tree with n external nodes (and consequently n 1 internal nodes). Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
12 Applications of PTs Understanding genetic relatedness of species Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
13 Applications of PTs Understanding genetic relatedness of species Understanding the underlying evolutionary process Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
14 Applications of PTs Understanding genetic relatedness of species Understanding the underlying evolutionary process Predicting possible future outcomes Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
15 Applications of PTs Understanding genetic relatedness of species Understanding the underlying evolutionary process Predicting possible future outcomes Testing appropriateness of random models Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
16 Applications of PTs Understanding genetic relatedness of species Understanding the underlying evolutionary process Predicting possible future outcomes Testing appropriateness of random models Making conservation decisions in genetics Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
17 Applications of PTs Understanding genetic relatedness of species Understanding the underlying evolutionary process Predicting possible future outcomes Testing appropriateness of random models Making conservation decisions in genetics Modeling the group formation process of social animals Etc. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
18 Yule-Harding model or Kingman Coalescent Example: time species Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
19 Yule-Harding model or Kingman Coalescent Example: Random Model 1: At every time point, two yellow nodes uniformly coalescent. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
20 Yule-Harding model or Kingman Coalescent Example: Random Model 1: At every time point, two yellow nodes uniformly coalescent. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
21 Yule-Harding model or Kingman Coalescent Example: Random Model 1: At every time point, two yellow nodes uniformly coalescent. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
22 Yule-Harding model or Kingman Coalescent Example: Random Model 1: At every time point, two yellow nodes uniformly coalescent. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
23 Yule-Harding model or Kingman Coalescent Example: Random Model 1: At every time point, two yellow nodes uniformly coalescent. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
24 Yule-Harding model or Kingman Coalescent Example: Random Model 1: At every time point, two yellow nodes uniformly coalescent. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
25 Yule-Harding model or Increasing Binary Trees Example: Random Model 2: At every time point, a yellow node is replaced by a cherry. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
26 Yule-Harding model or Increasing Binary Trees Example: Random Model 2: At every time point, a yellow node is replaced by a cherry. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
27 Yule-Harding model or Increasing Binary Trees Example: Random Model 2: At every time point, a yellow node is replaced by a cherry. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
28 Yule-Harding model or Increasing Binary Trees Example: Random Model 2: At every time point, a yellow node is replaced by a cherry. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
29 Yule-Harding model or Increasing Binary Trees Example: Random Model 2: At every time point, a yellow node is replaced by a cherry. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
30 Yule-Harding model or Increasing Binary Trees Example: Random Model 2: At every time point, a yellow node is replaced by a cherry. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
31 Yule-Harding model or Increasing Binary Trees Example: Random Model 2: At every time point, a yellow node is replaced by a cherry. Random model 1 and random model 2 are the same. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
32 Binary Search Trees Example: Input 1, 4, 2, 5, 3 Random Model 3: All permutations of the input are equally equally. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
33 Binary Search Trees Example: Input 1, 4, 2, 5, 3 1 Random Model 3: All permutations of the input are equally equally. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
34 Binary Search Trees Example: Input 1, 4, 2, 5, Random Model 3: All permutations of the input are equally equally. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
35 Binary Search Trees Example: Input 1, 4, 2, 5, Random Model 3: All permutations of the input are equally equally. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
36 Binary Search Trees Example: Input 1, 4, 2, 5, 3 1 Random Model 3: All permutations of the input are equally equally. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
37 Binary Search Trees Example: Input 1, 4, 2, 5, 3 1 Random Model 3: All permutations of the input are equally equally. 3 Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
38 Binary Search Trees Example: Input 1, 4, 2, 5, 3 1 Random Model 3: All permutations of the input are equally equally. 3 Again the same as random model 1 and random model 2. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
39 Patterns in PTs k-pronged node (McKenzie and Steel 2000; Rosenberg 2006): Node with an induced subtree of size k. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
40 Patterns in PTs k-pronged node (McKenzie and Steel 2000; Rosenberg 2006): Node with an induced subtree of size k. k-caterpillar (Rosenberg 2006): Induced subtree of size k with an internal node which is descendent of all other internal nodes. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
41 Patterns in PTs k-pronged node (McKenzie and Steel 2000; Rosenberg 2006): Node with an induced subtree of size k. k-caterpillar (Rosenberg 2006): Induced subtree of size k with an internal node which is descendent of all other internal nodes. Node with minimal clade size k 3 (Blum and François (2005)): Node with induced subtree of size k and either right or left subtree is an external node. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
42 Mean and variance of k-pronged nodes X n,k = # of k-pronged nodes in random PT of size n. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
43 Mean and variance of k-pronged nodes X n,k = # of k-pronged nodes in random PT of size n. Rosenberg 2006: we have µ n,k := E(X n,k ) = 2n, (n > k) k(k + 1) and 2(4k 2 3k 4)(k 1)n k(k + 1) 2, if n > 2k; (2k 1)(2k + 1) σn,k 2 := Var(X 2(5k 7)(k 1) n,k) = (k + 1) 2, if n = 2k; (2k 1) 2(k 2 + k 2n)n k 2 (k + 1) 2, if 2k > n > k. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
44 Mean and variance of k-pronged nodes X n,k = # of k-pronged nodes in random PT of size n. Rosenberg 2006: we have µ n,k := E(X n,k ) = 2n, (n > k) k(k + 1) and 2(4k 2 3k 4)(k 1)n k(k + 1) 2, if n > 2k; (2k 1)(2k + 1) σn,k 2 := Var(X 2(5k 7)(k 1) n,k) = (k + 1) 2, if n = 2k; (2k 1) 2(k 2 + k 2n)n k 2 (k + 1) 2, if 2k > n > k. This result + central limit theorem also obtained by Devroye in 1991! Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
45 The phase change Question: What happens for k as n? Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
46 The phase change Question: What happens for k as n? Theorem (Feng, Mahmoud, Panholzer; 2008) (i) (Normal range) Let k = o ( n). Then, X n,k µ n,k σ n,k (ii) (Poisson range) Let k c n. Then, X n,k d N (0, 1). d Po(2c 2 ). (iii) (Degenerate range) Let k < n and n = o(k). Then, L X 1 n,k 0. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
47 A General Framework A pattern of size k is a set of induced subtrees of size k. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
48 A General Framework A pattern of size k is a set of induced subtrees of size k. X n,k = # of occurrence of a pattern of size k in random PT of size n. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
49 A General Framework A pattern of size k is a set of induced subtrees of size k. X n,k = # of occurrence of a pattern of size k in random PT of size n. We have, d X n,k = XIn,k + Xn I, n,k where X k,k = Bernoulli(p k ), X In,k and Xn I n,k are conditionally independent given I n, and I n = Unif{1,..., n 1} Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
50 A General Framework A pattern of size k is a set of induced subtrees of size k. X n,k = # of occurrence of a pattern of size k in random PT of size n. We have, d X n,k = XIn,k + Xn I, n,k where X k,k = Bernoulli(p k ), X In,k and Xn I n,k are conditionally independent given I n, and I n = Unif{1,..., n 1} Here, p k shape parameter 1 # of k-pronged nodes 2/(k 1) # of nodes with minimal clade size k 2 k 2 /(k 1)! # of k caterpillars Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
51 Mean Value and Variance We have and for n > 2k. µ n,k = 2p kn k(k + 1), (n > k), σ 2 n,k = 2(4k3 + 4k 2 k 1 (11k 2 5)p k )p k n k(k + 1) 2 (2k 1)(2k + 1) Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
52 Mean Value and Variance We have and for n > 2k. µ n,k = 2p kn k(k + 1), (n > k), σ 2 n,k = 2(4k3 + 4k 2 k 1 (11k 2 5)p k )p k n k(k + 1) 2 (2k 1)(2k + 1) Note that for n > 2k and k. µ n,k σ 2 n,k 2p kn k 2 Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
53 Poisson Approximation Theorem (Chang and F.; 2010) Let k < n and k. Then, d T V (X n,k, Po(µ n,k )) = { O (p k /k), if µ n,k 1; O (p k /k µ n,k ), if µ n,k < 1. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
54 Poisson Approximation Theorem (Chang and F.; 2010) Let k < n and k. Then, d T V (X n,k, Po(µ n,k )) = { O (p k /k), if µ n,k 1; O (p k /k µ n,k ), if µ n,k < 1. Recently re-proved by Holmgren and Janson with Stein s method: C. Holmgren and S. Janson (2015). Limit laws for functions of fringe trees for binary search trees and recursive trees, Electronic J. Probability, 20:4, Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
55 Berry-Esseen bound and LLT Theorem (Chang and F.; 2010) For µ n,k, ( sup P X n,k µ n,k σ n,k x R ) < x Φ(x) = O ( ) k. pk n Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
56 Berry-Esseen bound and LLT Theorem (Chang and F.; 2010) For µ n,k, ( sup P X n,k µ n,k σ n,k x R ) < x Φ(x) = O ( ) k. pk n Theorem (Chang and F.; 2010) For µ n,k, ( P (X n,k = µ n,k + xσn,k 2 ) = e x2 /2 (1 + O (1 + x 3 ) 2πσn,k k pk n )), uniformly in x = o((p k n) 1/6 /k 1/3 ). Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
57 Short Summary We have proved: Explicit expressions for mean and variance. Local limit theorem + Berry-Esseen bound for the largest possible range of k. Poisson approximation + rate whenever k. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
58 Short Summary We have proved: Explicit expressions for mean and variance. Local limit theorem + Berry-Esseen bound for the largest possible range of k. Poisson approximation + rate whenever k. So, for k n the number of occurrences of the pattern can be approximated by the normal distribution, whereas for the remaining range a Poisson random variable should be used. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
59 Short Summary We have proved: Explicit expressions for mean and variance. Local limit theorem + Berry-Esseen bound for the largest possible range of k. Poisson approximation + rate whenever k. So, for k n the number of occurrences of the pattern can be approximated by the normal distribution, whereas for the remaining range a Poisson random variable should be used. In particular, the phase change occurs much earlier than predicted by Feng, Mahmoud, and Panholzer. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
60 Lloyd Shapley Lloyd Shapley (1923-) Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
61 Lloyd Shapley Shapley value: Measure of importance of each player in a cooperative game Lloyd Shapley (1923-) Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
62 Lloyd Shapley Shapley value: Measure of importance of each player in a cooperative game recently used as prioritization tool of taxa in phylogenetics Lloyd Shapley (1923-) Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
63 Shapley Value and Modified Shapley Value T... PT; a... taxon (=leaf) of T. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
64 Shapley Value and Modified Shapley Value T... PT; a... taxon (=leaf) of T. Shapley value SV T (a): SV T (a) = 1 ( S 1)!(n S )!(PD T (S) PD T (S \ {a})). n! S,a S Modified Shapley value SV T (a): SV T (a) = 1 ( S 1)!(n S )!(PD T (S) PD T (S \ {a})). n! S 2,a S PD(S) is the size of the ancestor of S. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
65 Fair Proportion Index Fair proportion index FP T (a): FP T (a) = e 1 D e with D e the number of taxa below e. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
66 Fair Proportion Index Fair proportion index FP T (a): FP T (a) = e 1 D e with D e the number of taxa below e. Selected (somehow arbitrarily) by Zoological Society of London for EDGE of Existence conservation program! Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
67 Fair Proportion Index Fair proportion index FP T (a): FP T (a) = e 1 D e with D e the number of taxa below e. Selected (somehow arbitrarily) by Zoological Society of London for EDGE of Existence conservation program! FP n = fair proportion index of random taxon in random PT of size n: { 1 j FP n (I n = j) = + FP j, with probability j/n; 1 n j + FP n j, with probability (n j)/n, where I n = Unif{1,..., n 1}. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
68 Strong Correlation between SV and FP Hartmann (2013): Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
69 SV = FP Assume a is in left subtree T l and T l = j. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
70 SV = FP Assume a is in left subtree T l and T l = j. Lemma We have, and SV T (a) = 1 j + SV T l (a) FP T (a) = 1 j + FP T l (a). Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
71 SV = FP Assume a is in left subtree T l and T l = j. Lemma We have, and SV T (a) = 1 j + SV T l (a) FP T (a) = 1 j + FP T l (a). Theorem (F. and Jin; 2015) We have, SV T (a) = FP T (a). Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
72 SV and FP (i) D T (a) = depth of a in T : FP T (a) = SV T (a) = SV T (a) + D T (a) n. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
73 SV and FP (i) D T (a) = depth of a in T : FP T (a) = SV T (a) = SV T (a) + D T (a) n. D n = depth of random taxon in random PT of size n. Lemma We have, Var(FP n ) = 10 6H (2) n 1 6 n 4 10 π2 n2 Var(D n ) = 2H n 4H n (2) log n, where H n = 1 j n 1/j and H(2) n = 1 j n 1/j2. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
74 SV and FP (ii) Lemma We have, Cov(FP n, D n ) = 4H (2) n n + 4 n 2 2π Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
75 SV and FP (ii) Lemma We have, Cov(FP n, D n ) = 4H (2) n n + 4 n 2 2π Theorem (F. and Jin; 2015) The correlation coefficient ρ( SV n, FP n ) of modified Shapley value and fair proportion index tends to 1, i.e., lim ρ( SV n, FP n ) = 1. n Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
76 Short Summary Shapley value was introduced for the purpose of making conservation decisions in genetics. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
77 Short Summary Shapley value was introduced for the purpose of making conservation decisions in genetics. Fair proportion index was used by Zoological Society of London but its biodiversity value was unclear. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
78 Short Summary Shapley value was introduced for the purpose of making conservation decisions in genetics. Fair proportion index was used by Zoological Society of London but its biodiversity value was unclear. Strong correlation between modified Shapley value and fair proportion index was observed by Hartmann and others. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
79 Short Summary Shapley value was introduced for the purpose of making conservation decisions in genetics. Fair proportion index was used by Zoological Society of London but its biodiversity value was unclear. Strong correlation between modified Shapley value and fair proportion index was observed by Hartmann and others. We proved that Shapley value equals fair proportion index. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
80 Short Summary Shapley value was introduced for the purpose of making conservation decisions in genetics. Fair proportion index was used by Zoological Society of London but its biodiversity value was unclear. Strong correlation between modified Shapley value and fair proportion index was observed by Hartmann and others. We proved that Shapley value equals fair proportion index. We showed that correlation coefficient of modified Shapley value and fair proportion index tends to 1 in Yule-Harding model and other random models. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
81 Social Animals Consider n animals of a class of social animals. Construct a random PT with leaves representing the animals. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
82 Social Animals Consider n animals of a class of social animals. Construct a random PT with leaves representing the animals. Describes the genetic relatedness of the animals. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
83 Animal Groups Durand, Blum and François (2007): Groups are formed more likely by animals which are genetically related. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
84 Animal Groups Durand, Blum and François (2007): Groups are formed more likely by animals which are genetically related. neutral model. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
85 Animal Groups Durand, Blum and François (2007): Groups are formed more likely by animals which are genetically related. neutral model. Clade of a leaf: All leafs of the tree rooted at the parent. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
86 Animal Groups Durand, Blum and François (2007): Groups are formed more likely by animals which are genetically related. neutral model. Clade of a leaf: All leafs of the tree rooted at the parent. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
87 Animal Groups Durand, Blum and François (2007): Groups are formed more likely by animals which are genetically related. neutral model. # of groups # of maximal clades Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
88 Animal Groups Durand, Blum and François (2007): Groups are formed more likely by animals which are genetically related. neutral model. # of groups # of maximal clades 2 Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
89 # of Groups X n = # of groups under the Yule Harding model Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
90 # of Groups X n = # of groups under the Yule Harding model We have, { d 1, if I n = 1 or I n = n 1, X n = X In + Xn I n, otherwise, where I n = Uniform{1,..., n 1} is the # of animals in the left subtree and X n is an independent copy of X n. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
91 # of Groups X n = # of groups under the Yule Harding model We have, { d 1, if I n = 1 or I n = n 1, X n = X In + Xn I n, otherwise, where I n = Uniform{1,..., n 1} is the # of animals in the left subtree and X n is an independent copy of X n. Theorem (Durand and François; 2010) We have, E(X n ) an (a := 1 e 2 4 ). Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
92 Comparison with Real-life Data Durand, Blum and François (2007) presented the following data: Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
93 Yi-Wen s Thesis (2012) Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
94 Variance and SLLN Theorem (Lee; 2012) We have, Var(X n ) (1 e 2 ) 2 n log n = 4a 2 n log n. 4 Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
95 Variance and SLLN Theorem (Lee; 2012) We have, Var(X n ) (1 e 2 ) 2 n log n = 4a 2 n log n. 4 Theorem (Lee; 2012) We have, ( ) P lim X n n E(X n ) 1 = 0 = 1. For SLLN, X n is constructed on the same probability space via the tree evolution process underlying the Yule-Harding model. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
96 Method of Moments Theorem Assume that E(X k n) E(X k ) for all k 1 and that X is uniquely characterised by its sequence of moments. Then, X n d X. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
97 Method of Moments Theorem Assume that E(X k n) E(X k ) for all k 1 and that X is uniquely characterised by its sequence of moments. Then, X n d X. Many sufficient conditions for X being uniquely characterized by its moments are known, e.g., E(X k ) zk k! k 1 has a positive radius of convergence. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
98 Higher Moments Theorem (Lee; 2012) For all k 3, E(X n E(X n )) k ( 1) k 2k k 2 ak n k 1. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
99 Higher Moments Theorem (Lee; 2012) For all k 3, E(X n E(X n )) k ( 1) k 2k k 2 ak n k 1. This implies that all moments larger than two of X n E(X n ) Var(Xn ) tend to infinity! Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
100 Higher Moments Theorem (Lee; 2012) For all k 3, E(X n E(X n )) k ( 1) k 2k k 2 ak n k 1. This implies that all moments larger than two of X n E(X n ) Var(Xn ) tend to infinity! Question: Is there a limit distribution? Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
101 Random Recursive Trees Unordered, rooted trees. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
102 Random Recursive Trees Unordered, rooted trees. Uniformly choose one of the nodes and attach a child. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
103 Random Recursive Trees Unordered, rooted trees. Uniformly choose one of the nodes and attach a child. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
104 Random Recursive Trees Unordered, rooted trees. Uniformly choose one of the nodes and attach a child. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
105 Random Recursive Trees Unordered, rooted trees. Uniformly choose one of the nodes and attach a child. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
106 Random Recursive Trees Unordered, rooted trees. Uniformly choose one of the nodes and attach a child. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
107 Random Recursive Trees Unordered, rooted trees. Uniformly choose one of the nodes and attach a child. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
108 Random Recursive Trees Unordered, rooted trees. Uniformly choose one of the nodes and attach a child. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
109 Cutting Down Random Recursive Trees Meir and Moon (1974): Randomly pick an edge and remove it; retain the tree containing the root. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
110 Cutting Down Random Recursive Trees Meir and Moon (1974): Randomly pick an edge and remove it; retain the tree containing the root. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
111 Cutting Down Random Recursive Trees Meir and Moon (1974): Randomly pick an edge and remove it; retain the tree containing the root. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
112 Cutting Down Random Recursive Trees Meir and Moon (1974): Randomly pick an edge and remove it; retain the tree containing the root. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
113 Cutting Down Random Recursive Trees Meir and Moon (1974): Randomly pick an edge and remove it; retain the tree containing the root. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
114 Cutting Down Random Recursive Trees Meir and Moon (1974): Randomly pick an edge and remove it; retain the tree containing the root. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
115 Cutting Down Random Recursive Trees Meir and Moon (1974): Randomly pick an edge and remove it; retain the tree containing the root. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
116 Cutting Down Random Recursive Trees Meir and Moon (1974): Randomly pick an edge and remove it; retain the tree containing the root. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
117 Cutting Down Random Recursive Trees Meir and Moon (1974): Randomly pick an edge and remove it; retain the tree containing the root. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
118 Cutting Down Random Recursive Trees Meir and Moon (1974): Randomly pick an edge and remove it; retain the tree containing the root. Y n = number of steps until tree is destroyed = number of edges cut = 4. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
119 Mean, Variance and Higher Moments Theorem (Panholzer; 2004) We have, and for k 2 E(Y n ) n log n E(Y n E(Y n )) k ( 1)k k(k 1) n k log k+1 n. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
120 Mean, Variance and Higher Moments Theorem (Panholzer; 2004) We have, and for k 2 E(Y n ) n log n E(Y n E(Y n )) k Thus, again the limit law of ( 1)k k(k 1) n k log k+1 n. Y n E(Y n ) Var(Yn ) cannot obtained from the method of moments! Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
121 Limit Law Theorem (Drmota, Iksanov, Moehle, Roessler; 2009) We have, with log 2 n n Y n log n log log n d Y E(e iλy ) = e iλ log λ π λ /2. The law of Y is spectrally negative stable with index of stability 1. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
122 Limit Law Theorem (Drmota, Iksanov, Moehle, Roessler; 2009) We have, with log 2 n n Y n log n log log n d Y E(e iλy ) = e iλ log λ π λ /2. The law of Y is spectrally negative stable with index of stability 1. Different proofs of this result exist. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
123 Limit Law of X n Theorem (Drmota, F., Lee; 2014) We have, X n E(X n ) Var(Xn )/2 d N(0, 1). Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
124 Limit Law of X n Theorem (Drmota, F., Lee; 2014) We have, X n E(X n ) Var(Xn )/2 d N(0, 1). For the proof, we use singularity perturbation analysis. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
125 Limit Law of X n Theorem (Drmota, F., Lee; 2014) We have, X n E(X n ) Var(Xn )/2 d N(0, 1). For the proof, we use singularity perturbation analysis. A probabilistic proof explaining the curious normalization was given recently. S. Janson (2015). Maximal clades in random binary search trees, Electronic J. Combinatorics, 22:1, paper 31. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
126 Some Ideas of the Proof (i) Set X(y, z) = n 2 E ( e yxn) z n. Then, This is a Riccati DE. z z X(y, z) = X(y, z) + X2 (y, z) + e y z 2 2ey z 3 1 z. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
127 Some Ideas of the Proof (i) Set X(y, z) = n 2 E ( e yxn) z n. Then, This is a Riccati DE. z z X(y, z) = X(y, z) + X2 (y, z) + e y z 2 2ey z 3 1 z. Set Then, X(y, z) = X(y, z). z z X(y, z) = X 2 (y, z) + e y 1 + z 1 z. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
128 Some Ideas of the Proof (ii) Set Then, This is Whittaker s DE. X(y, z) = V (y, z) V (y, z). V (y, z) + e y 1 + z V (y, z) = 0. 1 z Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
129 Some Ideas of the Proof (ii) Set Then, This is Whittaker s DE. Solution is given by V (y, z) = M e y/2,1/2 X(y, z) = V (y, z) V (y, z). V (y, z) + e y 1 + z V (y, z) = 0. 1 z ( ) ( ) 2e y/2 (z 1) + c(y)w e y/2,1/2 2e y/2 (z 1), where c(y) = ( e y/2 1 ) M e y/2 +1,1/2 ( 2e y/2 ) ( W ) e y/2 +1,1/2 2e y/2. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
130 Some Ideas of the Proof (iii) Lemma V (y, z) is analytic in = {z C : z < 1 + δ}\ {branch cut from 1 to } for all y < η. Moreover, V (y, z) has only one (simple) zero with z 0 (y) = 1 ay + 2a 2 y 2 log y + O(y 2 ). Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
131 Some Ideas of the Proof (iii) Lemma V (y, z) is analytic in = {z C : z < 1 + δ}\ {branch cut from 1 to } for all y < η. Moreover, V (y, z) has only one (simple) zero with z 0 (y) = 1 ay + 2a 2 y 2 log y + O(y 2 ). δ z 0 (y) 1 Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
132 Some Ideas of the Proof (iv) Let y = it/(2a n log n). Then, E ( e yxn) = 1 2πi C X(y, z) z n+1 dz. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
133 Some Ideas of the Proof (iv) Let y = it/(2a n log n). Then, E ( e yxn) = 1 2πi C X(y, z) z n+1 dz. Lemma We have, E ( e yxn) ( log = z 0 (y) n 3 ) n + O. n Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
134 Some Ideas of the Proof (iv) Let y = it/(2a n log n). Then, E ( e yxn) = 1 2πi C X(y, z) z n+1 dz. Lemma We have, E ( e yxn) ( log = z 0 (y) n 3 ) n + O. n This together with the expansion of z 0 (y) yields E ( e yxn) ( ) ( it n = exp 2 log n t2 1 + O 4 ( log log n log n )). Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
135 Summary Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
136 Summary We gave a comprehensive study of pattern occurrences in random phylogenetic trees. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
137 Summary We gave a comprehensive study of pattern occurrences in random phylogenetic trees. We explained the strong correlation between Shapley value and fair proportion index in phylogenetic trees. This is of great interest for people working in biodiversity. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
138 Summary We gave a comprehensive study of pattern occurrences in random phylogenetic trees. We explained the strong correlation between Shapley value and fair proportion index in phylogenetic trees. This is of great interest for people working in biodiversity. We found a curious central limit theorem for the number of groups formed by social animals under the neutral model. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
139 Summary We gave a comprehensive study of pattern occurrences in random phylogenetic trees. We explained the strong correlation between Shapley value and fair proportion index in phylogenetic trees. This is of great interest for people working in biodiversity. We found a curious central limit theorem for the number of groups formed by social animals under the neutral model. Analytic combinatorics is useful in studying mathematical problems for random phylogenetic trees. We expect many more applications. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
140 Summary We gave a comprehensive study of pattern occurrences in random phylogenetic trees. We explained the strong correlation between Shapley value and fair proportion index in phylogenetic trees. This is of great interest for people working in biodiversity. We found a curious central limit theorem for the number of groups formed by social animals under the neutral model. Analytic combinatorics is useful in studying mathematical problems for random phylogenetic trees. We expect many more applications. More input from biologists is needed to make our results more relevant from the point of view of applications. Michael Fuchs (NCTU) Phylogenetic Trees (PTs) June 30th, / 45
On the number of subtrees on the fringe of random trees (partly joint with Huilan Chang)
On the number of subtrees on the fringe of random trees (partly joint with Huilan Chang) Michael Fuchs Institute of Applied Mathematics National Chiao Tung University Hsinchu, Taiwan April 16th, 2008 Michael
More informationThe Subtree Size Profile of Plane-oriented Recursive Trees
The Subtree Size Profile of Plane-oriented Recursive Trees Michael Fuchs Department of Applied Mathematics National Chiao Tung University Hsinchu, Taiwan ANALCO11, January 22nd, 2011 Michael Fuchs (NCTU)
More informationMAXIMAL CLADES IN RANDOM BINARY SEARCH TREES
MAXIMAL CLADES IN RANDOM BINARY SEARCH TREES SVANTE JANSON Abstract. We study maximal clades in random phylogenetic trees with the Yule Harding model or, equivalently, in binary search trees. We use probabilistic
More informationThe Subtree Size Profile of Bucket Recursive Trees
Iranian Journal of Mathematical Sciences and Informatics Vol., No. (206, pp - DOI: 0.7508/ijmsi.206.0.00 The Subtree Size Profile of Bucket Recursive Trees Ramin Kazemi Department of Statistics, Imam Khomeini
More informationAlmost giant clusters for percolation on large trees
for percolation on large trees Institut für Mathematik Universität Zürich Erdős-Rényi random graph model in supercritical regime G n = complete graph with n vertices Bond percolation with parameter p(n)
More informationarxiv: v1 [q-bio.pe] 13 Aug 2015
Noname manuscript No. (will be inserted by the editor) On joint subtree distributions under two evolutionary models Taoyang Wu Kwok Pui Choi arxiv:1508.03139v1 [q-bio.pe] 13 Aug 2015 August 14, 2015 Abstract
More informationOn joint subtree distributions under two evolutionary models
Noname manuscript No. (will be inserted by the editor) On joint subtree distributions under two evolutionary models Taoyang Wu Kwok Pui Choi November 5, 2015 Abstract In population and evolutionary biology,
More informationRandom trees and branching processes
Random trees and branching processes Svante Janson IMS Medallion Lecture 12 th Vilnius Conference and 2018 IMS Annual Meeting Vilnius, 5 July, 2018 Part I. Galton Watson trees Let ξ be a random variable
More informationThe Subtree Size Profile of Plane-oriented Recursive Trees
The Subtree Size Profile of Plane-oriente Recursive Trees Michael FUCHS Department of Applie Mathematics National Chiao Tung University Hsinchu, 3, Taiwan Email: mfuchs@math.nctu.eu.tw Abstract In this
More informationDynamics of the evolving Bolthausen-Sznitman coalescent. by Jason Schweinsberg University of California at San Diego.
Dynamics of the evolving Bolthausen-Sznitman coalescent by Jason Schweinsberg University of California at San Diego Outline of Talk 1. The Moran model and Kingman s coalescent 2. The evolving Kingman s
More informationPROTECTED NODES AND FRINGE SUBTREES IN SOME RANDOM TREES
PROTECTED NODES AND FRINGE SUBTREES IN SOME RANDOM TREES LUC DEVROYE AND SVANTE JANSON Abstract. We study protected nodes in various classes of random rooted trees by putting them in the general context
More informationDependence between Path Lengths and Size in Random Trees (joint with H.-H. Chern, H.-K. Hwang and R. Neininger)
Dependence between Path Lengths and Size in Random Trees (joint with H.-H. Chern, H.-K. Hwang and R. Neininger) Michael Fuchs Institute of Applied Mathematics National Chiao Tung University Hsinchu, Taiwan
More informationCutting edges at random in large recursive trees
Cutting edges at random in large recursive trees arxiv:1406.2238v1 [math.pr] 9 Jun 2014 Erich Baur and Jean Bertoin ENS Lyon and Universität Zürich February 21, 2018 Abstract We comment on old and new
More informationk-protected VERTICES IN BINARY SEARCH TREES
k-protected VERTICES IN BINARY SEARCH TREES MIKLÓS BÓNA Abstract. We show that for every k, the probability that a randomly selected vertex of a random binary search tree on n nodes is at distance k from
More informationDISTRIBUTIONS OF CHERRIES FOR TWO MODELS OF TREES
DISTRIBUTIONS OF CHERRIES FOR TWO MODELS OF TREES ANDY McKENZIE and MIKE STEEL Biomathematics Research Centre University of Canterbury Private Bag 4800 Christchurch, New Zealand No. 177 May, 1999 Distributions
More informationProtected nodes and fringe subtrees in some random trees
Electron. Commun. Probab. 19 (214), no. 6, 1 1. DOI: 1.1214/ECP.v19-348 ISSN: 183-589X ELECTRONIC COMMUNICATIONS in PROBABILITY Protected nodes and fringe subtrees in some random trees Luc Devroye Svante
More informationON COMPOUND POISSON POPULATION MODELS
ON COMPOUND POISSON POPULATION MODELS Martin Möhle, University of Tübingen (joint work with Thierry Huillet, Université de Cergy-Pontoise) Workshop on Probability, Population Genetics and Evolution Centre
More informationThe expected value of the squared euclidean cophenetic metric under the Yule and the uniform models
The expected value of the squared euclidean cophenetic metric under the Yule and the uniform models Gabriel Cardona, Arnau Mir, Francesc Rosselló Department of Mathematics and Computer Science, University
More informationThe range of tree-indexed random walk
The range of tree-indexed random walk Jean-François Le Gall, Shen Lin Institut universitaire de France et Université Paris-Sud Orsay Erdös Centennial Conference July 2013 Jean-François Le Gall (Université
More informationRECOVERING NORMAL NETWORKS FROM SHORTEST INTER-TAXA DISTANCE INFORMATION
RECOVERING NORMAL NETWORKS FROM SHORTEST INTER-TAXA DISTANCE INFORMATION MAGNUS BORDEWICH, KATHARINA T. HUBER, VINCENT MOULTON, AND CHARLES SEMPLE Abstract. Phylogenetic networks are a type of leaf-labelled,
More informationChapter 11. Min Cut Min Cut Problem Definition Some Definitions. By Sariel Har-Peled, December 10, Version: 1.
Chapter 11 Min Cut By Sariel Har-Peled, December 10, 013 1 Version: 1.0 I built on the sand And it tumbled down, I built on a rock And it tumbled down. Now when I build, I shall begin With the smoke from
More informationAsymptotic distribution of two-protected nodes in ternary search trees
Asymptotic distribution of two-protected nodes in ternary search trees Cecilia Holmgren Svante Janson March 2, 204; revised October 5, 204 Abstract We study protected nodes in m-ary search trees, by putting
More informationarxiv: v1 [math.pr] 21 Mar 2014
Asymptotic distribution of two-protected nodes in ternary search trees Cecilia Holmgren Svante Janson March 2, 24 arxiv:4.557v [math.pr] 2 Mar 24 Abstract We study protected nodes in m-ary search trees,
More informationOn Symmetries of Non-Plane Trees in a Non-Uniform Model
On Symmetries of Non-Plane Trees in a Non-Uniform Model Jacek Cichoń Abram Magner Wojciech Spankowski Krystof Turowski October 3, 26 Abstract Binary trees come in two varieties: plane trees, often simply
More informationCENTRAL LIMIT THEOREMS FOR ADDITIVE TREE PARAMETERS
CENRAL LIMI HEOREMS FOR ADDIIVE REE PARAMEERS SEPHAN WAGNER DEDICAED O HE MEMORY OF PHILIPPE FLAJOLE Abstract. We call a tree parameter additive if it can be determined recursively as the sum of the parameter
More informationAsymptotic and Exact Poissonized Variance in the Analysis of Random Digital Trees (joint with Hsien-Kuei Hwang and Vytas Zacharovas)
Asymptotic and Exact Poissonized Variance in the Analysis of Random Digital Trees (joint with Hsien-Kuei Hwang and Vytas Zacharovas) Michael Fuchs Institute of Applied Mathematics National Chiao Tung University
More informationON THE DISTRIBUTION OF SUBTREE ORDERS OF A TREE. 1. Introduction
ON THE DISTRIBUTION OF SUBTREE ORDERS OF A TREE DIMBINAINA RALAIVAOSAONA AND STEPHAN WAGNER Abstract. We investigate the distribution of the number of vertices of a randomly chosen subtree of a tree. Specifically,
More informationANALYSIS OF SOME TREE STATISTICS
ANALYSIS OF SOME TREE STATISTICS ALOIS PANHOLZER Abstract. In this survey we present results and the used proof techniques concerning the analysis of several parameters in frequently used tree families.
More informationWorkshop III: Evolutionary Genomics
Identifying Species Trees from Gene Trees Elizabeth S. Allman University of Alaska IPAM Los Angeles, CA November 17, 2011 Workshop III: Evolutionary Genomics Collaborators The work in today s talk is joint
More informationModern Discrete Probability Branching processes
Modern Discrete Probability IV - Branching processes Review Sébastien Roch UW Madison Mathematics November 15, 2014 1 Basic definitions 2 3 4 Galton-Watson branching processes I Definition A Galton-Watson
More informationarxiv: v1 [math.co] 12 Apr 2018
COMBINATORIAL ANALYSIS OF GROWTH MODELS FOR SERIES-PARALLEL NETWORKS MARKUS KUBA AND ALOIS PANHOLZER arxiv:804.0550v [math.co] Apr 08 Abstract. We give combinatorial descriptions of two stochastic growth
More informationarxiv: v1 [math.co] 1 Aug 2018
REDUCING SIMPLY GENERATED TREES BY ITERATIVE LEAF CUTTING BENJAMIN HACKL, CLEMENS HEUBERGER, AND STEPHAN WAGNER arxiv:1808.00363v1 [math.co] 1 Aug 2018 ABSTRACT. We consider a procedure to reduce simply
More informationMod-φ convergence I: examples and probabilistic estimates
Mod-φ convergence I: examples and probabilistic estimates Valentin Féray (joint work with Pierre-Loïc Méliot and Ashkan Nikeghbali) Institut für Mathematik, Universität Zürich Summer school in Villa Volpi,
More informationAlmost sure asymptotics for the random binary search tree
AofA 10 DMTCS proc. AM, 2010, 565 576 Almost sure asymptotics for the rom binary search tree Matthew I. Roberts Laboratoire de Probabilités et Modèles Aléatoires, Université Paris VI Case courrier 188,
More information14 Branching processes
4 BRANCHING PROCESSES 6 4 Branching processes In this chapter we will consider a rom model for population growth in the absence of spatial or any other resource constraints. So, consider a population of
More informationApproximate Counting via the Poisson-Laplace-Mellin Method (joint with Chung-Kuei Lee and Helmut Prodinger)
Approximate Counting via the Poisson-Laplace-Mellin Method (joint with Chung-Kuei Lee and Helmut Prodinger) Michael Fuchs Department of Applied Mathematics National Chiao Tung University Hsinchu, Taiwan
More informationPaul-Eugène Parent. March 12th, Department of Mathematics and Statistics University of Ottawa. MAT 3121: Complex Analysis I
Paul-Eugène Parent Department of Mathematics and Statistics University of Ottawa March 12th, 2014 Outline 1 Holomorphic power Series Proposition Let f (z) = a n (z z o ) n be the holomorphic function defined
More informationForty years of tree enumeration
October 24, 206 Canada John Moon (Alberta) Counting labelled trees Canadian mathematical monographs What is our subject? Asymptotic enumeration or Probabilistic combinatorics? What is our subject? Asymptotic
More informationLecture 2: Discrete Probability Distributions
Lecture 2: Discrete Probability Distributions IB Paper 7: Probability and Statistics Carl Edward Rasmussen Department of Engineering, University of Cambridge February 1st, 2011 Rasmussen (CUED) Lecture
More informationROSENA R.X. DU AND HELMUT PRODINGER. Dedicated to Philippe Flajolet ( )
ON PROTECTED NODES IN DIGITAL SEARCH TREES ROSENA R.X. DU AND HELMUT PRODINGER Dedicated to Philippe Flajolet (98 2) Abstract. Recently, 2-protected nodes were studied in the context of ordered trees and
More informationHard-Core Model on Random Graphs
Hard-Core Model on Random Graphs Antar Bandyopadhyay Theoretical Statistics and Mathematics Unit Seminar Theoretical Statistics and Mathematics Unit Indian Statistical Institute, New Delhi Centre New Delhi,
More informationFRINGE TREES, CRUMP MODE JAGERS BRANCHING PROCESSES AND m-ary SEARCH TREES
FRINGE TREES, CRUMP MODE JAGERS BRANCHING PROCESSES AND m-ary SEARCH TREES CECILIA HOLMGREN AND SVANTE JANSON Abstract. This survey studies asymptotics of random fringe trees and extended fringe trees
More informationThe Moran Process as a Markov Chain on Leaf-labeled Trees
The Moran Process as a Markov Chain on Leaf-labeled Trees David J. Aldous University of California Department of Statistics 367 Evans Hall # 3860 Berkeley CA 94720-3860 aldous@stat.berkeley.edu http://www.stat.berkeley.edu/users/aldous
More informationarxiv: v1 [q-bio.pe] 1 Jun 2014
THE MOST PARSIMONIOUS TREE FOR RANDOM DATA MAREIKE FISCHER, MICHELLE GALLA, LINA HERBST AND MIKE STEEL arxiv:46.27v [q-bio.pe] Jun 24 Abstract. Applying a method to reconstruct a phylogenetic tree from
More informationNotes by Zvi Rosen. Thanks to Alyssa Palfreyman for supplements.
Lecture: Hélène Barcelo Analytic Combinatorics ECCO 202, Bogotá Notes by Zvi Rosen. Thanks to Alyssa Palfreyman for supplements.. Tuesday, June 2, 202 Combinatorics is the study of finite structures that
More informationMidterm Exam 1 Solution
EECS 126 Probability and Random Processes University of California, Berkeley: Fall 2015 Kannan Ramchandran September 22, 2015 Midterm Exam 1 Solution Last name First name SID Name of student on your left:
More informationn! (k 1)!(n k)! = F (X) U(0, 1). (x, y) = n(n 1) ( F (y) F (x) ) n 2
Order statistics Ex. 4. (*. Let independent variables X,..., X n have U(0, distribution. Show that for every x (0,, we have P ( X ( < x and P ( X (n > x as n. Ex. 4.2 (**. By using induction or otherwise,
More informationThe expansion of random regular graphs
The expansion of random regular graphs David Ellis Introduction Our aim is now to show that for any d 3, almost all d-regular graphs on {1, 2,..., n} have edge-expansion ratio at least c d d (if nd is
More informationSusceptible-Infective-Removed Epidemics and Erdős-Rényi random
Susceptible-Infective-Removed Epidemics and Erdős-Rényi random graphs MSR-Inria Joint Centre October 13, 2015 SIR epidemics: the Reed-Frost model Individuals i [n] when infected, attempt to infect all
More informationOn the number of matchings of a tree.
On the number of matchings of a tree. Stephan G. Wagner Department of Mathematics, Graz University of Technology, Steyrergasse 30, A-800 Graz, Austria Abstract In a paper of Klazar, several counting examples
More informationNetworks: Lectures 9 & 10 Random graphs
Networks: Lectures 9 & 10 Random graphs Heather A Harrington Mathematical Institute University of Oxford HT 2017 What you re in for Week 1: Introduction and basic concepts Week 2: Small worlds Week 3:
More information5. Analytic Combinatorics
ANALYTIC COMBINATORICS P A R T O N E 5. Analytic Combinatorics http://aofa.cs.princeton.edu Analytic combinatorics is a calculus for the quantitative study of large combinatorial structures. Features:
More informationProbabilistic analysis of the asymmetric digital search trees
Int. J. Nonlinear Anal. Appl. 6 2015 No. 2, 161-173 ISSN: 2008-6822 electronic http://dx.doi.org/10.22075/ijnaa.2015.266 Probabilistic analysis of the asymmetric digital search trees R. Kazemi a,, M. Q.
More informationRandomized Algorithms III Min Cut
Chapter 11 Randomized Algorithms III Min Cut CS 57: Algorithms, Fall 01 October 1, 01 11.1 Min Cut 11.1.1 Problem Definition 11. Min cut 11..0.1 Min cut G = V, E): undirected graph, n vertices, m edges.
More informationOn the distribution of subtree orders of a tree
Also available at http://amc-journal.eu ISSN 855-3966 printed edn., ISSN 855-3974 electronic edn. ARS MATHEMATICA CONTEMPORANEA 4 208 29 56 On the distribution of subtree orders of a tree Abstract Dimbinaina
More informationGibbs distributions for random partitions generated by a fragmentation process
Gibbs distributions for random partitions generated by a fragmentation process Nathanaël Berestycki and Jim Pitman University of British Columbia, and Department of Statistics, U.C. Berkeley July 24, 2006
More informationMath Bootcamp 2012 Miscellaneous
Math Bootcamp 202 Miscellaneous Factorial, combination and permutation The factorial of a positive integer n denoted by n!, is the product of all positive integers less than or equal to n. Define 0! =.
More informationUniversity of Regina. Lecture Notes. Michael Kozdron
University of Regina Statistics 252 Mathematical Statistics Lecture Notes Winter 2005 Michael Kozdron kozdron@math.uregina.ca www.math.uregina.ca/ kozdron Contents 1 The Basic Idea of Statistics: Estimating
More informationarxiv: v1 [cs.cc] 9 Oct 2014
Satisfying ternary permutation constraints by multiple linear orders or phylogenetic trees Leo van Iersel, Steven Kelk, Nela Lekić, Simone Linz May 7, 08 arxiv:40.7v [cs.cc] 9 Oct 04 Abstract A ternary
More informationRANDOM CUTTING AND RECORDS IN DETERMINISTIC AND RANDOM TREES
RANDOM CUTTING AND RECORDS IN DETERMINISTIC AND RANDOM TREES SVANTE JANSON Abstract. We study random cutting down of a rooted tree and show that the number of cuts is equal (in distribution) to the number
More informationarxiv:math.pr/ v1 17 May 2004
Probabilistic Analysis for Randomized Game Tree Evaluation Tämur Ali Khan and Ralph Neininger arxiv:math.pr/0405322 v1 17 May 2004 ABSTRACT: We give a probabilistic analysis for the randomized game tree
More informationFrom Individual-based Population Models to Lineage-based Models of Phylogenies
From Individual-based Population Models to Lineage-based Models of Phylogenies Amaury Lambert (joint works with G. Achaz, H.K. Alexander, R.S. Etienne, N. Lartillot, H. Morlon, T.L. Parsons, T. Stadler)
More information18.175: Lecture 17 Poisson random variables
18.175: Lecture 17 Poisson random variables Scott Sheffield MIT 1 Outline More on random walks and local CLT Poisson random variable convergence Extend CLT idea to stable random variables 2 Outline More
More informationMaximum Agreement Subtrees
Maximum Agreement Subtrees Seth Sullivant North Carolina State University March 24, 2018 Seth Sullivant (NCSU) Maximum Agreement Subtrees March 24, 2018 1 / 23 Phylogenetics Problem Given a collection
More informationOn finding a minimum spanning tree in a network with random weights
On finding a minimum spanning tree in a network with random weights Colin McDiarmid Department of Statistics University of Oxford Harold S. Stone IBM T.J. Watson Center 30 August 1996 Theodore Johnson
More informationFINAL EXAM: Monday 8-10am
ECE 30: Probabilistic Methods in Electrical and Computer Engineering Fall 016 Instructor: Prof. A. R. Reibman FINAL EXAM: Monday 8-10am Fall 016, TTh 3-4:15pm (December 1, 016) This is a closed book exam.
More informationGENERATING FUNCTIONS AND CENTRAL LIMIT THEOREMS
GENERATING FUNCTIONS AND CENTRAL LIMIT THEOREMS Michael Drmota Institute of Discrete Mathematics and Geometry Vienna University of Technology A 1040 Wien, Austria michael.drmota@tuwien.ac.at http://www.dmg.tuwien.ac.at/drmota/
More informationThe Wright-Fisher Model and Genetic Drift
The Wright-Fisher Model and Genetic Drift January 22, 2015 1 1 Hardy-Weinberg Equilibrium Our goal is to understand the dynamics of allele and genotype frequencies in an infinite, randomlymating population
More informationDiscrete Distributions Chapter 6
Discrete Distributions Chapter 6 Negative Binomial Distribution section 6.3 Consider k r, r +,... independent Bernoulli trials with probability of success in one trial being p. Let the random variable
More informationCombinatorial/probabilistic analysis of a class of search-tree functionals
Combinatorial/probabilistic analysis of a class of search-tree functionals Jim Fill jimfill@jhu.edu http://www.mts.jhu.edu/ fill/ Mathematical Sciences The Johns Hopkins University Joint Statistical Meetings;
More informationarxiv: v1 [math.co] 29 Nov 2018
ON THE INDUCIBILITY OF SMALL TREES AUDACE A. V. DOSSOU-OLORY AND STEPHAN WAGNER arxiv:1811.1010v1 [math.co] 9 Nov 018 Abstract. The quantity that captures the asymptotic value of the maximum number of
More informationHow should we organize the diversity of animal life?
How should we organize the diversity of animal life? The difference between Taxonomy Linneaus, and Cladistics Darwin What are phylogenies? How do we read them? How do we estimate them? Classification (Taxonomy)
More informationAdditive tree functionals with small toll functions and subtrees of random trees
Additive tree functionals with small toll functions and subtrees of random trees Stephan Wagner To cite this version: Stephan Wagner. Additive tree functionals with small toll functions and subtrees of
More informationGrowing and Destroying Catalan Stanley Trees
Discrete Mathematics and Theoretical Computer Science DMTCS vol. 20:1, 2018, #11 Growing and Destroying Catalan Stanley Trees Benjamin Hackl 1 Helmut Prodinger 2 arxiv:1704.03734v3 [math.co] 26 Feb 2018
More informationThis exam is closed book and closed notes. (You will have access to a copy of the Table of Common Distributions given in the back of the text.
TEST #3 STA 536 December, 00 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. You will have access to a copy
More informationDiscrete Random Variable
Discrete Random Variable Outcome of a random experiment need not to be a number. We are generally interested in some measurement or numerical attribute of the outcome, rather than the outcome itself. n
More informationPermutations with Inversions
1 2 3 47 6 23 11 Journal of Integer Sequences, Vol. 4 2001, Article 01.2.4 Permutations with Inversions Barbara H. Margolius Cleveland State University Cleveland, Ohio 44115 Email address: b.margolius@csuohio.edu
More informationTrees, Tanglegrams, and Tangled Chains
Trees, Tanglegrams, and Tangled Chains Sara Billey University of Washington Based on joint work with: Matjaž Konvalinka and Frederick (Erick) Matsen IV arxiv:1507.04976 LSU Math Student Colloquium, February
More informationA note on a scaling limit of successive approximation for a differential equation with moving singularity.
A note on a scaling limit of successive approximation for a differential equation with moving singularity. Tetsuya Hattori and Hiroyuki Ochiai ABSTRACT We study a method of successive approximation to
More informationThe space requirement of m-ary search trees: distributional asymptotics for m 27
1 The space requirement of m-ary search trees: distributional asymptotics for m 7 James Allen Fill 1 and Nevin Kapur 1 Applied Mathematics and Statistics, The Johns Hopkins University, 3400 N. Charles
More informationTheory of Stochastic Processes 3. Generating functions and their applications
Theory of Stochastic Processes 3. Generating functions and their applications Tomonari Sei sei@mist.i.u-tokyo.ac.jp Department of Mathematical Informatics, University of Tokyo April 20, 2017 http://www.stat.t.u-tokyo.ac.jp/~sei/lec.html
More informationRandom Boolean expressions
Random Boolean expressions Danièle Gardy PRiSM, Univ. Versailles-Saint Quentin and CNRS UMR 8144. Daniele.Gardy@prism.uvsq.fr We examine how we can define several probability distributions on the set of
More informationAverage Case Analysis of QuickSort and Insertion Tree Height using Incompressibility
Average Case Analysis of QuickSort and Insertion Tree Height using Incompressibility Tao Jiang, Ming Li, Brendan Lucier September 26, 2005 Abstract In this paper we study the Kolmogorov Complexity of a
More informationn! (k 1)!(n k)! = F (X) U(0, 1). (x, y) = n(n 1) ( F (y) F (x) ) n 2
Order statistics Ex. 4.1 (*. Let independent variables X 1,..., X n have U(0, 1 distribution. Show that for every x (0, 1, we have P ( X (1 < x 1 and P ( X (n > x 1 as n. Ex. 4.2 (**. By using induction
More informationRandomized algorithm
Tutorial 4 Joyce 2009-11-24 Outline Solution to Midterm Question 1 Question 2 Question 1 Question 2 Question 3 Question 4 Question 5 Solution to Midterm Solution to Midterm Solution to Midterm Question
More informationAnalysis of Approximate Quickselect and Related Problems
Analysis of Approximate Quickselect and Related Problems Conrado Martínez Univ. Politècnica Catalunya Joint work with A. Panholzer and H. Prodinger LIP6, Paris, April 2009 Introduction Quickselect finds
More informationPart IA Probability. Definitions. Based on lectures by R. Weber Notes taken by Dexter Chua. Lent 2015
Part IA Probability Definitions Based on lectures by R. Weber Notes taken by Dexter Chua Lent 2015 These notes are not endorsed by the lecturers, and I have modified them (often significantly) after lectures.
More informationExtremal Statistics on Non-Crossing Configurations
Extremal Statistics on Non-Crossing Configurations Michael Drmota Anna de Mier Marc Noy Abstract We analye extremal statistics in non-crossing configurations on the n vertices of a convex polygon. We prove
More information4 Sums of Independent Random Variables
4 Sums of Independent Random Variables Standing Assumptions: Assume throughout this section that (,F,P) is a fixed probability space and that X 1, X 2, X 3,... are independent real-valued random variables
More informationFE 5204 Stochastic Differential Equations
Instructor: Jim Zhu e-mail:zhu@wmich.edu http://homepages.wmich.edu/ zhu/ January 20, 2009 Preliminaries for dealing with continuous random processes. Brownian motions. Our main reference for this lecture
More informationExpectation MATH Expectation. Benjamin V.C. Collins, James A. Swenson MATH 2730
MATH 2730 Expectation Benjamin V.C. Collins James A. Swenson Average value Expectation Definition If (S, P) is a sample space, then any function with domain S is called a random variable. Idea Pick a real-valued
More informationPhylogeny and systematics. Why are these disciplines important in evolutionary biology and how are they related to each other?
Phylogeny and systematics Why are these disciplines important in evolutionary biology and how are they related to each other? Phylogeny and systematics Phylogeny: the evolutionary history of a species
More informationProbabilities of Evolutionary Trees under a Rate-Varying Model of Speciation
Probabilities of Evolutionary Trees under a Rate-Varying Model of Speciation Mike Steel Biomathematics Research Centre University of Canterbury, Private Bag 4800 Christchurch, New Zealand No. 67 December,
More informationMath 510 midterm 3 answers
Math 51 midterm 3 answers Problem 1 (1 pts) Suppose X and Y are independent exponential random variables both with parameter λ 1. Find the probability that Y < 7X. P (Y < 7X) 7x 7x f(x, y) dy dx e x e
More informationProbability and Measure
Chapter 4 Probability and Measure 4.1 Introduction In this chapter we will examine probability theory from the measure theoretic perspective. The realisation that measure theory is the foundation of probability
More informationLecture 17: Trees and Merge Sort 10:00 AM, Oct 15, 2018
CS17 Integrated Introduction to Computer Science Klein Contents Lecture 17: Trees and Merge Sort 10:00 AM, Oct 15, 2018 1 Tree definitions 1 2 Analysis of mergesort using a binary tree 1 3 Analysis of
More informationBandits, Experts, and Games
Bandits, Experts, and Games CMSC 858G Fall 2016 University of Maryland Intro to Probability* Alex Slivkins Microsoft Research NYC * Many of the slides adopted from Ron Jin and Mohammad Hajiaghayi Outline
More informationRandomized Sorting Algorithms Quick sort can be converted to a randomized algorithm by picking the pivot element randomly. In this case we can show th
CSE 3500 Algorithms and Complexity Fall 2016 Lecture 10: September 29, 2016 Quick sort: Average Run Time In the last lecture we started analyzing the expected run time of quick sort. Let X = k 1, k 2,...,
More informationReconstruction of certain phylogenetic networks from their tree-average distances
Reconstruction of certain phylogenetic networks from their tree-average distances Stephen J. Willson Department of Mathematics Iowa State University Ames, IA 50011 USA swillson@iastate.edu October 10,
More information6 Introduction to Population Genetics
70 Grundlagen der Bioinformatik, SoSe 11, D. Huson, May 19, 2011 6 Introduction to Population Genetics This chapter is based on: J. Hein, M.H. Schierup and C. Wuif, Gene genealogies, variation and evolution,
More information