Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Functional data analysis: Application to daily observation of COVID-19 prevalence in France

  • Received: 15 November 2021 Revised: 20 December 2021 Accepted: 26 December 2021 Published: 06 January 2022
  • MSC : 62R10, 62P10

  • In this paper we use the technique of functional data analysis to model daily hospitalized, deceased, Intensive Care Unit (ICU) cases and return home patient numbers along the COVID-19 outbreak, considered as functional data across different departments in France while our response variables are numbers of vaccinations, deaths, infected, recovered and tests in France. These sets of data were considered before and after vaccination started in France. After smoothing our data set, analysis based on functional principal components method was performed. Then, a clustering using k-means techniques was done to understand the dynamics of the pandemic in different French departments according to their geographical location on France map. We also performed canonical correlations analysis between variables. Finally, we made some predictions to assess the accuracy of the method using functional linear regression models.

    Citation: Kayode Oshinubi, Firas Ibrahim, Mustapha Rachdi, Jacques Demongeot. Functional data analysis: Application to daily observation of COVID-19 prevalence in France[J]. AIMS Mathematics, 2022, 7(4): 5347-5385. doi: 10.3934/math.2022298

    Related Papers:

    [1] Saudia Jabeen, Bandar Bin-Mohsin, Muhammad Aslam Noor, Khalida Inayat Noor . Inertial projection methods for solving general quasi-variational inequalities. AIMS Mathematics, 2021, 6(2): 1075-1086. doi: 10.3934/math.2021064
    [2] Bancha Panyanak, Chainarong Khunpanuk, Nattawut Pholasa, Nuttapol Pakkaranang . A novel class of forward-backward explicit iterative algorithms using inertial techniques to solve variational inequality problems with quasi-monotone operators. AIMS Mathematics, 2023, 8(4): 9692-9715. doi: 10.3934/math.2023489
    [3] Feng Ma, Bangjie Li, Zeyan Wang, Yaxiong Li, Lefei Pan . A prediction-correction based proximal method for monotone variational inequalities with linear constraints. AIMS Mathematics, 2023, 8(8): 18295-18313. doi: 10.3934/math.2023930
    [4] Doaa Filali, Mohammad Dilshad, Mohammad Akram . Generalized variational inclusion: graph convergence and dynamical system approach. AIMS Mathematics, 2024, 9(9): 24525-24545. doi: 10.3934/math.20241194
    [5] Muhammad Aslam Noor, Khalida Inayat Noor, Bandar B. Mohsen . Some new classes of general quasi variational inequalities. AIMS Mathematics, 2021, 6(6): 6406-6421. doi: 10.3934/math.2021376
    [6] Nagendra Singh, Sunil Kumar Sharma, Akhlad Iqbal, Shahid Ali . On relationships between vector variational inequalities and optimization problems using convexificators on the Hadamard manifold. AIMS Mathematics, 2025, 10(3): 5612-5630. doi: 10.3934/math.2025259
    [7] Xiuzhi Yang, G. Farid, Waqas Nazeer, Muhammad Yussouf, Yu-Ming Chu, Chunfa Dong . Fractional generalized Hadamard and Fejér-Hadamard inequalities for m-convex functions. AIMS Mathematics, 2020, 5(6): 6325-6340. doi: 10.3934/math.2020407
    [8] Shahram Rezapour, Maryam Iqbal, Afshan Batool, Sina Etemad, Thongchai Botmart . A new modified iterative scheme for finding common fixed points in Banach spaces: application in variational inequality problems. AIMS Mathematics, 2023, 8(3): 5980-5997. doi: 10.3934/math.2023301
    [9] Sevda Sezer, Zeynep Eken . The Hermite-Hadamard type inequalities for quasi $ p $-convex functions. AIMS Mathematics, 2023, 8(5): 10435-10452. doi: 10.3934/math.2023529
    [10] Jurancy Ereú, Luz E. Marchan, Liliana Pérez, Henry Rojas . Some results on the space of bounded second κ-variation functions. AIMS Mathematics, 2023, 8(9): 21872-21892. doi: 10.3934/math.20231115
  • In this paper we use the technique of functional data analysis to model daily hospitalized, deceased, Intensive Care Unit (ICU) cases and return home patient numbers along the COVID-19 outbreak, considered as functional data across different departments in France while our response variables are numbers of vaccinations, deaths, infected, recovered and tests in France. These sets of data were considered before and after vaccination started in France. After smoothing our data set, analysis based on functional principal components method was performed. Then, a clustering using k-means techniques was done to understand the dynamics of the pandemic in different French departments according to their geographical location on France map. We also performed canonical correlations analysis between variables. Finally, we made some predictions to assess the accuracy of the method using functional linear regression models.



    In 2014, Li et al. [1] suggested a new class of inverse mixed variational inequality in Hilbert spaces that has simple problem of traffic network equilibrium control, market equilibrium issues as applications in economics and telecommunication network problems. The concept of gap function plays an important role in the development of iterative algorithms, an evaluation of their convergence properties and useful stopping rules for iterative algorithms, see [2,3,4,5]. Error bounds are very important and useful because they provide a measure of the distance between a solution set and a feasible arbitrary point. Solodov [6] developed some merit functions associated with a generalized mixed variational inequality, and used those functions to achieve mixed variational error limits. Aussel et al. [7] introduced a new inverse quasi-variational inequality (IQVI), obtained local (global) error bounds for IQVI in terms of certain gap functions to demonstrate the applicability of IQVI, and provided an example of road pricing problems, also see [8,9]. Sun and Chai [10] introduced regularized gap functions for generalized vector variation inequalities (GVVI) and obtained GVVI error bounds for regularized gap functions. Wu and Huang [11] implemented generalized f-projection operators to deal with mixed variational inequality. Using the generalized f-projection operator, Li and Li [12] investigated a restricted mixed set-valued variational inequality in Hilbert spaces and proposed four merit functions for the restricted mixed set valued variational inequality and obtained error bounds through these functions.

    Our goal in this paper is to present a problem of generalized vector inverse quasi-variational inequality problems. They propose three gap functions, the residual gap function, the regularized gap function, and the global gap function, and obtain error bounds for generalized vector inverse quasi-variational inequality problem using these gap functions and generalized f-projection operator under the monotonicity and Lipschitz continuity of underlying mappings.

    Throughout this article, R+ denotes the set of non-negative real numbers, 0 denotes the origins of all finite dimensional spaces, and , denotes the norms and the inner products in finite dimensional spaces, respectively. Let Ω,F,P:RnRn be the set-valued mappings with nonempty closed convex values, Ni:Rn×RnRn(i=1,2,,m) be the bi-mappings, B:RnRn be the single-valued mappings, and fi:RnR(i=1,2,,m) be real-valued convex functions. We put

    f=(f1,f2,,fm),N(,)=(N1(,),N2(,),,Nm(,)),

    and for any x,wRn,

    N(x,x),w=(N1(x,x),w,N2(x,x),w,,Nn(x,x),w).

    In this paper, we consider the following generalized vector inverse quasi-variational inequality for finding ˉxΩ(ˉx), ˉuF(ˉx) and ˉvP(ˉx) such that

    N(ˉu,ˉv),yB(ˉx)+f(y)f(B(ˉx))intRm+,yΩ(ˉx), (2.1)

    and solution set is denoted by .

    Special cases:

    (i) If P is a zero mapping and N(,)=N(), then (2.1) reduces to the following problem for finding ˉxΩ(ˉx) and ˉuF(ˉx) such that

    N(ˉu),yB(ˉx)+f(y)f(B(ˉx))intRm+,yΩ(ˉx), (2.2)

    studied in[13] and solution set is denoted by 1.

    (ii) If F is single valued mapping, then (2.2) reduces to the following vector inverse mixed quasi-variational inequality for finding ˉxΩ(ˉx) such that

    N(ˉx),yB(ˉx)+f(y)f(B(ˉx))intRm+,yΩ(ˉx), (2.3)

    studied in [14] and solution set is denoted by 2.

    (iii) If CRn is a nonempty closed and convex subset, B(x)=x and Ω(x)=C for all xRn, then (2.3) collapses to the following generalized vector variational inequality for finding ˉxC such that

    N(ˉx),yx+f(y)f(ˉx)intRm+,yC, (2.4)

    which is considered in [10].

    (iv) If f(x)=0 for all xRn, then (2.4) reduces to vector variational inequality for finding ˉxC such that

    N(ˉx),yxintRm+,yC, (2.5)

    studied in [15].

    (v) If Rm+=R+ then (2.5) reduces to variational inequality for finding ˉxC such that

    N(ˉx),yx0,yC, (2.6)

    studied in [16].

    Definition 2.1 [7] Let G:RnRn and g:RnRn be two maps.

    (i) (G,g) is said to be a strongly monotone if there exists a constant μg>0 such that

    G(y)G(x),g(y)g(x)μgyx2,x,yRn;

    (ii) g is said to be Lg-Lipschitz continuous if there exists a constant Lg>0 such that

    g(x)g(y)Lgxy,x,yRn.

    For any fixed γ>0, let G:RnטΩ(,+] be a function defined as follows:

    G(φ,x)=x22φ,x+φ2+2γf(x),φRn,x˜Ω, (2.7)

    where ˜ΩRn is a nonempty closed and convex subset, and f:RnR is convex.

    Definition 2.2 [11] We say that f˜Ω:Rn2˜Ω is a generalized f-projection operator if

    f˜Ωφ={w˜Ω:G(φ,w)=infy˜ΩG(φ,y)},φRn.

    Remark 2.3 If f(x)=0 for all x˜Ω, then the generalized f-projection operator f˜Ω is equivalent to the following metric projection operator:

    P˜Ω(φ)={w˜Ω:wφ=infy˜Ωyφ},φRn.

    Lemma 2.4 [1,11] The following statements hold:

    (i) For any given φRn,f˜Ωφ is nonempty and single-valued;

    (ii) For any given φRn,x=f˜Ωφ if and only if

    xφ,yx+γf(y)γf(x)0,y˜Ω;

    (iii) f˜Ω:RnΩ is nonexpansive, that is,

    f˜Ωxf˜Ωyxy,x,yRn.

    Lemma 2.5 [17] Let m be a positive number, BRn be a nonempty subset such that

    dmfor alldB.

    Let Ω:RnRn be a set-valued mapping such that, for each xRn, Ω(x) is a closed convex set, and let f:RnR be a convex function on Rn. Assume that

    (i) there exists a constant τ>0 such that

    D(Ω(x),Ω(y))τxy,x,yRn,

    where D(,) is a Hausdorff metric defined on Rn;

    (ii) 0wRnΩ(w);

    (iii) f is -Lipschitz continuous on Rn. Then there exists a constant κ=6τ(m+γ) such that

    fΩ(x)zfΩ(x)zκxy,x,yRn,zB.

    Lemma 2.6 A function r:RnR is said to be a gap function for a generalized vector inverse quasi-variational inequality on a set ˜SRn if it satisfies the following properties:

    (i) r(x)0for anyx˜S; (ii) r(ˉx)=0,ˉx˜S if and only if ˉx is a solution of (2.1).

    Definition 2.7 Let B:RnRn be the single-valued mapping and N:Rn×RnRn be a bi-mapping.

    (i) (N,B) is said to be a strongly monotone with respect to the first argument of N and B, if there exists a constant μB>0 such that

    N(y,)N(x,),B(y)B(x)μByx2,x,yRn;

    (ii) (N,B) is said to be a relaxed monotone with respect to the second argument of N and B, if there exists a constant ζB>0 such that

    N(,y)N(,x),B(y)B(x)ζByx2,x,yRn;

    (iii) N is said to be σ-Lipschitz continuous with respect to the first argument with constant σ>0 and -Lipschitz continuous with respect to the second argument with constant >0 such that

    N(x,ˉx)N(y,ˉy)σxy+ˉxˉy,x,ˉx,y,ˉyRn.

    (iv) B is said to be -Lipschitz continuous if there exists a constant >0 such that

    B(x)B(y)xy,x,yRn.

    Example 2.8 The variational inequality (2.6) can be solved by transforming it into an equivalent optimization problem for the so-called merit function r(;τ):X=RnR{+} defined by

    r(x;τ)=sup{N(ˉx),yxXτˉxx2X|xC} for ˉxC,

    where τ is a nonnegative parameter. If X is finite dimensional, the function r(;0) is usually called the gap function for τ=0, and the function r(;τ) for τ>0 is called the regularized gap function.

    Example 2.9 Assume that N:RnRn be a given mapping and C a closed convex set in Rn. Let and be given scalar satisfying >>0 then (2.6) has a D-gap function if

    N(x)=N(x)N(x),xRn

    where D stands for difference.

    In this section, we discuss the residual gap function for generalized vector inverse quasi-variational inequality problem by using the strong monotonicity, relaxed monotonicity, Hausdorff Lipschitz continuity and prove error bounds related to the residual gap function. We define the residual gap function for (2.1) as follows:

    rγ(x)=min1im{B(x)fiΩ(x)[B(x)γNi(u,v)]},xRn,uF(x),vP(x),γ>0. (3.1)

    Theorem 3.1 Suppose that F,P:RnRn are set-valued mappings and Ni:Rn×RnRn(i=1,2,,m) are the bi-mappings. Assume that B:RnRn is single-valued mapping, then for any γ>0,rγ(x) is a gap function for (2.1) on Rn.

    Proof. For any xRn,

    rγ(x)0.

    On the other side, if

    rγ(ˉx)=0,

    then there exists 0i0m such that

    B(ˉx)=fi0Ω(ˉx)[B(ˉx)γNi0(ˉu,ˉv)],ˉuF(ˉx),ˉvP(ˉx).

    From Lemma 2.4, we have

    B(ˉx)[B(ˉx)γNi0(ˉu,ˉv)],yB(ˉx)+γf(y)γf(B(ˉx))0,yΩ(ˉx),ˉuF(ˉx),ˉvP(ˉx)

    and

    Ni0(ˉu,ˉv),yB(ˉx)+f(y)f(B(ˉx))0,yΩ(ˉx),ˉuF(ˉx),ˉvP(ˉx).

    It gives that

    N(ˉu,ˉv),yB(ˉx)+f(y)f(B(ˉx))intRm+,yΩ(ˉx),ˉuF(ˉx),ˉvP(ˉx).

    Thus, ˉx is a solution of (2.1).

    Conversely, if ˉx is a solution of (2.1), there exists 1i0m such that

    Ni0(ˉu,ˉv),yB(ˉx)+fi0(y)fi0(B(ˉx))0,yΩ(ˉx),ˉuF(ˉx),ˉvP(ˉx).

    By using the Lemma 2.4, we have

    B(ˉx)=fi0Ω(ˉx)[B(ˉx)γNi0(ˉu,ˉv)],ˉuF(ˉx),ˉvP(ˉx).

    This means that

    rγ(ˉx)=min1im{B(ˉx)fiΩ(ˉx)[B(ˉx)γNi(ˉu,ˉv)]}=0.

    The proof is completed.

    Next we will give the residual gap function rγ, error bounds for (2.1).

    Theorem 3.2 Let F,P:RnRn be D-ϑF-Lipschitz continuous and D-ϱP-Lipschitz continuous mappings, respectively. Let Ni:Rn×RnRn(i=1,2,,m) be σi-Lipschitz continuous with respect to the first argument and i-Lipschits continuous with respect to the second argument, and B:RnRn be -Lipschitz continuous, and (Ni,B) be strongly monotone with respect to the first argument of Ni and B with positive constant μBi, and relaxed monotone with respect to the second argument of Ni and B with positive constant ζBi. Let

    mi=1(i).

    Assume that there exists κi(0,μBiζBiσiϑF+ϱPi) such that

    fiΩ(x)zfiΩ(y)zκixy,x,yRn,uF(x),vP(x),z{ww=B(x)γNi(u,v)}. (3.2)

    Then, for any xRn and μBi>ζBi+κi(σiϑF+iϱP),

    γ>κiμBiζBiκi(σiϑF+iϱP),
    d(x,)γ(σiϑFiϱP)+γ(μBiζBiκi(σiϑF+iϱP))κirγ(x),

    where

    d(x,)=infˉxxˉx

    denotes the distance between the point x and the solution set .

    Proof. Since

    mi=1(i).

    Let ˉxΩ(ˉx) be the solution of (2.1) and thus for any i{1,,m}, we have

    Ni(ˉu,ˉv),yB(ˉx)+fi(y)fi(B(ˉx))0,yΩ(ˉx),ˉuF(ˉx),ˉvP(ˉx). (3.3)

    From the definition of fiΩ(ˉx)[B(x)γNi(u,v)], and Lemma 2.4, we have

    fiΩ(ˉx)[B(x)γNi(u,v)](B(x)γNi(u,v)),yfiΩ(ˉx)[B(x)γNi(u,v)]+γfi(y)γfi(fiΩ(ˉx)[B(x)γNi(u,v)])0,yΩ(ˉx),uF(x),vP(x). (3.4)

    Since

    ˉxmi=1(i),andB(ˉx)Ω(ˉx).

    Replacing y by B(ˉx) in (3.4), we get

    fiΩ(ˉx)[B(x)γNi(u,v)](B(x)γNi(u,v)),B(ˉx)fiΩ(ˉx)[B(x)γNi(u,v)]+γfi(B(ˉx))γfi(fiΩ(ˉx)[B(x)γNi(u,v)])0,uF(x),vP(x). (3.5)

    Since

    fiΩ(ˉx)[B(x)γNi(u,v)]Ω(ˉx),

    from (3.3), it follows that

    γNi(ˉu,ˉv),fiΩ(ˉx)[B(x)γNi(u,v)]B(ˉx)+γfi(fiΩ(ˉx)[B(x)γNi(u,v)])γfi(B(ˉx))0. (3.6)

    Utilizing (3.5) and (3.6), we have

    γNi(ˉu,ˉv)γNi(u,v)fiΩ(ˉx)[B(x)γNi(u,v)]+B(x),fiΩ(ˉx)[B(x)γNi(u,v)]B(ˉx)0,

    which implies that

    γNi(ˉu,ˉv)γNi(u,v),fiΩ(ˉx)[B(x)γNi(u,v)]B(x)γNi(ˉu,ˉv)γNi(u,v),B(ˉx)B(x)
    +B(x)fiΩ(ˉx)[B(x)γNi(u,v)],fiΩ(ˉx)[B(x)γNi(u,v)]B(x)
    +B(x)fiΩ(ˉx)[B(x)γNi(u,v)],B(x)B(ˉx)0.

    Since F is D-ϑF-Lipschitz continuous, P is D-ϱP-Lipschits continuous and Ni is σi-Lipschitz continuous with respect to the first argument and i-Lipschitz continuous with respect to the second argument, we have

    ˉuuD(F(ˉx),F(x))ϑFˉxx;ˉvvD(P(ˉx),P(x))ϱPˉxx;Ni(ˉu,ˉv)Ni(u,v)σiˉuu+iˉvv. (3.7)

    Again, for i=1,2,,m, (Ni,B) are strongly monotone with respect to the first argument of Ni and B with a positive constant μBi,, and relaxed monotone with respect to the second argument of Ni and B with a positive constant ζBi, we have

    γNi(ˉu,ˉv)γNi(u,v),fiΩ(ˉx)[B(x)γNi(u,v)]B(x)B(x)fiΩ(ˉu)[B(x)γNi(u,v)]2
    +B(x)fiΩ(ˉx)[B(x)γNi(u,v)],B(x)B(ˉx)γμBixˉx2γζBixˉx2.

    By adding fiΩ(x)[B(x)γNi(u,v)] and using the Cauchy-Schwarz inequality along with the triangular inequality, we have

    γNi(ˉu,ˉv)γNi(u,v){fiΩ(ˉx)[B(x)γNi(u,v)]fiΩ(x)[B(x)γNi(u,v)]
    +fiΩ(x)[B(x)γNi(u,v)]B(x)}
    +B(x)B(ˉx){B(x)fiΩ(x)[B(x)γNi(u,v)]+fiΩ(x)[B(x)γNi(u,v)]
    fiΩ(ˉx)[B(x)γNi(u,v)]}γμBixˉx2γζBixˉx2.

    Using the (3.7) and condition (3.2), we have

    (σiϑF+iϱP)γˉxx{κiˉxx+fiΩ(x)[B(x)γNi(u,v)]B(x)}
    +xˉx{B(x)fiΩ(x)[B(x)γNi(u,v)]+κixˉx}γ(μBiζBi)xˉx2.

    Hence, for any xRn and i{1,2,,m}, μBi>ζBi+κi(σiϑF+iϱP),

    γ>κiμBiζBiκi(σiϑF+iϱP),

    we have

    xˉxγ(σiϑF+iϱP)+γ(μBiζBiκi(σiϑF+iϱP))κiB(x)fiΩ(x)[B(x)γNi(u,v)],uF(x),vP(x).

    This implies

    xˉxγ(σiϑF+iϱP)+γ(μBiζBiκi(σiϑF+iϱP))κimin1im{B(x)fiΩ(x)[B(x)γNi(u,v)]}

    which means that

    d(x,)xˉxγ(σiϑF+iϱP)+γ(μBiζBiκi(σiϑF+iϱP))κirγ(x).

    The proof is completed.

    The regularized gap function for (2.1) is defined for all xRn as follows:

    ϕγ(x)=min1imsupyΩ(x),uF(x),vP(x){Ni(u,v),B(x)y+fi(B(x))fi(y)12γB(x)y2}

    where γ>0 is a parameter.

    Lemma 4.1 We have

    ϕγ(x)=min1im{Ni(u,v),Riγ(x)+fi(B(x))fi(B(x)Riγ(x))12γRiγ(x)2}, (4.1)

    where

    Riγ(x)=B(x)fiΩ(x)[B(x)γNi(u,v)],xRn,uF(x),vP(x)

    and if

    xB1(Ω)

    and

    B1(Ω)={ξRn|B(ξ)Ω(ξ)},

    then

    ϕγ(x)12γrγ(x)2. (4.2)

    Proof. For given xRn,uF(x),vP(x) and i{1,2,,m}, set

    ψi(x,y)=Ni(u,v),B(x)y+fi(B(x))fi(y)12γB(x)y2,yRn.

    Consider the following problem:

    gi(x)=maxyΩ(x)ψi(x,y).

    Since ψi(x,) is a strongly concave function and Ω(x) is nonempty closed convex, the above optimization problem has a unique solution zΩ(x). Evoking the condition of optimality at z, we get

    0Ni(u,v)+fi(z)+1γ(zB(x))+NΩ(x)(z),

    where NΩ(x)(z) is the normal cone at z to Ω(x) and fi(z) denotes the subdifferential of fi at z. Therefore,

    z(B(x)γNi(u,v)),yz+γfi(y)γfi(z)0,yΩ(x),uF(x),vP(x)

    and so

    z=fiΩ(x)[B(x)γNi(u,v)],uF(x),vP(x).

    Hence gi(x) can be rewritten as

    gi(x)=Ni(u,v),B(x)fiΩ(x)[B(x)γNi(u,v)]+fi(B(x))fi(fiΩ(x)[B(x)γNi(u,v)])
    12γB(x)fiΩ(x)[B(x)γNi(u,v)]2,uF(x),vP(x).

    Letting

    Riγ(x)=B(x)fiΩ(x)[B(x)γNi(u,v)],uF(x),vP(x),

    we get

    gi(x)=Ni(u,v),Riγ(x)+fi(B(x))fi(B(x)Riγ(x))12γRiγ(x)2,uF(x),vP(x), (4.3)
    (4.4)

    and so

    ϕγ(x)=min1im{Ni(u,v),Riγ(x)+fi(B(x))fi(B(x)Riγ(x))12γRiγ(x)2}.

    From the definition of projection fiΩ(x)[B(x)γNi(u,v)], we have

    fiΩ(x)[B(x)γNi(u,v)]B(x)+γNi(u,v),yfiΩ(x)[B(x)γNi(u,v)]+γfi(y)γfi(fiΩ(x)[B(x)γNi(u,v)])0,uF(x),vP(x). (4.5)

    For any xB1(Ω), we have

    B(x)Ω(x).

    Therefore, putting y=B(x) in (4.5), we get

    γNi(u,v)Riγ(x),Riγ(x)+γfi(B(x))γfi(B(x)Riγ(x))0,uF(x),vP(x),

    that is,

    Ni(u,v),Riγ(x)+fi(B(x))fi(B(x)Riγ(x))1γRiγ(x),Riγ(x)=1γRiγ(x)2. (4.6)

    From the definition of rγ(x) and (4.1), we get

    ϕγ(x)12γrγ(x)2.

    The proof is completed.

    Theorem 4.2 For γ>0,ϕγ is a gap function for (2.1) on the set

    B1(Ω)={ξRn|B(ξ)Ω(ξ)}.

    Proof. From the definition of ϕγ, we have

    ϕγ(x)min1im{Ni(u,v),B(x)y+fi(B(x))fi(y)12γB(x)y2},for allyΩ(x),uF(x),vP(x). (4.7)

    Therefore, for any xB1(Ω), putting y=B(x) in (4.7), we have

    ϕγ(x)0.

    Suppose that ˉxB1(ξ) with ϕγ(ˉx)=0. From (4.2), it follows that

    rγ(ˉx)=0,

    which implies that ˉx is the solution of (2.1).

    Conversely, if ˉx is a solution of (2.1), there exists 1i0m such that

    Ni0(ˉu,ˉv),B(ˉx)y+fi0(B(ˉx))fi0(y)0,yΩ(ˉx),ˉuF(ˉx),ˉvP(ˉx),

    which means that

    min1im{supyΩ(ˉx),ˉuF(ˉx),ˉvP(ˉx){Ni(ˉu,ˉv),B(ˉx)y+fi(B(ˉx))fi(y)12γB(ˉx)y2}}0.

    Thus,

    ϕγ(ˉx)0.

    The preceding claim leads to

    ϕγ(ˉx)0

    and it implies that

    ϕγ(ˉx)=0.

    The proof is completed.

    Since ϕγ can act as a gap function for (2.1), according to Theorem 4.2, investigating the error bound properties that can be obtained with ϕγ is interesting. The following corollary is obtained directly by Theorem 3.2 and (3.5).

    Corollary 4.3 Let F,P:RnRn be D-ϑF-Lipschitz continuous and D-ϱP-Lipschitz continuous mappings, respectively. Let Ni:Rn×RnRn(i=1,2,,m) be σi-Lipschitz continuous with respect to the first argument and i-Lipschitz continuous with respect to the second argument, B:RnRn be -Lipschitz continuous, and (Ni,B) be strongly monotone with respect to the first argument of N and B with respect to the constant μBi>0, and relaxed monotone with respect to the second argument of N and B with respect to the constant ζBi>0. Let

    mi=1(i).

    Assume that there exists κi(0,μBiζBiϑFσi+iϱP) such that

    fiΩ(x)zfiΩ(y)zκixy,x,yRn,uF(x),vP(x)z{ww=B(x)γNi(u,v)}.

    Then, for any xB1(Ω) and any

    γ>κiμBiζBiκi(ϑFσi+iϱP),
    d(x,)γ(ϑFσi+iϱP)+γ(μBiζBiκi(ϑFσi+iϱP))κi2γϕγ(x).

    The regularized gap function ϕγ does not provide global error bounds for (2.1) on Rn. In this section, we first discuss the D-gap function, see [6] for (2.1), which gives Rn the global error bound for (2.1).

    For (2.1) with >>0, the D-gap function is defined as follows:

    G(x)=min1im{supyΩ(x),uF(x),vP(x){Ni(u,v),B(x)y+fi(B(x))fi(y)12B(x)y2}supyΩ(x)uF(x),vP(x){Ni(u,v),B(x)y+fi(B(x))fi(y)12B(x)y2}}.

    From (4.1), we know G can be rewritten as

    G(x)=min1im{Ni(u,v),Ri(x)+fi(B(x))fi(B(x)Ri(x))12Ri(x)2(Ni(u,v),Ri(x)+fi(B(x))fi(B(x)Ri(x))12Ri(x)2)},

    where

    Ri(x)=B(x)fiΩ(x)[B(x)Ni(u,v)]

    and

    Ri(x)=B(x)fiΩ(x)[B(x)Ni(u,v)],xRn,uF(x),vP(x).

    Theorem 5.1 For any xRn,>>0, we have

    12(11)r2(x)G(x)12(11)r2(x). (5.1)

    Proof. From the definition of G(x), it follows that

    G(x)=min1im{Ni(u,v),Ri(x)Ri(x)fi(B(x)Ri(x))12Ri(x)2+fi(B(x)Ri(x))+12Ri(x)2},uF(x),vP(x).

    For any given i{1,2,,m}, we set

    gi(x)=Ni(u,v),Ri(x)Ri(x)fi(B(x)Ri(x))12Ri(x)2+fi(B(x)Ri(x))+12Ri(x)2,uF(x),vP(x). (5.2)

    From fiΩ(x)[B(x)Ni(u,v)]Ω(x), by Lemma 2.4, we know

    fiΩ(x)[B(x)Ni(u,v)](B(x)Ni(u,v)),fiΩ(x)[B(x)Ni(u,v)]fiΩ(x)[B(x)Ni(u,v)]
    +fi(fiΩ(x)[B(x)Ni(u,v)])fi(fiΩ(x)[B(x)Ni(u,v)])0,uF(x),vP(x)

    which means that

    Ni(u,v)Ri(x),Ri(x)Ri(x)+fi(B(x)Ri(x))fi(B(x)Ri(x))0. (5.3)

    Combining (5.2) and (5.3), we get

    gi(x)1Ri(x),Ri(x)Ri(x)12Ri(x)2+12Ri(x)2=12Ri(x)Ri(x)2+12(11)Ri(x)2. (5.4)

    Since

    fiΩ(x)[B(x)Ni(u,v)]Ω(x),

    from Lemma 2.4, we have

    fiΩ(x)[B(x)Ni(u,v)](B(x)Ni(u,v)),fiΩ(x)[B(x)Ni(u,v)]fiΩ(x)[B(x)Ni(u,v)]
    +fi(fiΩ(x)[B(x)Ni(u,v)])fi(fiΩ(x)[B(x)Ni(u,v)])0,uF(x),vP(x).

    Hence

    Ni(u,v)Ri(x),Ri(x)Ri(x)+fi(B(x)Ri(x))
    fi(B(x)Ri(x))0,uF(x),vP(x)

    and so

    1Ri(x),Ri(x)Ri(x)Ni(u,v),Ri(x)Ri(x)fi(B(x)Ri(x))+fi(B(x)Ri(x)).

    It will require and (5.3),

    gi(x)1Ri(x),Ri(x)Ri(x)12Ri(x)2+12Ri(x)2=12Ri(x)Ri(x)2+12(11)Ri(x)2. (5.5)

    From (5.4) and (5.5), for any i{1,2,,m}, we get

    12(11)Ri(x)2gi(x)12(11)Ri(x)2.

    Hence

    12(11)min1im{Ri(x)2}min1im{gi(x)}12(11)min1im{Ri(x)2},

    and so

    12(11)r2(x)G(x)12(11)r2(x).

    The proof is completed.

    Now we are in position to prove that G in the set Rn is a global gap function for (2.1).

    Theorem 5.2 For 0<<, G is a gap function for (2.1) on Rn.

    Proof. From (5.2), we have

    G(x)0,xRn.

    Suppose that ˉxRn with

    G(ˉx)=0,

    then (5.2) implies that

    r(ˉx)=0.

    From Theorem 3.1, we know ˉx is a solution of (2.1).

    Conversely, if ˉx is a solution of (2.1), than from Theorem 3.1, it follows that

    r(ˉx)=0.

    Obviously, (5.2) shows that

    G(ˉx)=0.

    The proof is completed.

    Use Theorem 3.2 and (5.2), we immediately get a global error bound in the set Rn for (2.1).

    Corollary 5.3 Let F,P:RnRn be D-ϑF-Lipschitz continuous and D-ϱP-Lipschitz continuous mappings, respectively. Let Ni:Rn×RnRn(i=1,2,,m) be σi-Lipschitz continuous with respect to the first argument and i-Lipschitz continuous with respect to the second argument, and B:RnRn be -Lipschitz continuous. Let (Ni,B) be the strongly monotone with respect to the first argument of Ni and B with constant μBi and relaxed monotone with respect to the second argument of N and B with modulus ζBi. Let

    mi=1(i).

    Assume that there exists κi(0,μBiζBiϑFσi+iϱP) such that

    fiΩ(x)zfiΩ(y)zκixy,x,yRn,uF(x),vP(x),z{ww=B(x)Ni(u,v)}.

    Then, for any xRn and

    >κiμBiζBiκi(ϑFσi+iϱP),
    d(x,i)(ϑFσi+ϱPi)+(μBiζBiκi(ϑFσi+ϱPi))κi2G(x).

    One of the traditional approaches to evaluating a variational inequality (VI) and its variants is to turn into an analogous optimization problem by notion of a gap function. In addition, gap functions play a pivotal role in deriving the so-called error bounds that provide a measure of the distances between the solution set and feasible arbitrary point. Motivated and inspired by the researches going on in this direction, the main purpose of this paper is to further study the generalized vector inverse quasi-variational inequality problem (1.2) and to obtain error bounds in terms of the residual gap function, the regularized gap function, and the global gap function by utilizing the relaxed monotonicity and Hausdorff Lipschitz continuity. These error bounds provide effective estimated distances between an arbitrary feasible point and the solution set of (1.2).

    The authors are very grateful to the referees for their careful reading, comments and suggestions, which improved the presentation of this article.

    This work was supported by the Scientific Research Fund of Science and Technology Department of Sichuan Provincial (2018JY0340, 2018JY0334) and the Scientific Research Fund of SiChuan Provincial Education Department (16ZA0331).

    The authors declare that they have no competing interests.



    [1] D. Bernoulli, Essai d'une nouvelle analyse de la mortalité causée par la petite vérole, et des avantages de l'inoculation pour la prévenir, Histoire de l'Acad., Roy. Sci. (Paris) avec Mem, 1760, 1–45.
    [2] D. A. Henderson, The eradication of smallpox-An overview of the past, present, and future, Vaccine, 29 (2011), D7–D9. https://doi.org/10.1016/j.vaccine.2011.06.080 doi: 10.1016/j.vaccine.2011.06.080
    [3] D. Wujastyk, Medicine in India, In: J. van Alphen, A. Aris, F. Meyer, M. de Fraeye, Oriental medicine: An illustrated guide to the Asian arts of healing, London: Serindia Publications, 1995, 19–38.
    [4] A. M. Silverstein, A history of immunology, 2 Eds., London: Academic Press, 2009,293.
    [5] L. S. Benjamin, L Melville, Lady Mary Wortley Montagu, her life and letters (1689–1762), Hutchinson, London, 1925.
    [6] R. Ross, An application of the theory of probabilities to the study of a priori pathometry-part I, Proc. R. Soc. Ser. A, 92 (1916), 204–230. https://doi.org/10.1098/rspa.1916.0007 doi: 10.1098/rspa.1916.0007
    [7] A. G. McKendrick, Applications of mathematics to medical problems, Proc. Edinburgh Math. Soc., 44 (1925), 98–130. https://doi.org/10.1017/S0013091500034428 doi: 10.1017/S0013091500034428
    [8] J. Gaudart, O. Touré, N. Dessay, A. L. Dicko, S. Ranque, L. Forest, et al., Modelling malaria incidence with environmental dependency in a locality of Sudanese savannah area, Mali. Malaria J., 8 (2009), 61. https://doi.org/10.1186/1475-2875-8-61 doi: 10.1186/1475-2875-8-61
    [9] J. Gaudart, M. Ghassani, J. Mintsa, M. Rachdi, J. Waku, J. Demongeot, Demography and diffusion in epidemics: Malaria and black death spread, Acta Biotheor., 58 (2010), 277–305. https://doi.org/10.1007/s10441-010-9103-z doi: 10.1007/s10441-010-9103-z
    [10] J. Demongeot, J. Gaudart, A. Lontos, E. Promayon, J. Mintsa, M. Rachdi, Zero-diffusion domains in reaction-diffusion morphogenetic and epidemiologic processes, Int. J. Bifurcation Chaos, 22 (2012), 1250028. https://doi.org/10.1142/S0218127412500289 doi: 10.1142/S0218127412500289
    [11] J. Demongeot, J. Gaudart, J. Mintsa, M. Rachdi, Demography in epidemics modelling, Commun. Pure Appl. Anal., 11 (2012), 61–82. http://dx.doi.org/10.3934/cpaa.2012.11.61 doi: 10.3934/cpaa.2012.11.61
    [12] Z. Liu, P. Magal, O. Seydi, G. Webb, Understanding unreported cases in the COVID-19 epidemic outbreak in Wuhan, China, and importance of major public health interventions, Biology, 9(2020), 50. https://doi.org/10.3390/biology9030050 doi: 10.3390/biology9030050
    [13] J. Demongeot, Q. Griette, P. Magal, SI epidemic model applied to COVID-19 data in mainland China, Royal Soc. Open Sci., 7 (2020), 201878. https://doi.org/10.1098/rsos.201878 doi: 10.1098/rsos.201878
    [14] Z. Liu, P. Magal, O. Seydi, G. Webb, Predicting the cumulative number of cases for the COVID-19 epidemic in China from early data, Math. Biosci. Eng., 17 (2020), 3040–3051. https://doi.org/10.3934/mbe.2020172 doi: 10.3934/mbe.2020172
    [15] P. Magal, O. Seydi, G. Webb, Y. Wu, A model of vaccination for Dengue in the Philippines 2016–2018, Front. Appl. Math. Stat., 7 (2021), 760259. https://doi.org/10.3389/fams.2021.760259 doi: 10.3389/fams.2021.760259
    [16] K. Oshinubi, M. Rachdi, J. Demongeot, Modelling of COVID-19 pandemic vis-à-vis some socioeconomic factors, Front. Appl. Math. Stat., 7 (2021), 786983.
    [17] COVID-19 coronavirus pandemic, 2021. Available from: https://www.worldometers.info/coronavirus.
    [18] Données hospitalières relatives à l'épidémie de COVID-19, 2021. Available from: https://www.data.gouv.fr/fr/datasets/donnees-hospitalieres-relatives-a-lepidemie-de-Covid-19.
    [19] Live COVID-19 vaccination tracker. Available from: https://covidvax.live/location/fra.
    [20] F. Ferraty, P. Vieu, Nonparametric functional data analysis, New York: Springer, 2006. https://doi.org/10.1007/0-387-36620-2
    [21] J. D. Tucker, Functional component analysis and regression using elastic methods, PhD. Thesis, Florida State University, 2014.
    [22] J. D. Tucker, W. Wu, A. Srivastava, Generative models for functional data using phase and amplitude separation, Comput. Stat. Data Anal., 61 (2013), 50–66. https://doi.org/10.1016/j.csda.2012.12.001 doi: 10.1016/j.csda.2012.12.001
    [23] J. O. Ramsay, B. W. Silverman, Applied functional data analysis: Methods and case studies, New York: Springer, 2002. https://doi.org/10.1007/b98886
    [24] A. Srivastava, E. P. Klassen, Functional data and elastic registration, In: Functional and shape data analysis, New York: Springer, 2016, 73–123. https://doi.org/10.1007/978-1-4939-4020-2_4
    [25] J. O. Ramsay, G. Hooker, S. Graves, Functional data analysis with R and MATLAB, New York: Springer, 2009. https://doi.org/10.1007/978-0-387-98185-7
    [26] C. Tang, T. Wang, P. Zhang, Functional data analysis: An application to COVID-19 data in the United States, arXiv. Available from: https://arXiv.org/abs/2009.08363.
    [27] C. Acal, M. Escabias, A. M. Aguilera, M. J. Valderrama, COVID-19 data imputation by multiple function-on-function principal component regression, Mathematics, 9 (2021), 1237. https://doi.org/10.3390/math9111237 doi: 10.3390/math9111237
    [28] T. Boschi, J. Di Iorio, L. Testa, M. A. Cremona, F. Chiaromonte, Functional data analysis characterizes the shapes of the first COVID-19 epidemic wave in Italy, Sci. Rep., 11 (2021), 17054. https://doi.org/10.1038/s41598-021-95866-y doi: 10.1038/s41598-021-95866-y
    [29] Q. Griette, J. Demongeot, P. Magal, A robust phenomenological approach to investigate COVID-19 data for France, Math. Appl. Sci. Eng., 2 (2021), 149–218. https://doi.org/10.5206/mase/14031 doi: 10.5206/mase/14031
    [30] Q. Griette, J. Demongeot, P. Magal, What can we learn from COVID-19 data by using epidemic models with unidentied infectious cases, Math. Biosci. Eng., 2 (2021), 149–160. http://dx.doi.org/10.2139/ssrn.3868852 doi: 10.2139/ssrn.3868852
    [31] J. Gaudart, J. Landier, L. Huiart, E. Legendre, L. Lehot, M. K. Bendiane, et al., Factors associated with spatial heterogeneity of Covid-19 in France: A nationwide ecological study, Lancet Public Health, 6(2021), 222–231. https://doi.org/10.1016/s2468-2667(21)00006-2 doi: 10.1016/S2468-2667(21)00006-2
    [32] O. D. Ilie, R. O. Cojocariu, A. Ciobica, S. I. Timofte, I. Mavroudis, B. Doroftei, Forecasting the spreading of COVID-19 across nine countries from Europe, Asia, and the American continents using the ARIMA models, Microorganisms, 8 (2020), 1158. https://doi.org/10.3390/microorganisms8081158 doi: 10.3390/microorganisms8081158
    [33] J. Stojanovic, V. G. Boucher, J. Boyle, J. Enticott, K. L. Lavoie, S. L. Bacon, COVID-19 is not the flu: Four graphs from four countries, Front. Public Health, 2021, 628479. https://doi.org/10.3389/fpubh.2021.628479 doi: 10.3389/fpubh.2021.628479
    [34] C. Carroll, S. Bhattacharjee, Y. Chen, P. Dubey, J. Fan, A. Gajardo, et al., Time dynamics of COVID-19, Sci. Rep., 10 (2020), 21040. https://doi.org/10.1038/s41598-020-77709-4 doi: 10.1038/s41598-020-77709-4
    [35] A. Srivastava, G. Chowell, Modeling study: Characterizing the spatial heterogeneity of the COVID-19 pandemic through shape analysis of epidemic curves, Res. Square, 2021, 1–27. https://doi.org/10.21203/rs.3.rs-223226/v1 doi: 10.21203/rs.3.rs-223226/v1
    [36] J. Demongeot, Y. Flet-Berliac, H. Seligmann, Temperature decreases spread parameters of the new COVID-19 cases dynamics, Biology, 9 (2020), 94. https://doi.org/10.3390/biology9050094 doi: 10.3390/biology9050094
    [37] H. Seligmann, S. Iggui, M. Rachdi, N. Vuillerme, J. Demongeot, Inverted covariate effects for mutated 2nd vs 1st wave COVID-19: High temperature spread biased for young, Biology, 9 (2020), 226. https://doi.org/10.1101/2020.07.12.20151878 doi: 10.3390/biology9080226
    [38] S. Soubeyrand, J. Demongeot, L. Roques, Towards unified and real-time analyses of outbreaks at country-level during pandemics, One Health, 11 (2020), 100187. https://doi.org/10.1016/j.onehlt.2020.100187 doi: 10.1016/j.onehlt.2020.100187
    [39] J. Demongeot, H. Seligmann, SARS-CoV-2 and miRNA-like inhibition power, Med. Hypotheses, 144 (2020), 110245. https://doi.org/10.1016/j.mehy.2020.110245 doi: 10.1016/j.mehy.2020.110245
    [40] H. Seligmann, N. Vuillerme, J. Demongeot, Unpredictable, counter-intuitive geoclimatic and demographic correlations of COVID-19 spread rates, Biology, 10 (2021), 623. https://doi.org/10.3390/biology10070623 doi: 10.3390/biology10070623
    [41] K. Oshinubi, F. Al-Awadhi, M. Rachdi, J. Demongeot, Data analysis and forecasting of COVID-19 pandemic in Kuwait, MedRxiv, 2021, 1–17. https://doi.org/10.1101/2021.07.24.21261059 doi: 10.1101/2021.07.24.21261059
    [42] J. Demongeot, K. Oshinubi, M. Rachdi, L. Hobbad, M. Alahiane, S. Iggui, et al., The application of ARIMA model to analyze COVID-19 incidence pattern in several countries, J. Math. Comput. Sci., 12 (2022), 1–23. https://doi.org/10.28919/jmcs/6541 doi: 10.28919/jmcs/6541
    [43] K. Oshinubi, M. Rachdi, J. Demongeot, Analysis of reproduction number R0 of COVID-19 using current health expenditure as gross domestic product percentage (CHE/GDP) across countries, Healthcare, 9 (2021), 1247. https://doi.org/10.3390/healthcare9101247 doi: 10.3390/healthcare9101247
    [44] J. Demongeot, K. Oshinubi, M. Rachdi, H. Seligmann, F. Thuderoz, J. Waku, Estimation of daily reproduction rates in COVID-19 outbreak, MedRxiv, 9 (2021), 109. https://doi.org/10.1101/2020.12.30.20249010 doi: 10.1101/2020.12.30.20249010
    [45] J. Demongeot, A. Laksaci, F. Madani, M. Rachdi, Functional data: Local linear estimation of the conditional density and its application, Statistics, 47 (2013), 26–44. https://doi.org/10.1080/02331888.2011.568117 doi: 10.1080/02331888.2011.568117
    [46] M. Rachdi, A. Laksaci, J. Demongeot, A. Abdali, F. Madani, Theoretical and practical aspects on the quadratic error in the local linear estimation of the conditional density for functional data, Comput. Stat. Data Anal., 73 (2014), 53–68. https://doi.org/10.1016/j.csda.2013.11.011 doi: 10.1016/j.csda.2013.11.011
    [47] J. Demongeot, A. Laksaci, M. Rachdi, S. Rahmani, On the local linear modelization of the conditional distribution for functional data, Sankhya A, 76 (2014), 328–355. https://doi.org/10.1007/s13171-013-0050-z doi: 10.1007/s13171-013-0050-z
    [48] J. Demongeot, A. Hamie, A. Laksaci, M. Rachdi, Relative-error prediction in nonparametric functional statistics: Theory and practice, J. Multivar. Anal., 146 (2016), 261–268. https://doi.org/10.1016/j.jmva.2015.09.019 doi: 10.1016/j.jmva.2015.09.019
    [49] J. Demongeot, A. Laksaci, A. Naceri, M. Rachdi, Local linear regression modelization when all variables are curves, Stat. Probab. Lett., 121 (2017), 37–44. https://doi.org/10.1016/j.spl.2016.09.021 doi: 10.1016/j.spl.2016.09.021
    [50] A. Belkis, J. Demongeot, A. Laksaci, M. Rachdi, Functional data analysis: Estimation of the relative error in functional regression under random left-truncation, J. Nonparametr. Stat., 30 (2018), 472–490. https://doi.org/10.1080/10485252.2018.1438609 doi: 10.1080/10485252.2018.1438609
    [51] A. Henien, L. Ait-Hennani, J. Demongeot, A. Laksaci, M. Rachdi, Heteroscedasticity test when the covariables are functionals, C. R. Math., 356 (2018), 571–574. https://doi.org/10.1016/j.crma.2018.02.010 doi: 10.1016/j.crma.2018.02.010
    [52] J. Demongeot, O. Hansen, H. Hessami, A. S. Jannot, J. Mintsa, M. Rachdi, et al., Random modelling of contagious diseases, Acta Biotheor., 61 (2013), 141–172. https://doi.org/10.1007/s10441-013-9176-6 doi: 10.1007/s10441-013-9176-6
    [53] C. J. Rhodes, L. Demetrius, Evolutionary entropy determines invasion success in emergent epidemics, PLoS One, 5 (2010), e12951. https://doi.org/10.1371/journal.pone.0012951 doi: 10.1371/journal.pone.0012951
    [54] S. Triambak, D. P. Mahapatra, A random walk Monte Carlo simulation study of COVID-19-like infection spread, Physica A: Stat. Mech. Appl., 574 (2021), 126014. https://doi.org/10.1016/j.physa.2021.126014 doi: 10.1016/j.physa.2021.126014
    [55] Wikipedia, Available online: https://www.wikipedia.org/wiki/Departments_of_France.
  • This article has been cited by:

    1. Aviv Gibali, , Error bounds and gap functions for various variational type problems, 2021, 115, 1578-7303, 10.1007/s13398-021-01066-8
    2. Soumitra Dey, Simeon Reich, A dynamical system for solving inverse quasi-variational inequalities, 2023, 0233-1934, 1, 10.1080/02331934.2023.2173525
    3. Safeera Batool, Muhammad Aslam Noor, Khalida Inayat Noor, Merit functions for absolute value variational inequalities, 2021, 6, 2473-6988, 12133, 10.3934/math.2021704
    4. Ting Liu, Wensheng Jia, Samuel Nicolay, An Approximation Theorem and Generic Convergence of Solutions of Inverse Quasivariational Inequality Problems, 2022, 2022, 2314-8888, 1, 10.1155/2022/2539961
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(4272) PDF downloads(333) Cited by(21)

Figures and Tables

Figures(16)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog