
Supertree methods are tree reconstruction techniques that combine several smaller gene trees (possibly on different sets of species) to build a larger species tree. The question of interest is whether the reconstructed supertree converges to the true species tree as the number of gene trees increases (that is, the consistency of supertree methods). In this paper, we are particularly interested in the convergence rate of the maximum likelihood supertree. Previous studies on the maximum likelihood supertree approach often formulate the question of interest as a discrete problem and focus on reconstructing the correct topology of the species tree. Aiming to reconstruct both the topology and the branch lengths of the species tree, we propose an analytic approach for analyzing the convergence of the maximum likelihood supertree method. Specifically, we consider each tree as one point of a metric space and prove that the distance between the maximum likelihood supertree and the species tree converges to zero at a polynomial rate under some mild conditions. We further verify these conditions for the popular exponential error model of gene trees.
Citation: Vu Dinh, Lam Si Tung Ho. Convergence of maximum likelihood supertree reconstruction[J]. AIMS Mathematics, 2021, 6(8): 8854-8867. doi: 10.3934/math.2021513
[1] | Ahmad Mohammed Alghamdi, Sadek Gala, Maria Alessandra Ragusa . A regularity criterion of weak solutions to the 3D Boussinesq equations. AIMS Mathematics, 2017, 2(3): 451-457. doi: 10.3934/Math.2017.2.451 |
[2] | Xinli Wang, Haiyang Yu, Tianfeng Wu . Global well-posedness and optimal decay rates for the n-D incompressible Boussinesq equations with fractional dissipation and thermal diffusion. AIMS Mathematics, 2024, 9(12): 34863-34885. doi: 10.3934/math.20241660 |
[3] | Feng Cheng . On the dissipative solutions for the inviscid Boussinesq equations. AIMS Mathematics, 2020, 5(4): 2869-2876. doi: 10.3934/math.2020184 |
[4] | Xiaochun Sun, Yulian Wu, Gaoting Xu . Global well-posedness for the 3D rotating Boussinesq equations in variable exponent Fourier-Besov spaces. AIMS Mathematics, 2023, 8(11): 27065-27079. doi: 10.3934/math.20231385 |
[5] | Ru Bai, Tiantian Chen, Sen Liu . Global stability solution of the 2D incompressible anisotropic magneto-micropolar fluid equations. AIMS Mathematics, 2022, 7(12): 20627-20644. doi: 10.3934/math.20221131 |
[6] | Sen Ming, Jiayi Du, Yaxian Ma . The Cauchy problem for coupled system of the generalized Camassa-Holm equations. AIMS Mathematics, 2022, 7(8): 14738-14755. doi: 10.3934/math.2022810 |
[7] | Tariq Mahmood, Zhaoyang Shang . Blow-up criterion for incompressible nematic type liquid crystal equations in three-dimensional space. AIMS Mathematics, 2020, 5(2): 746-765. doi: 10.3934/math.2020051 |
[8] | Ailing Ban . Asymptotic behavior of non-autonomous stochastic Boussinesq lattice system. AIMS Mathematics, 2025, 10(1): 839-857. doi: 10.3934/math.2025040 |
[9] | Shang Mengmeng . Large time behavior framework for the time-increasing weak solutions of bipolar hydrodynamic model of semiconductors. AIMS Mathematics, 2017, 2(1): 102-110. doi: 10.3934/Math.2017.1.102 |
[10] | Yasir Nadeem Anjam . The qualitative analysis of solution of the Stokes and Navier-Stokes system in non-smooth domains with weighted Sobolev spaces. AIMS Mathematics, 2021, 6(6): 5647-5674. doi: 10.3934/math.2021334 |
Supertree methods are tree reconstruction techniques that combine several smaller gene trees (possibly on different sets of species) to build a larger species tree. The question of interest is whether the reconstructed supertree converges to the true species tree as the number of gene trees increases (that is, the consistency of supertree methods). In this paper, we are particularly interested in the convergence rate of the maximum likelihood supertree. Previous studies on the maximum likelihood supertree approach often formulate the question of interest as a discrete problem and focus on reconstructing the correct topology of the species tree. Aiming to reconstruct both the topology and the branch lengths of the species tree, we propose an analytic approach for analyzing the convergence of the maximum likelihood supertree method. Specifically, we consider each tree as one point of a metric space and prove that the distance between the maximum likelihood supertree and the species tree converges to zero at a polynomial rate under some mild conditions. We further verify these conditions for the popular exponential error model of gene trees.
Assume that ℧ is a non-empty, closed and convex subset of a real Banach space Ω. The normalized duality mapping J:Ω→2Ω∗ (Ω∗ is the dual space of Ω) is defined by
J(ν)={ϰ∈Ω∗:⟨ν,ϰ⟩=‖ν‖2=‖ϰ‖2}. |
A mapping Z:℧→℧ is called nonexpansive if
‖Zν−Z℘‖≤‖ν−℘‖, ∀ν,℘∈℧. |
Here, we denote ∇(Z) by the set of fixed points of Z, that is ∇(Z)={ν∈℧:Zν=ν} and we consider ∇(Z)≠∅.
One of the strong contributions with important applications to fixed point theory is Banach contraction principle, which states:
Theorem 1.1. [1] Every contraction mapping Z:ℏ→ℏdefined on a complete metric space (ℏ,σ) has a unique fixedpoint, where σ is the distance that describes the mapping Z, i.e.,
σ(Zν,Z℘)≤μσ(ν,℘),∀ν,℘∈ℏ,μ<1. | (1.1) |
Moreover, for arbitrary ν0∈ℏ, the sequence {νm}created by
νm+1=Zνm,m≥0, | (1.2) |
converges strongly to the unique fixed point.
It should be noted that a mapping Z verifying (1.1) is called a strict contraction and if μ=1 in (1.2), then it is called nonexpansive. The iterative sequence (1.2) is due to Picard [2]. For the iterative formula, it was observed that if the condition μ<1 on the operator Z is weakened to μ=1, the sequence {νn} defined by (1.2) may fail to converge to a fixed point of Z. To overcome this shortcoming, Krasnoselskii [3], replaced Picard iteration formula by the following formula:
ν0∈ℏ, νm+1=12(νm+Zνm), m≥0. |
He proved that the iterative sequence converges to the fixed point.
In the light of [3], the successful iterative method presented and known as Krasnoselskii-Mann iterative scheme and formulated as follows:
ν0∈℧, νm+1=(1−ηm)νm+ηmZνm, m≥0, | (1.3) |
where {ηm} is a sequence of non-negative real numbers in (0,1). It was observed that via the stipulation ∇(Z)≠∅ and mild assumptions forced on {ηm}, the sequence {νm} generated by (1.3) converges weakly to a fixed point of Z.
Krasnoselskii-Mann algorithm is one of many successful iteration schemes for approximating fixed points of nonexpansive mappings. It provides a unified framework for many algorithms in various disciplines, so the following approach is important.
Theorem 1.2. [4] Assume that Z is a nonexpansive mapping on a realHilbert space ℸ and ∇(Z)≠∅. Then the sequence {νn} maked by (1.3) converges weakly to a fixed point of Z, provided that ηm∈[0,1] and ∞∑m=0ηm=∞.
It should be remarked that all previous contributions on Krasnoselskii-Mann algorithm for nonexpansive mappings have only weak convergence even in a real HS, see [4]. Further, Krasnoselskii-Mann algorithm was generalized by Yang and Zhao [5]. They introduced some important theorems about it and they called their theorems KM theorems.
Bruck [6] noted that the importance of studying nonexpansive mappings lies in two main reasons:
i) nonexpansive mappings are closely related to the monotonicity methods that were updated in the early 1960s and constitute one of the first classes of nonlinear mapping to be treated using the fixed point technique by studying the exact geometrical properties of the basic Banach spaces rather than the compactness properties.
ii) nonexpansive mappings assignments in applications appear as transitional parameters for initial value problems of differential inclusions in the form 0∈dμdτ+Υ(τ)μ=0, where a set-valued operators {Υ(τ)}, are accretive or minimally continuous and dissipative.
In nonlinear mapping theory and its applications, building fixed point for nonexpansive mappings assignments is a very important topic, especially, in signal processing and image recovery, see [7,8,9]. Study of Krasnoselskii-Mann iterative procedures to approximate fixed points of nonexpansive mappings assignments and fixed points of some of their generalizations and approximate zeros of operators of accretive-type has become more widespread and prosperous over the past thirty years or so, for further clarification, we would like to guide the reader to [10,11,12,13,14].
Very recently, a new form for Mann's algorithm is proposed by Bot et al. [15] to overcome the deficiency described before, and he described it as follows: let φ0 be arbitrary in ℸ, for all m≥0,
νm+1=ηmνm+1+ζm(Z(ηmνm)−ηmνm), | (1.4) |
they showed that the iterative sequence {νn} generated by (1.4) is strongly convergent via suitable assumptions for {ηm} and {ζm}. A sequence {ζm} in (1.2) has an effective role in acceleration, it called Tikhonov regularization sequence. Many theoretical and numerical discussions to study strong convergence using Tikhonov regularization methodology have been presented by [18,19,20].
Recently newer types of algorithms have been developed and introduced, such as the inertia algorithm first introduced by Polyak [18]. He used an inertial extrapolation methodology for minimizing a smooth convex function. It is worth noting that, these simple changes affected positively in the performance and effectiveness of these algorithms.
After adopting this concept, researchers were able to implicate additional terms to the inertial algorithm to delve into the study of many vital applications, for example, but not limited to, inertial extragradient algorithms [19,20,21], inertial projection algorithms [22,23,24,25,26], inertial Mann algorithms [27] and inertial forward-backward splitting algorithms [28,29]. There is no doubt that these algorithms are significantly faster than the inertial algorithms.
Based on the above work, in the present manuscript, Ishikawa algorithm has been developed by adding the term of the inertial to obtain an advanced algorithm, called a modified inertial Ishikawa algorithm. A strong convergence using the proposed algorithm is also discussed in a real uniformly convex Banach space with uniformly Gâteaux differentiable norm. Moreover, as an application, we find zeros of accretive mappings. Our results generalize and extend many of the older findings in this regard. Finally, two numerical experiments are given to illustrate the behavior of the purposed algorithm.
Assume that Ω is a real normed linear space and assume ℵ={ν∈Ω:‖ν‖=1}. We say that Ω have a Gâteaux differentiable norm if the limit below exists for all ν,℘∈ℵ,
limτ→∞‖ν+τ℘‖−‖ν‖τ, |
and Ω is called smooth. Furthermore, we say that Ω has a uniformly Gâteaux differentiable norm, if for any ℘∈ℵ the limit is attained uniformly for ν∈ℵ. Also, Ω is called uniformly smooth if the limit exists uniformly for (ν,℘)∈ℵ. It is obvious that any duality mapping on Ω is a single-valued if Ω is smooth and if Ω has a uniformly Gâteaux differentiable norm then the duality mapping is norm-to-weak∗ uniformly continuous on bounded subsets of Ω.
Suppose that Δ is a non-empty, closed, convex and bounded subset of a real Banach space Ω and let d(Δ)=sup{‖ν−℘‖, ν,℘∈Δ} refer to the diameter of Δ, and for w(Δ)=inf{w(ν,℘), ν∈Δ} refer to the Chebyshev radius of Δ relative to itself, where for ν∈Δ, w(ν,Δ)=sup{‖ν−℘‖, ℘∈Δ}. The normal structure coefficient N(Ω) of Ω introduced by Bynum [30] as follows: Let Δ be a non-empty, closed, convex, and bounded subset of Ω, then N(Ω) is defined by N(Ω)=inf{d(Δ)w(Δ):d(Δ)>0}. If N(Ω)>1, then the space Ω has a uniform normal structure. It should be noted that, every space with a uniform normal structure is reflexive, this implies that all uniformly convex and uniformly smooth Banach spaces have a uniform normal structure, for example, see, [11,31].
The lemmas below are very important in the sequel.
Lemma 2.1. [32] Let Ω be a real uniformly convex Banach space.For arbitrary u>0, assume that ℵu(0)={ν∈Ω:‖ν‖≤u} and α∈[0,1]. Thenthere is a continuous strictly increasing convex function r:[0,2u]→R, r(0)=0 so that the inequality below holds
‖αν+(1−α)℘‖2≤α‖ν‖2+(1−α)‖℘‖2−α(1−α)r(‖ν−℘‖). |
Lemma 2.2. Let Ω be a real normed linear space, then for all ν,℘∈Ω, j(ν+℘)∈J(ν+℘), we have
‖ν+℘‖2≤‖ν‖2+2⟨℘,j(ν+℘)⟩. |
Lemma 2.3. [33] Let Ω be a uniformly convex Banach space, Δ be a non-empty, closed and convex, subset of Ω and Z:Δ→Δ be a nonexpansive mapping with a fixed point.Suppose that the sequence {νm} in Δ is so that νm⇀ν and νm−Zνm⟶℘. Then ν−Zν=℘.
Lemma 2.4. [31] Assume that Ω is a Banach space with uniformnormal structure, Δ is a nononexpansive mapping bounded subset of Ω and Z:Δ→Δ is a uniformly L−Lipschitzianmapping with L<N(Ω)12. Consider there is a non-emptybounded closed convex subset ℜ of Δ with the property (D) below:
ν∈ℜ⇒ϖw(ν)∈ℜ. |
Then Z has a fixed point in Δ.
Note: ϖw(ν) here is the ϖ−limit set of Z at ν, that is, the set {℘∈Ω:y=weak ϖ−limZnjν,for some nj→∞}.
Lemma 2.5. [34] Assume that (ν0,ν1,ν2,...)∈l∞, is so that δmνm≤0 for all Banach limits δ. If lim supm→∞(νm+1−νm)≤0, then lim supm→∞νm≤0.
Lemma 2.6. [35] Suppose that {en} is a sequence of non-negativereal numbers verifying the inequality below
em+1≤(1−cm)em+fmσm+πm,m≥1, |
if
● {cm}⊂[0,1], ∑cm=∞;
● lim supm→∞σm≤0;
● for each m≥0, πm≥0, ∑πm<∞.
Then, limm→∞em=0.
Under mild conditions, in this section, we shall discuss the strong convergence of a modified inertial Ishikawa algorithm for nonexpansive mappings.
Theorem 3.1. Let Ω be a real uniformly convex Banach space withuniformly Gâteaux differentiable norm. Suppose that Z:Ω→Ω is a nonexpansive mapping so that ∇(Z)≠∅. Consider the hypotheses below hold:
(i) limm→∞σm=0, ∞∑m=1σm=∞, σm∈(0,1), ρm∈[ℓ1,ℓ2]⊂(0,1),
(ii) limm→∞πm=0, πm∈(0,1) and ∞∑m=0πm‖νm−νm−1‖<∞.
Set ν0,ν1 arbitrary. Let the sequence {νm} createditeratively by
{ℏm=νm+πm(νm−νm−1),℘m=(1−σm)ℏm,νm+1=(1−ρm)℘m+ρmZ℘m,m≥1. | (3.1) |
Then the sequence {νm} converges strongly to a point in ∇(Z).
Proof. For any d∈∇(Z), by (3.1), we have
‖νm+1−d‖=‖(1−ρm)(℘m−d)+ρm(Z℘m−d)‖≤(1−ρm)‖℘m−d‖+ρm‖Z℘m−d‖=(1−ρm)‖℘m−d‖+ρm‖Z℘m−Zd‖≤(1−ρm)‖℘m−d‖+ρm‖℘m−d‖≤‖℘m−d‖=‖(1−σm)ℏm−d‖≤(1−σm)‖ℏm−d‖+σm‖d‖≤(1−σm)‖(νm−d)+πm(νm−νm−1)‖+σm‖d‖≤(1−σm)‖νm−d‖+(1−σm)πm‖νm−νm−1‖+σm‖d‖≤max{‖νm−d‖,‖νm−νm−1‖,‖d‖}. |
By mathematical induction, it is easy to see that
‖νm−d‖≤max{‖ν1−d‖,‖ν1−ν0‖,‖d‖}, ∀m≥0. |
This shows that {νm} is bounded and also are {ℏm} and {℘m}.
Based on Lemmas 2.1, 2.2 and Algorithm (3.1), one sees that
‖νm+1−d‖2=‖(1−ρm)(℘m−d)+ρm(Z℘m−d)‖2≤(1−ρm)‖℘m−d‖2+ρm‖Z℘m−d‖2−ρm(1−ρm)r(‖Z℘m−℘m‖)≤(1−ρm)‖℘m−d‖2+ρm‖℘m−d‖2−ρm(1−ρm)r(‖Z℘m−℘m‖)=‖℘m−d‖2−ρm(1−ρm)r(‖Z℘m−℘m‖)≤‖ℏm−d‖2+2σm⟨ℏm−d,j(℘m−d)⟩−ρm(1−ρm)r(‖Z℘m−℘m‖)≤‖νm−d‖2+2πm⟨νm−d,j(ℏm−d)⟩+2σm⟨ℏm−d,j(℘m−d)⟩−ρm(1−ρm)r(‖Z℘m−℘m‖). |
On the other hand, one can write
ρm(1−ρm)r(‖Z℘m−℘m‖)≤‖νm−d‖2−‖νm+1−d‖2+2πm⟨νm−d,j(ℏm−d)⟩+2σm⟨ℏm−d,j(℘m−d)⟩. | (3.2) |
The boundedness of {νm}, {ℏm} and {℘m} leads to there are constants Λ1, Λ2>0 so that
⟨νm−d,j(ℏm−d)⟩≤Λ1 and ⟨ℏm−d,j(℘m−d)⟩≤Λ2 for all m≥1. | (3.3) |
Applying (3.3) in (3.2), we have
ρm(1−ρm)r(‖Z℘m−℘m‖)≤‖νm−d‖2−‖νm+1−d‖2+2πmΛ1+2σmΛ2. | (3.4) |
In order to obtain the strong convergence, we discuss the following cases:
Case (a).If the sequence {‖νm−d‖} is monotonically decreasing, then {‖νm−d‖} is convergent. It is easy to see that
‖νm+1−d‖2−‖νm−d‖2→0, |
as m→∞, this leads to directly by (3.4),
ρm(1−ρm)r(‖Z℘m−℘m‖)→0. |
By the property of r and since ρm∈[ℓ1,ℓ2]⊂(0,1), we have
‖Z℘m−℘m‖→0. | (3.5) |
Combining (3.1) and (3.5), we have
‖νm+1−℘m‖=ρm(Z℘m−℘m)→0. | (3.6) |
It follows from (3.1) and condition (i) that
‖℘m−ℏm‖=σm‖ℏm‖→0. | (3.7) |
By condition (ii), we get
‖ℏm−νm‖=πm‖νm−νm−1‖→0. | (3.8) |
Based on (3.7) and (3.8), we can write
‖℘m−νm‖≤‖℘m−ℏm‖+‖ℏm−νm‖→0. | (3.9) |
Using (3.6) and (3.9), we have
‖νm+1−νm‖≤‖νm+1−℘m‖+‖℘m−νm‖→0 as m→∞. |
From (3.5), (3.7) and (3.8), we get
‖Zνm−νm‖≤‖Zνm−Z℘m‖+‖Z℘m−℘m‖+‖νm−℘m‖≤2‖νm−℘m‖+‖Z℘m−℘m‖≤2(‖℘m−ℏm‖+‖ℏm−νm‖)+‖Z℘m−℘m‖→0. |
Since {νm} is bounded, then there exists the subsequence {νmb}⊂{νm} so that it converges weakly to d∈Ω. Furthermore, Lemma 2.3 implies that d∈∇(Z).
Now, we shall show that
lim supm→∞⟨−d,j(℘m−d)⟩≤0. |
For this, define a map χ:Ω→R by
χ(ν)=δm‖℘m−ν‖2, ∀ν∈Ω. |
Then, χ(ν)→∞ as ‖ν‖→∞, χ is convex and continuous. As Ω is reflexive, then there is ℘∗∈Ω so that χ(℘∗)=mina∈Ωχ(a). Thus, the set ˆℜ≠∅, where ˆℜ is defined as
ˆℜ={ν∈Ω:χ(ν)=mina∈Ωχ(a)}. |
Again, since limm→∞‖Z℘m−℘m‖=0, then by induction, we can see that limm→∞‖Zn℘m−℘m‖=0 for all n≥1. Hence, from Lemma 2.4, if ν∈ℜ and ℘=ϖ−limj→∞Znjν, then from weak lower semi-continuity of χ and limm→∞‖Z℘m−℘m‖=0, we get (since limm→∞‖Z℘m−℘m‖=0 implies limm→∞‖Zn℘m−℘m‖=0, n≥1, this is easily proved by induction)
χ(℘)≤lim infj→∞χ(Znjν)≤lim supn→∞χ(Znν)=lim supn→∞(δm‖℘m−Znν‖2)=lim supn→∞(δm‖℘m−Z℘m+Z℘m−Znν‖2)≤lim supn→∞(δm‖Z℘m−Znν‖2)≤lim supn→∞(δm‖℘m−ν‖2)=χ(ν)=infa∈Ωχ(a). |
Thus, ℘∗∈ˆℜ. Therefore by Lemma 2.4, Z has a fixed point in ˆℜ and so ˆℜ∩∇(Z)≠∅. As a special case without losing the general case, suppose that ℘∗=d∈ˆℜ∩∇(Z). Consider τ∈(0,1). Then it is easy to see that χ(d)≤χ(d−τd) with the helping of Lemma 2.2, one sees that
‖℘m−d+τd‖2≤‖℘m−d‖2+2τ⟨d,j(℘m−d+τd)⟩, |
by the properties of χ, we can write
1δmχ(d−τd)≤1δmχ(d)+2τ⟨d,j(℘m−d+τd)⟩, |
By arranging the above inequality, we have
2τδm⟨−d,j(℘m−d+τd)⟩≤χ(d)−χ(d−τd)≤0. |
This leads to
δm⟨−d,j(℘m−d+τd)⟩≤0. |
Moreover,
δm⟨−d,j(℘m−d)⟩≤δm⟨−d,j(℘m−d)−j(℘m−d+τd)⟩+δm⟨−d,j(℘m−d+τd)⟩≤δm⟨−d,j(℘m−d)−j(℘m−d+τd)⟩. | (3.10) |
Since the normalized duality mapping is norm-to-weak∗ uniformly continuous on bounded subsets of Ω, then we have, as τ→0 and for fixed n,
⟨−d,j(℘m−d)−j(℘m−d+τd)⟩≤⟨−d,j(℘m−d)⟩−⟨−d,j(℘m−d+τd)⟩→0. |
Thus, for each ϵ>0, there is ςϵ>0 so that for all τ∈(0,ςϵ),
⟨−d,j(℘m−d)⟩−⟨−d,j(℘m−d+τd)⟩<ϵ. |
Hence,
δm⟨−d,j(℘m−d)⟩−δm⟨−d,j(℘m−d+τd)⟩≤ϵ. |
Because ϵ is an arbitrary, then by (3.10), one can obtain
δm⟨−d,j(℘m−d)⟩≤0. |
By triangle inequality, we have
‖℘m+1−℘m‖≤‖℘m+1−νm+1‖+‖νm+1−℘m‖. |
According to (3.6) and (3.9), we get
limm→∞‖℘m+1−℘m‖=0. |
Since the normalized duality mapping is norm-to-weak∗ uniformly continuous on bounded subsets of Ω, then we have
limm→∞(⟨−d,j(℘m−d)⟩−⟨−d,j(℘m+1−d)⟩)=0. |
It follows from Lemma 2.5 that
lim supm→∞⟨−d,j(℘m−d)⟩≤0. |
Ultimately, from (3.1), Stipulation (ii) and Lemma 2.2, we have
‖νm+1−d‖2=‖(1−ρm)(℘m−d)+ρm(Z℘m−d)‖2≤(1−ρm)‖℘m−d‖2+ρm‖Z℘m−d‖2≤‖℘m−d‖2=‖(1−σm)(ℏm−d)−σmd‖2=(1−σm)‖ℏm−d‖2+2σm⟨−d,j(℘m−d)⟩≤(1−σm)‖(νm−d)+πm(νm−νm−1)‖2+2σm⟨−d,j(℘m−d)⟩≤(1−σm)‖νm−d‖2+2πm⟨νm−νm−1,j(ℏm−d)⟩+2σm⟨−d,j(℘m−d)⟩=(1−σm)‖νm−d‖2+2σm⟨−d,j(℘m−d)⟩. | (3.11) |
Applying Lemma 2.6, we conclude that, {νm}→d∈∇(Z).
Case (b).If the sequence {‖νm−d‖} is not monotonically decreasing. Put Ξm=‖νm−d‖2 and assume that Π:N→N is a mapping defined by
Π(m)=max{ℏ∈N:ℏ≤m, Ξℏ≤Ξℏ+1}. |
Obviously, Π is a non-decreasing sequence so that limm→∞Π(m)=∞ and ΞΠ(m)≤ΞΠ(m)+1 for m≥m0 (for some m0 large enough). Based on (3.4), one sees that
ρΠ(m)(1−ρΠ(m))r(‖Z℘Π(m)−℘Π(m)‖)≤‖νΠ(m)−d‖2−‖νΠ(m)+1−d‖2+2πΠ(m)Λ1+2σΠ(m)Λ2=ΞΠ(m)−ΞΠ(m)+1+2πΠ(m)Λ1+2σΠ(m)Λ2≤2πΠ(m)Λ1+2σΠ(m)Λ2→0 as m→∞. |
Furthermore, we get
‖Z℘Π(m)−℘Π(m)‖→0 as m→∞. |
By following the same scenario in Case (a) we can prove that {νΠ(m)}⇀d as Π(m)→∞ and lim supΠ(m)→∞⟨−d,j(℘Π(m)−d)⟩≤0. For all m≥m0, we obtain by (3.11) that
0≤‖νΠ(m)+1−d‖2−‖νΠ(m)−d‖2≤σΠ(m)[2⟨−d,j(℘Π(m)−d)⟩−‖νΠ(m)−d‖2], |
this implies that
‖νΠ(m)−d‖2≤2⟨−d,j(℘Π(m)−d)⟩. |
Since lim supΠ(m)→∞⟨−d,j(℘Π(m)−d)⟩≤0, then we have after taking the limit as m→∞ in the above inequality,
limm→∞‖νΠ(m)−d‖2=0. |
Hence
limm→∞ΞΠ(m)=limm→∞ΞΠ(m)+1=0. |
Moreover, for all m≥m0, it is easy to notice that Ξm≤ΞΠ(m)+1 if m≠Π(m) (that is, Π(m)<m), since Ξi>Ξi+1 for Π(m)+1≤i≤m. As a result, for all m≥m0, we get
0≤Ξm≤max{ΞΠ(m),ΞΠ(m)+1}=ΞΠ(m)+1. |
Thus, limm→∞Ξm=0, this conclude that {νm} converges strongly to a point d. This finishes the proof.
Remark 3.2. (r1) Here, the results of Tan and Cho [36] are generalized from a real HS to a real uniformly convex Banach space with uniformly Gâteaux differentiable norm.
(r2) Because of the wide applications in most branches of mathematics and engineering for the problem of finding fixed points of nonexpansive mappings, it has attracted the attention of many researchers.
(r3) Since every uniformly smooth Banach space has uniformly G âteaux differentiable norm. Then, our theorem can be stated in a uniformly convex Banach space which is also uniformly smooth. (Corollary 3.3).
Corollary 3.3. Let Ω be a real uniformly convex Banach space which isalso uniformly smooth. Assume that Z:Ω→Ω is anonexpansive mapping so that ∇(Z)≠∅. Let {νm}be a sequence created iteratively by (3.1). Then the sequence {νm} converges strongly to a point in ∇(Z).
Let Ω be a real uniformly convex Banach space with uniformly Gâteaux differentiable norm. We say that a mapping Υ:D(Υ)→Ω that has a domain D(Υ) is accereative if there is j(ν−℘)∈J(ν−℘) such that
⟨j(ν−℘),Υν−Υ℘)⟩≥0, for ν,℘∈D(Υ). | (4.1) |
According to Inequality (4.1), Kato [37] introduced another definition of the accereative mapping as follows: A mapping Υ is called accereative if the inequality below holds
‖ν−℘‖≤‖ν−℘+b(Υν−Υ℘)‖ ∀s>0 and for each ν,℘∈D(Υ). | (4.2) |
We must recall that accerative operators are monotone if Ω is a Hilbert space. Moreover, if Υ is accretive and its range is R(I+eA)=Ω, for all e>0, then Υ is called m−accretive. Also, if ¯D(Υ)⊆R(I+eΥ) for all e>0, then Υ is said to satisfy the range condition, where ¯D(Υ) is the closure of the domain of Υ. Furthermore, if Υ is accerative [38], then the mapping JΥ:R(I+Υ)→D(Υ), which defined by JΥ=(I+Υ)−1 is a single-valued nonexpansive and ∇(JΥ)=N(Υ), where N(Υ)={ν∈D(Υ):Υν=0} and ∇(JΥ)={ν∈Ω:JΥν=ν}.
In 1967, the accerative operators are presented independently by Browder [39] and Kato [37]. The study of such mappings is extremely interesting because of their firm link with the existence theory for nonlinear equations of evolution in Banach spaces.
Accerative operators are heavily involved under a suitable Banach space in many physically significant problems where these problems can be formulated as an initial boundary value problem of the form
dμdτ+Υμ=0, μ(0)=μ0. | (4.3) |
There are several embedded models of evolution equations such as Schrö dinger, heat and wave equation [40]. Heavy work on the theory of accretive operators has been published by Browder [39] explains that if Υ is locally Lipschitzian and accretive on Ω, then Problem (4.3) has a solution. Also, under the same conditions and the existence result of (4.3), he proved that Υ is m− accretive and there is a solution to the equation below
Υμ=0. | (4.4) |
By Ray [40], Browder's results are elegantly and refined using fixed point theory of Caristi [41]. Martin [42] generalized the results of Browder by proving that in the space Ω. Problem (4.3) is solvable if Υ is continuous and accretive. Moreover, he showed that if Υ is continuous and accretive, then Υ is m−accretive. For more details about theorems for zeros of accretive operators see Browder [43] and Deimling [44].
It should be noted that, if μ is independent of τ in Eq (4.3), then dμdτ=0. Thus, Eq (4.3) reduces to (4.4) whose solution describes the stable or the equilibrium state of the problem created by (4.3). This in turn is very exciting in many elegant applications such as, to name but a few, economics, physics and ecology. As a result, strenuous efforts have been made to solve Eq (4.4) when Υ is accretive. Because Υ, in general, is nonlinear, there is no known way to find a close solution to this equation, and this is what made researchers interested in studying the fixed point and approximate iterative methods for zeros of m− accretive mappings. So it became a thriving area for research to the present time.
In this part, we involve the results of Theorem 3.1 to describe applications of the above results to finding zeros of accretive mappings. Recall, we assume that Ω is a real uniformly convex Banach space with uniformly Gâteaux differentiable norm and consider Υ:Ω→Ω is continuous and accretive mapping. We will find a solution to the equation:
find ν∈Ω so that Υν=0. | (4.5) |
Now the statements and proof of the theorem for finding the solution to Eq (4.5) are fit for presentation.
Theorem 4.1. Let Ω be a real uniformly convex Banach space withuniformly Gâteaux differentiable norm. Assume that Υ:Ω→Ω is a continuous and accretive mapping so that N(Υ)≠∅. Let the sequence {νm} createditeratively by ν0,ν1∈Ω,
{ℏm=νm+πm(νm−νm−1),℘m=(1−σm)ℏm,νm+1=(1−ρm)℘m+ρmJΥ℘m,m≥1, |
where JΥ=(I+Υ)−1. Then the sequence {νm}converges strongly to a point in N(Υ), provided that theassumptions below hold:
(i) limm→∞σm=0, ∞∑m=1σm=∞, σm∈(0,1), ρm∈[ℓ1,ℓ2]⊂(0,1).
(ii) ∞∑m=0πm‖νm−νm−1‖<∞.
Proof. Based on the results of Martin [42,43] and Cioranescu [38], Υ is m−accretive. This shows that JΥ=(I+Υ)−1 is nonexpansive and ∇(JΥ)=N(Υ). Putting JΥ=Z in Theorem 3.1 and continuing with the same approach, we get the desired result.
Remark 4.2. ● The problem of finding zeros of accretive mappings in a real uniformly convex Banach space with uniformly Gâteaux differentiable norm given in (4.5) above gave us the motivation to extend the result of Tan and Cho [37] from Hilbert spaces to real uniformly convex Banach spaces with uniformly Gâteaux differentiable norm.
● If we set πm=0 in our algorithm (3.1), then, we have Ishikawa iterative scheme [45]. So, our results extend comparable results for approximating fixed point of nonexapnsive mappings, like the results of Tan and Xu [46]. Moreover, the obtained results here complement the results of Aoyama et al. [47], Chapter 16 of Chidume [11] and Theorem 5.4 of Berinde [10].
Now, we study the behavior of Algorithm (3.1) for approximating the fixed point by the following two experiments:
Example 5.1. Assume that Ω=R with the usual norm. Define a mapping Z:Ω→Ω by
Z(ν)=(5ν2−2ν+48)13,∀ν∈A, |
where the set A is defined by A={ν:0≤ν≤50}.
Experiment 1: In this experiment we have use different values of control parameter σm=1(km+2) for k=1,2,3,5,10. Also, consider σm=1(km+2), ρm=0.80, ν0=ν1=10, Dn=‖νn+1−νn‖,, we have Figures 1 and 2.
Experiment 2: In this experiment we have use different values of control parameter ρm=k for k=0.15,0.35,0.55,0.75,0.95.
Also, consider σm=1(2m+2), ρm=k, ν0=ν1=10, Dn=‖νn+1−νn‖, we have Figures 3 and 4.
Krasnoselskii-Mann iterative scheme is widely used in the solution of the fixed point equation which takes the shape Zx=x, where Z:℧→℧ is nonexpansive mapping and ℧ is a non-empty, closed and convex, subset of a Banach space Ω. This algorithm converges weakly to the fixed point of Z provided the underlying space Ω is a Hilbert space. It is interesting to address the apparent deficiency of the previous algorithm by building an algorithm that converges strongly to the fixed point of Z. For this purpose, in this manuscript, we introduce an inertial Krasnoselskii-Mann Algorithm (3.1) for nonexpansive mappings in a real uniformly convex Banach space with uniformly Gâteaux differentiable norm and prove that the proposed Algorithm (3.1) has strong convergence. Moreover, it should be noted that the evidence here differs from previous literature.
The authors express many thanks to the Editor-in-Chief, handling editor, and the reviewers for their outstanding comments that improve our paper.
The authors declare that they have no competing interests concerning the publication of this article.
[1] |
N. Amenta, M. Godwin, N. Postarnakevich, K. S. John, Approximating geodesic tree distance, Inform. Process. Lett., 103 (2007), 61-65. doi: 10.1016/j.ipl.2007.02.008
![]() |
[2] |
B. R. Baum, Combining trees as a way of combining data sets for phylogenetic inference, and the desirability of combining gene trees, Taxon, 41 (1992), 3-10. doi: 10.2307/1222480
![]() |
[3] |
M. S. Bayzid, T. Warnow, Naive binning improves phylogenomic analyses, Bioinformatics, 29 (2013), 2277-2284. doi: 10.1093/bioinformatics/btt394
![]() |
[4] |
L. J. Billera, S. P. Holmes, K. Vogtmann, Geometry of the space of phylogenetic trees, Adv. Appl. Math., 27 (2001), 733-767. doi: 10.1006/aama.2001.0759
![]() |
[5] |
D. Bryant, R. Bouckaert, J. Felsenstein, N. A. Rosenberg, A. RoyChoudhury, Inferring species trees directly from biallelic genetic markers: bypassing gene trees in a full coalescent analysis, Mol. Biol. Evol., 29 (2012), 1917-1932. doi: 10.1093/molbev/mss086
![]() |
[6] | J. Chakerian, S. Holmes, DISTORY: Distance between phylogenetic histories. R package version, 1 (2013). |
[7] |
J. Chifman, L. Kubatko, Quartet inference from SNP data under the coalescent model, Bioinformatics, 30 (2014), 3317-3324. doi: 10.1093/bioinformatics/btu530
![]() |
[8] | J. A. Cotton, M. Wilkinson, Majority-rule supertrees, Syst. biol., 56 (2007), 445-452. |
[9] | V. Dinh, L. S. T. Ho, M. A. Suchard, F. A. Matsen IV, Consistency and convergence rate of phylogenetic inference via regularization, Ann. Stat., 46 (2018), 1481. |
[10] |
J. Gatesy, M. S. Springer, Phylogenetic analysis at deep timescales: unreliable gene trees, bypassed hidden support, and the coalescence/concatalescence conundrum, Mol. Phylogenet. Evol., 80 (2014), 231-266. doi: 10.1016/j.ympev.2014.08.013
![]() |
[11] | J. Heled, A. J. Drummond, Bayesian inference of species trees from multilocus data, Mol. Biol. Evol., 27 (2009), 570-580. |
[12] |
W. Hoeffding, Probability inequalities for sums of bounded random variables, J. Am. Stat. Assoc., 58 (1963), 13-30. doi: 10.1080/01621459.1963.10500830
![]() |
[13] | S. Ji, J. Kollár, B. Shiffman, A global Łojasiewicz inequality for algebraic varieties, T. Am. Math. Soc., 329 (1992), 813-818. |
[14] |
L. S. Kubatko, B. C. Carstens, L. L. Knowles, STEM: species tree estimation using maximum likelihood for gene trees under coalescence, Bioinformatics, 25 (2009), 971-973. doi: 10.1093/bioinformatics/btp079
![]() |
[15] | M. K. Kuhner, J. Felsenstein, A simulation comparison of phylogeny algorithms under equal and unequal evolutionary rates, Mol. Biol. Evol., 11 (1994), 459-468. |
[16] |
B. R. Larget, S. K. Kotha, C. N. Dewey, C. Ané, BUCKy: gene tree/species tree reconciliation with bayesian concordance analysis, Bioinformatics, 26 (2010), 2910-2911. doi: 10.1093/bioinformatics/btq539
![]() |
[17] |
L. Liu, L. Yu, Estimating species trees from unrooted gene trees, Syst. Biol., 60 (2011), 661-667. doi: 10.1093/sysbio/syr027
![]() |
[18] |
L. Liu, L. Yu, S. V. Edwards, A maximum pseudo-likelihood approach for estimating species trees under the coalescent model, BMC Evol. Biol., 10 (2010), 302. doi: 10.1186/1471-2148-10-302
![]() |
[19] |
S. Mirarab, M. S. Bayzid, B. Boussau, T. Warnow, Statistical binning enables an accurate coalescent-based estimation of the avian tree, Science, 346 (2014), 1250463. doi: 10.1126/science.1250463
![]() |
[20] |
S. Mirarab, R. Reaz, M. S. Bayzid, T. Zimmermann, M. S. Swenson, T. Warnow, ASTRAL: genome-scale coalescent-based species tree estimation, Bioinformatics, 30 (2014), i541-i548. doi: 10.1093/bioinformatics/btu462
![]() |
[21] | E. Mossel, S. Roch, Incomplete lineage sorting: consistent phylogeny estimation from multiple loci, IEEE ACM T. Comput. Bi., 7 (2008), 166-171. |
[22] | S. Patel, R. T. Kimball, E. L. Braun, Error in phylogenetic estimation for bushes in the tree of life, Journal of Phylogenetics & Evolutionary Biology, (2013). |
[23] |
D. F. Robinson, Comparison of labeled trees with valency three, J. Comb. Theory B, 11 (1971), 105-119. doi: 10.1016/0095-8956(71)90020-7
![]() |
[24] |
S. Roch, M. Nute, T. Warnow, Long-branch attraction in species tree estimation: inconsistency of partitioned likelihood and topology-based summary methods, Syst. biol., 68 (2019), 281-297. doi: 10.1093/sysbio/syy061
![]() |
[25] |
S. Roch, T. Warnow, On the robustness to gene tree estimation error (or lack thereof) of coalescent-based species tree methods, Syst. Biol., 64 (2015), 663-676. doi: 10.1093/sysbio/syv016
![]() |
[26] |
A. Rokas, B. L. Williams, N. King, S. B. Carroll, Genome-scale approaches to resolving incongruence in molecular phylogenies, Nature, 425 (2003), 798-804. doi: 10.1038/nature02053
![]() |
[27] |
K. P. Schliep, phangorn: phylogenetic analysis in r, Bioinformatics, 27 (2011), 592-593. doi: 10.1093/bioinformatics/btq706
![]() |
[28] |
M. Steel, A. Rodrigo, Maximum likelihood supertrees, Syst. Biol., 57 (2008), 243-250. doi: 10.1080/10635150802033014
![]() |
[29] |
P. Vachaspati, T. Warnow, ASTRID: accurate species trees from internode distances, BMC genomics, 16 (2015), 1-13. doi: 10.1186/1471-2164-16-1
![]() |
1. | Oussama Melkemi, Global Regularity for the Inviscid Micropolar-Rayleigh-Bénard Convection System With Nonlinear Thermal Diffusion, 2025, 0921-7134, 10.1177/09217134241308453 |