Fixation Probabilities of Evolutionary Graphs Based on the Positions of New Appearing Mutants

Evolutionary graph theory is a nice measure to implement evolutionary dynamics on spatial structures of populations. To calculate the fixation probability is usually regarded as aMarkov chain process, which is affected by the number of the individuals, the fitness of the mutant, the game strategy, and the structure of the population. However the position of the new mutant is important to its fixation probability. Here the position of the new mutant is laid emphasis on. The method is put forward to calculate the fixation probability of an evolutionary graph (EG) of single level. Then for a class of bilevel EGs, their fixation probabilities are calculated and some propositions are discussed.The conclusion is obtained showing that the bilevel EG is more stable than the corresponding one-rooted EG.


Introduction
The idea of evolution can be dated back to the times of Darwin.In recent years, there is a growing interest in the investigation of evolutionary dynamics on spatial structures of populations; see Lieberman et al. [1], Ohtsuki et al. [2], Nowak et al. (2006), Nie and Zhang [3][4][5], Shakarian et al. [6], Broom et al. [7], Fu et al. [8], and the references mentioned therein.Evolutionary graph theory is a nice measure to implement evolutionary dynamics, in which evolutionary dynamics are arranged on graphs.Individuals in a population are posed on a graph, the weighted edges denoting reproductive rates which govern how often individuals place offspring into adjacent vertices.When a mutant appears in this population, its fixation probability is the probability that this mutant takes over the whole population.After a mutant appears in a population, this mutant may choose to leave the group.The fixation probability of an EG can be affected by the number of the individuals, the fitness of the mutant, the game strategy, and the structure of the population, where the fitness reflects the fit degree of the mutant.A population with lower fixation probability is consequently more stable than a population with higher fixation probability.The stability of a graph refers to the degree of stabilization.Calculating the fixation probability of an EG is a basic and key task.Many researchers have studied the problem.The evolution dynamics on an EG are treated as a discrete time Markov chain process.Some EGs are sensitive to the positions where the new mutants appear.Some papers as Lieberman et al. [1], Shakarian et al. [6], Nowak et al. (2006), Nie [3], and Zhang et al. [9,10] have noticed that the positions of the new appeared mutants in EGs are important to their fixation probabilities.For the directed line, only the mutant appearing at the origin vertex can take over the whole population.The fixation probabilities of some multilevel EGs depend on the appearance of the levels of the new mutants.In this paper, we lay emphasis on the position of the new arisen mutant and put forward a way to calculate the fixation probability of an EG and calculate them for the directed lines, directed circles, and a class of bilevel EGs which is discussed by Zhang et al. [9] and correct their results.
The paper is arranged as follows.Section 2 puts forward the method to calculate the fixation probabilities of singlelevel EGs and applies it to directed lines and directed circles.Section 3 gives the method to calculate the fixation probabilities of a kind of bilevel EGs.Section 4 is the conclusion.

The Fixation Probabilities of Single-Level EGs
The evolutionary graph was initially introduced in 2005 by Lieberman et al. [1] and significantly developed by Nie and Zhang [3][4][5] Proof.Let  be the event that a new mutant appears in the population and   ( = 1, 2, . . ., ) the event that the new mutant appears at the position of vertex .So  1 ,  2 , . . .,   are inconsistent and the following equation is true: Let  be the event that the new mutant is fixed; then we have the following equations: Since the new mutant can appear at any vertices of the population, the vertices are inconsistent, and the individuals are treated identically; the events  1 ,  2 , . . .,   are inconsistent such that  1 ,  2 , . . .,   are inconsistent.Therefore the second equation is satisfied, (  ) = 1/( = 1, 2, . . ., ), and the result is obtained.The proof is complete.
We have given the statement of the calculation of fixation probabilities of EGs.In the following we apply it to the directed lines and directed circles.Proposition 2. A directed line with  vertices named 1, 2, . . .,  is shown as follows: Suppose a new mutant with fitness  is introduced and all the residents are with fitness 1.Then the fixation probability of the new mutant is 1/.
Proof.It is obvious that the new mutant will occupy the whole population if it appears at position 1; the new mutant will disappear if it appears at any other positions.So we have  (1) = 1 and   = 0,  = 2, 3, . . ., .From (1) we have The proof is complete.
We have applied Theorem 1 to the directed lines and we have gotten the same result as Nowak [11] did.In the following we will apply Theorem 1 to directed circles.Proposition 3. A directed circle with  vertices named 1, 2, . . .,  is shown in Figure 1.Suppose a new mutant with fitness  is introduced and all the residents are with fitness 1.
The result in Proposition 3 is the same as that Nowak [11] obtained, but we lay emphasis on the position of the new mutant.
In the section, we gave the calculation of fixation probability which has a close relation with the vertex where the new mutant appears in Theorem 1 and applied it to the directed lines and directed circles.From the above discussion, we can get the conclusion that the directed line cannot be regarded a Markov chain process, for it is sensitive to the position where the new mutant is introduced, while the directed circle can be looked upon as a Markov chain process, for it is insensitive to the position.In the following section, a class of bilevel EGs will be discussed.

The Fixation Probabilities of a Class of Bilevel EGs
Nie [3] introduced bilevel EGs: the upper levels are stars whereas the lower levels are isothermal EGs.Nie showed that the fixation probabilities of the EGs on two levels are lower than those of the star structure, such that the EGs on two levels are more stable.Nie explained the phenomenon of autoeciousness.Zhang et al. [9] discussed a class of bilevel EGs: the upper levels are isothermal and the lower levels are one rooted.The bilevel EGs reflect hierarchical population structures.The work on bilevel EGs is extended by Zhang et al. [10] where the considered bilevel EGs have different fitness depending on the levels of the graphs.Zhang et al. [10] described several anecdotal applications of bilevel EGs, including the budgerigar of Oceania and the Aptenodytes forsteri (emperor penguins) of Antarctica.Zhou [12] took a step further toward a real-world application by using bilevel EGs and star EGs to examine the stability of various types of business forms.They model corporations with individual decisions (CIDs) as 1-level star graphs and multiperson decision corporations (MDCs) and stock corporations (SCs) are modeled as bilevel graphs.They find that, under reasonable conditions, MDCs have higher fixation probabilities than CIDs, which have higher fixation probabilities than SCs.Hence, through the lens of EGT, SCs represent the most stable organizational structure for business among these three types of business forms.The bilevel EGs and multilevel EGs are important.There are many works to be done in this region.
Here we give and correct the fixation probabilities of the bilevel EGs discussed in Zhang et al. [9].For the bilevel EG as shown in Figure 2, the vertices can be divided into two sets.One is named ; the other is named .The vertices in set  are called the leaders with number , while the vertices in set  are called the followers with number .The leaders form an isothermal .When the leaders (set ) are collapsed into a single vertex, we have a one-rooted EG.In the following theorem, we give the fixation probability of the described bilevel EG.
Theorem 4. For a bilevel EG, the upper level is isothermal with  vertices whereas the lower level is one-rooted with  A Figure 2: A bilevel EG where the leaders (set ) form an isothermal EG.When they are collapsed into a single vertex, we have a onerooted EG [6].
vertices.If a new mutant with fitness  is introduced and all the residents are with fitness 1 in the EG, its fixation probability is Proof.Let  be the event that the new mutant appears in set ,  the event that the new mutant appears in set  which is the opposite event of ,  the event that new mutant occupies set , and  the event that the new mutant occupies the whole vertices.If a mutant occupies the whole vertices, it must occupy all the vertices in set . Therefore we have  ⊆  such that  = .The fixation probability of a new mutant refers to the probability that the new mutant occupies the whole EG such that  = ().Since vertices are treated identically, the probability that the new mutant appears in set  is /( + ) and the probability that the new mutant appears in set  is /( + ).The proof is complete.
The fixation probability of the bilevel EG is obtained, which corrects the result of Zhang et al. [9].In the following we will discuss its property.Proposition 5. Let  be the fixation probability of the bilevel EG with  leaders and  followers and let  be the fitness of the new mutant.When all the leaders in the upper level are collapsed into a single vertex, a one-rooted EG with  + 1 vertices is obtained.The fixation probability of the one-rooted EG with  + 1 vertices is denoted by  1 .If the new mutant is advantageous, that is,  > 1, one can get  <  1 .
Proof.Since  1 is the fixation probability of a one-rooted EG with  + 1 vertices,  1 = 1/( + 1). therefore The proof is complete.
In this section, we give the fixation probability of a bilevel EG and discuss its proposition.We get the result that the bilevel EG has lower fixation probability than that of the corresponding one-rooted EG when the leaders are collapsed into one vertex.The same conclusion is obtained as that in Nie [3] and Zhang et al. [9].

Conclusion Remark
In this paper, we lay emphasis on the effect of the position of the new mutant on its fixation probability.We give the fixation probability of an EG of single level.The result is applied to directed lines and directed circles.For a directed line, the fixation probability of a new mutant has a close relation with the position where it appears.Only the mutant appearing at the "first" position will occupy the whole population while the mutant appearing at other positions will disappear.So the directed line is sensitive to the position and the evolutionary dynamics on the directed line cannot be regarded as a Markov chain process.For a directed circle, the fixation probability of the new mutant appearing at each position is identical.So the directed circle is insensitive to the position and the evolutionary dynamics on the directed circle can be regarded as a Markov chain process.
Then the idea is applied to a class of bilevel EGs.The upper levels are isothermal while the lower levels are one rooted.The fixation probabilities are given which correct the result in Zhang et al. [9].We compare the bilevel EG with the case where all the leaders in the upper level are collapsed into only one vertex, that is, the one-rooted EG.For an advantage mutant  > 1, we get the result that the bilevel EG has lower fixation probability such that the bilevel EG is more stable.
Zhou [12] modeled a CID as a star graph and an MDC as a bilevel EG when the lower level and the upper level are both isothermal.Here we would like to model a CID as a onerooted graph and an MDC as a bilevel EG when the upper level is isothermal and the lower level is one rooted.From the conclusion of Proposition 5 we know that the fixation probability of the bilevel graph is lower than that of the corresponding one-rooted graph such that the bilevel EG is more stable.The result, however, is accordant with Zhou [12].

Figure 1 :
Figure 1: A directed circle with  individuals.