
Survival of living tumor cells underlies many influences such as nutrient saturation, oxygen level, drug concentrations or mechanical forces. Data-supported mathematical modeling can be a powerful tool to get a better understanding of cell behavior in different settings. However, under consideration of numerous environmental factors mathematical modeling can get challenging. We present an approach to model the separate influences of each environmental quantity on the cells in a collective manner by introducing the "environmental stress level". It is an immeasurable auxiliary variable, which quantifies to what extent viable cells would get in a stressed state, if exposed to certain conditions. A high stress level can inhibit cell growth, promote cell death and influence cell movement. As a proof of concept, we compare two systems of ordinary differential equations, which model tumor cell dynamics under various nutrient saturations respectively with and without considering an environmental stress level. Particle-based Bayesian inversion methods are used to quantify uncertainties and calibrate unknown model parameters with time resolved measurements of in vitro populations of liver cancer cells. The calibration results of both models are compared and the quality of fit is quantified. While predictions of both models show good agreement with the data, there is indication that the model considering the stress level yields a better fitting. The proposed modeling approach offers a flexible and extendable framework for considering systems with additional environmental factors affecting the cell dynamics.
Citation: Sabrina Schönfeld, Alican Ozkan, Laura Scarabosio, Marissa Nichole Rylander, Christina Kuttler. Environmental stress level to model tumor cell growth and survival[J]. Mathematical Biosciences and Engineering, 2022, 19(6): 5509-5545. doi: 10.3934/mbe.2022258
[1] | Said Attaoui, Billal Bentata, Salim Bouzebda, Ali Laksaci . The strong consistency and asymptotic normality of the kernel estimator type in functional single index model in presence of censored data. AIMS Mathematics, 2024, 9(3): 7340-7371. doi: 10.3934/math.2024356 |
[2] | Yue Li, Yunyan Wang . Strong consistency of the nonparametric kernel estimator of the transition density for the second-order diffusion process. AIMS Mathematics, 2024, 9(7): 19015-19030. doi: 10.3934/math.2024925 |
[3] | Salim Bouzebda, Amel Nezzal . Uniform in number of neighbors consistency and weak convergence of kNN empirical conditional processes and kNN conditional U-processes involving functional mixing data. AIMS Mathematics, 2024, 9(2): 4427-4550. doi: 10.3934/math.2024218 |
[4] | Breix Michael Agua, Salim Bouzebda . Single index regression for locally stationary functional time series. AIMS Mathematics, 2024, 9(12): 36202-36258. doi: 10.3934/math.20241719 |
[5] | Fatimah Alshahrani, Wahiba Bouabsa, Ibrahim M. Almanjahie, Mohammed Kadi Attouch . Robust kernel regression function with uncertain scale parameter for high dimensional ergodic data using k-nearest neighbor estimation. AIMS Mathematics, 2023, 8(6): 13000-13023. doi: 10.3934/math.2023655 |
[6] | M. R. Irshad, S. Aswathy, R. Maya, Amer I. Al-Omari, Ghadah Alomani . A flexible model for bounded data with bathtub shaped hazard rate function and applications. AIMS Mathematics, 2024, 9(9): 24810-24831. doi: 10.3934/math.20241208 |
[7] | Salim Bouzebda . Weak convergence of the conditional single index U-statistics for locally stationary functional time series. AIMS Mathematics, 2024, 9(6): 14807-14898. doi: 10.3934/math.2024720 |
[8] | Gaosheng Liu, Yang Bai . Statistical inference in functional semiparametric spatial autoregressive model. AIMS Mathematics, 2021, 6(10): 10890-10906. doi: 10.3934/math.2021633 |
[9] | Qiang Zhao, Zhaodi Wang, Jingjing Wu, Xiuli Wang . Weighted expectile average estimation based on CBPS with responses missing at random. AIMS Mathematics, 2024, 9(8): 23088-23099. doi: 10.3934/math.20241122 |
[10] | Oussama Bouanani, Salim Bouzebda . Limit theorems for local polynomial estimation of regression for functional dependent data. AIMS Mathematics, 2024, 9(9): 23651-23691. doi: 10.3934/math.20241150 |
Survival of living tumor cells underlies many influences such as nutrient saturation, oxygen level, drug concentrations or mechanical forces. Data-supported mathematical modeling can be a powerful tool to get a better understanding of cell behavior in different settings. However, under consideration of numerous environmental factors mathematical modeling can get challenging. We present an approach to model the separate influences of each environmental quantity on the cells in a collective manner by introducing the "environmental stress level". It is an immeasurable auxiliary variable, which quantifies to what extent viable cells would get in a stressed state, if exposed to certain conditions. A high stress level can inhibit cell growth, promote cell death and influence cell movement. As a proof of concept, we compare two systems of ordinary differential equations, which model tumor cell dynamics under various nutrient saturations respectively with and without considering an environmental stress level. Particle-based Bayesian inversion methods are used to quantify uncertainties and calibrate unknown model parameters with time resolved measurements of in vitro populations of liver cancer cells. The calibration results of both models are compared and the quality of fit is quantified. While predictions of both models show good agreement with the data, there is indication that the model considering the stress level yields a better fitting. The proposed modeling approach offers a flexible and extendable framework for considering systems with additional environmental factors affecting the cell dynamics.
Numerous fields of study, including econometrics, epidemiology, environmental science, image analysis, oceanography, meteorology, geostatistics, etc., frequently create spatial data. Typically, these data are gathered in numerous fields and analyzed statistically at measurement sites. Consult [1,2,3,4,5], as well as the references contained within, to identify credible sources of references to the research literature in this field and to learn about certain statistical applications. We highlight that modeling a spatio-temporal interaction of financial data is crucial in financial risk management. The momentousness of this matter is motivated by the digitalization of financial transactions. Actually, with technological development, most financial transactions are carried out by the internet allowing for an increase in the impact of the spatial interaction between financial institutes in financial movements. On the other hand with these new instruments in financial operations, the spatial correlation of the financial time series data is not standard. It is related to many additional factors than the geographical position. Indeed, the spatio-temporal feature in financial time is impacted by the economic exchange between countries, the relationship between economic sectors, and the mobility between countries, among others. Thus, analyzing the spatio-temporal features of financial data is tough work. Motivated by this challenge, we combine the ideas of functional single index modeling with the recent financial risk techniques to provide a statistical model that allows us to fit the spatio-temporal feature of financial data in risk management. Noting that it is well-recognized that spatio-temporal modeling is a particular case of spatio-functional data analysis, the single index model is more appropriate in econometrics and financial areas. Thus, it will be very interesting to utilize recent developments in spatio-functional statistics to introduce a new financial risk model based on a single index structure. Indeed, the FSIR model is one of the key tools in econometrics as well as in financial time series data. In particular, for financial areas, this model is often used to reduce the high number of factors of an investment or to determine the principal assets in a given portfolio. From a theoretical point of view, the single index model belongs to the semiparametric family of models, which many authors have studied from practical and theoretical point of view, for instance, see [6,7]. We return to [8,9] for the first results in the vectorial explanatory case. In the functional setting, the authors of [10] propose the Nadaraya-Watson-kernel-estimator (NWKE) for the nonparametric part in a functional single index regression (FSIR) structure. They proved the consistency of the NWKE when the covariate pertains to the Hilbertian subspace. In the last few decades, the popularity of these models has increased. We cite, for instance, [11] who introduced a new estimation in single index modeling. They proposed a multi-index fitting aproach adaptable to linear projections for functional data and show that their method makes it possible to predict with polynomial convergence rates. Recently, [12] generalized the functional parametric regression model to partially linear functional single index models. More recent advances in functional single index models were obtained by [13]. They used the k-nearest neighbors algorithm to estimate the nonparametric link function of the FSIR. They elaborated the Borel-Cantelli consistency (BCC) of the NWKE using a quasi-associated dependence structure.
In this paper, we investigate conditional exectiles, which is based on least asymmetrically weighted squares estimation, which was adopted from the econometrics literature and is a fundamental statistical application tool. This method frequently employs the [14] concept of expectiles, the least-squares equivalent of the conventional quantiles. They were given this name because they resemble the quantiles of a random variable, but, unlike quantiles, they are based on a quadratic loss function, as in the case of the expectation; see [15,16] for more information. Since it is the only elicitable coherent risk measure, we refer to [17] and its references. Refer to [18] for applying the expectile regression in heteroscedasticity analysis. We refer you to the recent paper by [19,20] for further justification of the expectile model's application. For an overview of the use of expectile curves in regression analysis, see [21] and the extensive discussions of that paper, especially the research [22] for an evaluation of expectiles and [23] for a critical perspective. Despite their disparities in construction, quantiles and expectiles share similar characteristics. As demonstrated by the research of [24], the primary reason is that expectiles are identical to quantiles if the original distribution is transformed. Quantiles and expectiles, which comprise information about a random variable's complete distribution, are extensions of the median and mean, respectively. Expectiles are superior replacements for quantiles in a variety of pertinent applications. Motivating advantages include the fact that expectiles are more sensitive than quantiles to the magnitude of infrequent catastrophic losses and that they depend on both the tail realizations of the predictor and their probability, whereas quantiles depend only on the frequency of tail realizations. This sensitivity of expectiles to tail behavior enables more prudent and responsive risk management. Observe that the quantiles are not always adequate and can be criticized for being challenging to compute because the corresponding loss function is not continuously differentiable. Of course, these features increase the importance of the expectile in some specific areas, namely, in financial risk analysis. However, the robustness of the quantile is also a good advantage in some alternative areas, namely for the prediction issue. From a statistical point of view, there is a bothersome trade-off between sensitivity and robustness. Thus, it will be interesting to develop a bridge between the two phenomenons. This is the main motivation of the expectile regression as an alternative model to the quantile when the sensitivity is required. The functional expectile with regression was recently introduced by [20,25]. They demonstrated almost complete consistency with the rate and the limit distribution of NWKE of this functional model. Such results were obtained under the independence structure. The last contribution was extended to the mixing situation by [19]. The authors of [26] stated the BCC of the functional NWKE when the functional time is ergodic. The spatial structure was handled by [27]. It should be noted that this last situation is of great importance in practice. The framework of the present contribution concerns the nonparametric spatial statistic. The pioneer works of this domains are the reaearches [15,16,28,29,30,31], among other. This subject was investigated for the first time in functional statistics by [32]. They proved the BCC functional NWKE is the regression operator.
The third axe of this contribution concerns the spatio-functional data analysis (SFDA). Such a topic is relatively recent. It is classified by [33] as second-generation functional data analysis. Of course, the main difficulty of this topic in the analysis of spatial data comes from the fact that observations indexed in the multi-dimensional space do not have a linear order, as it happens for classical time series data. For a brief literature review in SFDA, we cite the first pioneer work [34]. They consider the problem of the functional linear regression estimation by the spline method when the observations are spatially correlated. Among the wide range of applications of SFDA we mention the hyperspatial imagery processing developed by [35]. We return to [36] for more advanced application of SFDA in the environmental area. The mathematical development of nonparametric spatio-functional data analysis was started by [37]. They state the asymptotic property of the kernel estimation of the density of a functional random field. We refer to [38] for spatial local linear estimation in functional data. Recently, the authors of [39] treated spatio-functional quantile regression using different techniques, both parametric and nonparametric. However, despite the great importance of the semi-parametric modeling of spatio-functional data in practice, this problem has not been explored yet. To the best of our knowledge, the present contribution is the first work in this direction. It is worth noting that, the functional data analysis (FDA) is actually in continuous progression. For an overview of the recent developments and trends in this area, we may refer the reader to some journal special issues devoted to this topic, such as [40,41,42,43,44,45] to cite a few.
Precisely, the main aim of the present work is to estimate the spatial expectile regression when the input and the output variables are linked with the FSIR structure. We construct an estimator of the nonparametric part of this link function using kernel weighting. The consistency with the rate of the constructed estimator is the main asymptotic result of this work. Moreover, the main novelty of this contribution is the treatment of this model under the FSIR structure. From a practical point of view, this kind of semiparametric model is very common in the econometric area. Such popularity is motivated by two important features. The first is its characteristic as an excellent reducer of the data dimension, and the second feature is the easy interpretation of the co-variability between the output variable and the regressor of the functional index included in these models. The estimability of the functional index is discussed using two cross-validation rules. The first is based on the weighted least squared error, and the second is obtained by the maximum likelihood function. A simulation investigation is conducted to compare these criteria. Finally, we highlight the great impact of the present contribution in financial time series data via a real-data application, allowing us to show the superiority of our spatial model over its competitors and to explore the spatial interaction in financial data. Recall that the flexibility of the additive model and the importance of the spatial interaction of the financial activities are the principal motivations to investigate the expectile by regression in spatio-functional SIM. Thus, our study allows us to digitalize the spatial co-movement in financial time series data using new algorithms adapted to high-frequency data observed over a thinner discretization grid. The new method proposed in this contribution constitutes an alternative approach to classical models based on quantile by regression and shortfall regression.
This paper is organized as follows: In the following sections, we introduce the model and the kernel estimator of the conditional expectile function. The main asymptotic results are stated in Section 3. Section 4.2 is devoted to the functional index θ's estimability issue. The applicability of considered rules will be examined in Section 5. The last section also includes a real world data application. Some concluding remarks are given in Section 6. Finally, the proofs of the main results are given in Appendix.
Set N∈IN∗ and assume (Xi,Yi), i∈INN as a F×IR-valued strictly stationary spatial process. (Xi,Yi), i∈INN is defined in (Ω,A,Pr), a given probability space. Specifically, F is a separable Hilbert space where ⟨⋅,⋅⟩ is the inner product and with an orthonormal basis denoted by {ep:p≥1}. We suppose that the process (Xi,Yi) is observed under a rectangular area
In={i=(i1,⋯,iN)∈INN:1≤ik=1,⋯,N≤nk=1,⋯,N}. |
We will write that the N-uplet in INN,
n=(n1,⋯,nN)→∞, |
if min{nk}→∞ and |nj/nk|<C for a constant 0<C<∞, and for all j, k is less than N. Typically, we explore the spatial interaction of
Zi=(Xi,Yi)i∈In∈N |
by assuming that there exists a non-increasing function φ(t) towards 0 as t→∞, such that, all subsets E′ and E of INN have finite cardinals
α(B(E′),B(E))=supB∈B(E′),C∈B(E)|Pr(B)Pr(C)−Pr(B∩C)|≤ψ(Card(E′),Card(E))φ(dist(E′,E)), | (2.1) |
where Card(E′) and Card(E) are the cardinality of E′ and E, B(E′) (resp. B(E)) is the Borel σ-algebra generated from {Zi: i∈E′} (resp. {Zi: i∈E}). The quantity dist(E′,E) represents the Euclidean metric between E′ and E. The function
ψ:IN2→IR+ |
is a positive symmetric function and nondecreasing such that, for all m′,m∈IN,
ψ(m′,m)≤Cmin(m′,m), | (2.2) |
for some C>0. Additionally, the function φ satisfies
∞∑i=1iδφ(i)<∞ for some δ>0. | (2.3) |
All these assumptions that characterize the spatio-functional framework of this study are standards. They are similar to those used by [16,46]. It should be pointed out that if N=1, then (Xi,Yi) is called strongly mixing (see [47] for discussion on mixing and examples). Furthermore, the spatial linear process [48] satisfies this mixing assumption (see [49], for more details on the required conditions). We refer to [1,50,51,52,53] for additional information on mixing coefficients for random fields.
Now, we assume that the behavior of Yi is linked to Xi through the FSIR model with functional θ, a fixed index in F. Therefore, Yi and Xi are linked by
IE[Yi∣Xi]=IE[Yi∣⟨θ,Xi⟩],θin F. | (2.4) |
The model identifiability has been examined by [10]. A sufficient condition of the FSIR-identifiability is:
(i) The link function r(x)=IE[Yi∣Xi=x] is differentiable.
(ii) The functional index θ is such that ⟨θ,e1⟩=1.
The goal of this paper is to estimate the qth expectile with regression of Y given X=z defined, for 0<q<1, by
ξq(θ,z)=argmint∈RΓq(Y,θ,z,t), |
where
Γq(Y,θ,z,t)=IE[q(Y−t)21(Y−t)>0∣⟨z,θ⟩]+IE[(1−q)(Y−t)21(Y−t)≤0∣⟨z,θ⟩]. | (2.5) |
It is worth noticing that (2.5) generalizes the conditional expectation of Y given ⟨z,θ⟩, which coincides with ξq(θ,z), specifically when q=1/2 and such as the conditional quantile generalizes the conditional median. On the other hand, (2.5) is similar to the conditional q-quantile of Y given X=x, which can be obtained by replacing (Y−t)2 by |Y−t| in (2.5). Hence, the name conditional q-expectile. By straightforward calculus, we demonstrate that ξq(θ,z) is a zero of
q1−q=G1(θ,z,t)G2(θ,z,t), |
where
G1(θ,z,t)=−IE[(Y−t)1(Y−t)≤0∣⟨z,θ⟩],G2(θ,z,t)=IE[(Y−t)1(Y−t)>0∣⟨z,θ⟩]. |
Since the function
G(θ,z,t):=G1(θ,z,t)G2(θ,z,t) |
is a non-decreasing function, see for example [20], the expectile ξq(θ,z) of order q is expressed as
ξq(θ,z)=inf{t∈R:G(θ,z,t)≥q1−q}. | (2.6) |
Finally, the spatial estimator of the qth expectile with the regression of Y given X is
^ξq(θ,z)=inf{t∈R/ˆG(θ,z,t)≥q1−q}, | (2.7) |
where
ˆG(θ,z,t)=^G1(θ,z,t)^G2(θ,z,t) |
with
^G1(θ,z,t)=∑i∈InK(h−1nINθ(x−Xi))1(Yi−t)≤0(Yi−t)∑i∈InK(h−1nINθ(x−Xi)) |
and
^G2(θ,z,t)=∑i∈InK(h−1nNθ(x−Xi))1(Yi−t)>0(Yi−t)∑i∈InK(h−1nNθ(x−Xi)), |
where
Nθ(z)=⟨z,θ⟩, |
K(⋅) is a kernel function, and hn is a sequence of positive real numbers tending to zero as n tends to infinity.
Recall that our aim is the establishment of the BCC of the estimator ^ξq(θ,z) to ξq(θ,z) under the strong mixing structure (2.1). To accomplish this goal, we define
B(θ,z,r):={x′∈F/|Nθ(x′−z)|≤r}forr>0. |
Next, we introduce the following conditions:
(H1) ∀r>0,Pr(X∈B(θ,z,r))=:ϕ(θ,z,r)>0. Additionally, ϕ(θ,z,r)⟶0 as r⟶0.
(H2) The random field (Xi,Yi)i∈N satisfies
{IE[|Yi|p|Nθ(Xi)]<C<∞, for some p>4,For all i≠j,IE[|YjYi||Nθ(Xi),Nθ(Xj)]≤C<∞,and 0<supi≠jPr[(Xi,Xj)∈B(θ,z,h)×B(θ,z,h)]≤C(ϕ(θ,z,hn))(a+1)/a,for some0<a<δN−1. |
(H3) The kernel function K(⋅) is supported in [0,1] and there exist C and C′>0, such that
C1(0,1)(⋅)≤K(⋅)≤C′1(0,1)(⋅). |
(H4) The functions Gl(z,⋅,⋅), for l=1,2, continuously-differentiable in IR, for all (t1,t2)∈IR and X1,X2∈F, we have
|Gl(θ,X1,t1)−Gl(θ,X2,t2)|≤C(‖X1−X2‖kl+|t1−t2|ςl), |
for some ςl,kl>0.
(H5) There exists (γn), a sequence of nonnegative real numbers, such that
{γ−1nϕ(1−β)/β(θ,z,hn)→0 with β=p−2p,∑nˆn((1+2ς1)/2ς1)γ−pn<∞,∑nˆn((1+2ς1)/2ς1)−δ/2Nγδ/Nnϕ−δ/2N(θ,z,hn)logδ/2Nˆn<∞, |
where
ˆn=N∏knk. |
The assumptions are standard. Specifically, hypothesis (H1) is stated for several functional processes. For the same special case, see [54] or [55]. These works evaluated (H1) over some Gaussian processes. Condition (H2) is added to get the same convergence rate as in the i.i.d. setting. The assumption (H3) it also standard in this context of functional statistics. In particular, its technical assumption is verified by numerous standard kernels, such as Epanechnikov kernel, beta kernel, triangular kernel, and others (see for the same example [56]). (H4) is a moderate regularity postulate imposed to state the bias term. The requirements in (H5) are reasonable technical conditions to simplify the proofs.
The BCC of ˆξq(z), when (Xi,Yi) satisfies (2.2) and (2.3), is given in the following theorem.
Theorem 1. Under hypotheses (H1)–(H5) and in addition, if we have
∂G(θ,z,ξq(z))∂t>0, |
then,
ˆξq(z)−ξq(z)=O(hk1n)+O(hk2n)+O(logˆnˆnϕ(θ,z,hn))1/2 |
in BCC-consistency-mode, as n→∞.
We point out that the structure of the obtained convergence rate keeps its usual form of the functional kernel smoothing approach in the sense that it is decomposed into its principal terms. The first one is the bias term expressed with respect to the degree of the smoothing assumption in (H4). The second term is the stochastic term, which is expressed with respect to the functional structure through the function ϕ(⋅) of the assumption (H1).
Proposition 3.1. Under conditions (H1)–(H5), we have
supt∈[ξq(z)−ϵ0,ξq(z)+ϵ0]|ˆG(θ,z,t)−G(θ,z,t)|=O(hk1n)+O(hk2n)+O(logˆnˆnϕ(θ,z,hn))1/2 |
in BCC-consistency-mode, as n→∞.
It is clear that the convergence rate of the spatial estimator ˆξq is comparable to the i.i.d case. This statement follows from the observation that, in our theoretical study, we seek to determine the appropriate conditions to get a good asymptotic property of the estimator, namely those reduce the convergence rate of the estimator. However, in practice, such an optimal situation is not usually available. Thus, it is very interesting to examine the effect of the spatial correlation of the data on the computationablity of the constructed estimator. Indeed, in the nonfunctional case, the spatial correlation is evaluated through the covariogram or the variogram function; see [57]. Alternatively, in functional statistics, we use the trace-variogram function (see [58]) to examine the spatial correlation. Specifically, we adopt the ideas of [32] by adding the spatial controller to the definition of the estimator. Indeed, we compute ^ξq(θ,Xk), for a new observation Xk in new site k∉In by replacing ^G1 and ^G2 in (2.7) with
^G1(θ,z,t)=∑i∈InK(h−1nINθ(x−Xi))1Wk(i)1(Yi−t)≤0(Yi−t)∑i∈InK(h−1nINθ(x−Xi))1Wk(i) |
and
^G2(θ,z,t)=∑i∈InK(h−1nNθ(x−Xi))1Wk(i)1(Yi−t)>0(Yi−t)∑i∈InK(h−1nNθ(x−Xi))1Wk(i), |
where Wk is a vicinity set of the fixed site k defined by
Wk={i,such that γ(i,k)≤ιn}, | (4.1) |
where γ is the trace-variogram function and ιn is a appropriate sequence of positive real numbers. We point out that the trace-variogram function γ is estimated empirically by
ˆγ(l,k)=12#Nl,k∑i,j∈Nl,kd(Xi,Xj) |
with
Nl,k={i,j∈Insuch that ‖i−j‖=‖l−k‖}, |
and #Nl,k is the cardinal of Nl,k. Observe that the use of the trace-variogram function also allows for the integration of the functional nature of the data. Furthermore, in the isotropic case where the dependence is related only to the distance between the locations, we can proceed with the vicinity set
Vk={i,such that distance(i,k)≤νn}, | (4.2) |
where νn is an appropriate sequence of positive real numbers. Of course, the distance here is the locating function between the different sites defined by the user.
This section is devoted to showing how we implement ˆξq(z) in practice. Naturally, the applicability of ˆξq(z) depends on the precision of the parameters utilized in the estimator. In this paragraph, we focus on the principal one, which is the single index θ. It is worth noting that the FSIR estimation has been developed by multiple authors, for instance, see [59,60]. However, in this study we will introduce and control the spatial structure of the data. Precisely, we will control this aspect over the three usual selector procedures mentioned below.
The single index model is widely employed in econometrics. It is usually used to reduce the number of factors in econometric data analysis. In mathematical statistics, this kind of model belongs to the family of additive models. Theoretically, these techniques are used to improve the convergence rate of the nonparametric approach. We integrate this vicinity subset of the previous paragraph in the least squares rule to select the best index as
ˆθ=argminθ∈Θ⊂F∑k∈In(Yk−ˆR−kθ(Xk))2, | (4.3) |
where ˆR−kθ is the leave-one-out estimator of the conditional expectation, which is defined by
ˆR−kθ(Xk)=∑j∈In−{k}K(h−1nNθ((Xk)−Xj))Yj1Wk(j)∑j∈In−{k}K(h−1nNθ((Xk)−Xj))1Wk(j). |
As discussed in the second section, this rule is justified by the fact that conditional expectation can be viewed as a particular case of the q-expectile with q=0.5. However, the rule (4.3) can be generalized for various orders q by taking
ˆθ=argminθ∈Θ⊂F∑k∈Inρq(Yk−ˆξ−kq(Xk)), | (4.4) |
where
ρq(s)=|q−1{s<0}|s2, |
and ˆξ−kq is the leave one-out estimator of ξq constructed in the same manner as ˆR−kθ.
An alternative approach based on the maximum likelihood method is used to select the best optimal index. Indeed, with this rule, the best index model is realized by maximizing the conditional likelihood function
θ=argmaxθ∈Θ⊂Ff(y∣Nθ(X)), |
where f(⋅∣⋅) is the density of Y conditioning on Nθ(X). So, the practical determination of the single index is found on the nonparametric estimation of f(⋅∣⋅). Once again, to explore the spatial correlation, we integrate the same vicinity subset in the conditional density function
ˆf(y|Xk)=h−1n∑i∈InK(h−1nNθ(x−Xi))K(h−1n(y−Yi))1Wk(i)∑i∈InK(h−1nNθ(x−Xi))1Wk(i). |
Therefore, we have
ˆθ=argmaxθ∈Θ⊂F1nn∑i=1logˆf(Yi∣ˆr−iθ,n(Xi)). | (4.5) |
This criterion was used by [61] in the nonfunctional case. The use of conditional density outperforms the first criterion (4.3) based on the conditional expectation. Since the conditional density is more informative than the conditional expectation, (4.5) seems to be more adequate than (4.3).
The aim of this subsection is the evaluation of the effectiveness of the proposed estimator using a finite sample size. The goal is to show how the spatial interaction impacts the choice the functional-index θ as well as the smoothing parameter hn. Particularly, we will check the effect of spatial dependency over the two cross-validation rules of the previous sections. For this purpose, we simulate spatio-functional data using the SFIM as follows:
Yi=rθ(Xi)+ϵi |
with
rθ(⋅)=r(Nθ(⋅)). |
The function r(⋅) denotes a nonparametric regression link and ϵi is a white noise spatial process that is supposed to be a Gaussian isotropic random field. The covariance function of this spatial process is
C(u)=e(−u22). |
For this experimental analysis, we generate the response variable by taking
r(z)=∫101x2(t)+1 |
and θ=e1 is the first element of the basis function of Karhunen-Loève decomposition. Thereby, θ is the eigenfunction that corresponds to the greatest eigenvalue of the covariance of the process (Xi)i. The latter is drawn from the formula
Xi=sin(πWit)+Witcos(Wit),t∈[0,1], |
where Wi is a random Weibull field with covariance that has an exponential function
C(u)=e(−ψu) for u≥0. |
The simulation result is displayed in Figure 1.
Recall that θ is unknown in practice. To estimate θ, we compare the two rules of the previous sections. Based on the two rules, we select the best functional index from Θ the finite subset defined by
Θ=Θn={θ∈F,θ=k0∑i=1ciei,‖θ‖=1,and ∃j∈1,…,k such that Nθ(ej)>0}, |
where (ci)i are some calibrated real constants that ensure model identifiability. Usually, we choose the (ci)i from {−1,0,1} with calibration. Next, for this computational study, we assume that the functional subset Θ belongs to the Hilbert subspace spanned by the finite basis functions of (ei)i=1,⋯,k0. This basis function constitutes the k0-eigenfunction associated with the k0 largest eigenvalue. For the sake of brevity, we have fixed k0=5. Thereafter, we use the same cross-validation rules to select the smoothing parameter hn. Specifically, we have
hn=argminargminh∈Hn∑k∈Inρq(Yk−ˆξ−kq(Xk)), |
where
Hn={a≥0:∑i∈In1B(z,a)(Xi)=k}, |
k∈{5,15,25,⋯,0.5ˆn}, and the ball is defined with respect to the L2-distance between the functional regressors. We simulate with quadratic kernel-defined as
K(t)=32(1−t2)1[0,1). |
The effectiveness of the estimator ˆξq is evaluated by computing the mean square errors (MSE)
MSE(p)=1ˆn∑j(ξq(Xj)−ˆξq(Xj))2, | (5.1) |
where the theoretical expectile regression ξq is obtained by using the routine code qenorm in the R-package VGAM.
The results are recorded in Table 1. It compares the SCVLSE-rule and SCVML-rule with an arbitrary selector method. The latter of these rules is obtained by dividing the optimal parameters of the SCVLSE-rule by 2. We emphasize our choice of utilizing the neighborhood set Wk, which enables the integration of the spatial component into the functional component through the trace-variogram function. The sequence ιn in the vicinity set Wk is selected among the quantile of order q of the estimator vector trace-variogram γ(k,i) given by the code trace.variog in the geofd package. The following table contains the MSE(q) for q=0.01, q=0.05, q=0.1, ˆn=200 and ψ=0.5,2, which describes the covariance function of the random Weibull field in the regressor. Unsurprisingly, the efficiency of ˆξq is heavily affected by the spatial correlation as well as the selection rule to choose the parameter involved in the computation of ˆξq, such as the functional index θ and the smoothing parameters hn. It is clear that the arbitrary way significantly destroys the estimation quality. On the other hand, the effect of the spatial correlation of the data also impacts the estimation quality. Indeed, it is clear that the estimation quality decreases with large values of ψ. To better illustrate this observation, we plot in Figure 2 theMSE value with respect to the values of ψ.
rule | value of ψ | q=0.01 | q=0.05 | q=0.1 |
SCVLSE-rule | 0.5 | 0.34 | 0.37 | 0.31 |
2 | 0.19 | 0.15 | 0.22 | |
SCVML-rule | 0.5 | 0.28 | 0.32 | 0.34 |
2 | 0.11 | 0.09 | 0.08 | |
Arbitrary -rule | 0.5 | 0.77 | 0.73 | 0.84 |
2 | 0.61 | 0.68 | 0.62 |
In the second illustration, we compare our approach to its competitive such as the parametric and the nonparametric methods. To conduct a fair comparison between the three algorithms, we regenerate the response variable using three different situations. Indeed, in addition to the initial data of the first illustration, we define
Yi=∫θ(t)(Xi)(t)dt+ϵi, linear case,Yi=r(Xi)+ϵi,purely nonparametric case, |
where θ, r, ϵi, and Xi are defined in the first setting. We were inspired by [27] to define the spatial functional linear expectile regression as
¯ξq(Xk)=⟨¨ξq,Xk⟩with¨ξq=M∑j=1^sj,q^υj, |
where the vector
(^s1,q^s2,q⋮^sM,q)=argminζ∈RM∑i∈Inρq(Yi−m∑j=1ζj⟨^υj,Xi⟩), |
and where the (ζj)j are the components of ζ in the basis function, and (^υj)j=1,⋯,M, for the M eigenfunctions associated with the M greatest eigenvalues of the spatial empirical version of the covariance operator
Γn(u)=1ˆn∑i∈In⟨Xi,u⟩Xi. |
For the spatial nonparametric functional expectile regression we adopt the estimator of [27], which is defined by estimating G1 and G2 as
~G1(θ,Xk,t)=∑i∈InK(h−1n‖Xk−Xi‖)1Wk(i)1(Yi−t)≤0(Yi−t)∑i∈InK(h−1n‖Xk−Xi‖)1Wk(i) |
and
~G2(θ,Xk,t)=∑i∈InK(h−1n‖Xk−Xi‖)1Wk(i)1(Yi−t)>0(Yi−t)∑i∈InK(h−1n‖Xk−Xi‖)1Wk(i). |
Using the same procedure as in the first setting, we select the sequence of the vicinity set (ιn), the smoothing parameter, the same metric, and the same kernel. We examine the performance of the three approaches ¯ξq, ~ξq and ^ξq using the MSE of (5.1) presented in Figures 3–5. Unsurprisingly, the semi-parametric estimator ^ξq(Xk) is more stable for the three situations. In the sense that its MSE has slow variability with respect the different situations. It is of order 0.33 in the linear case, 0.36 in the semi-parametric case and 0.38 of the purely nonparametric case. The MSEs of ¯ξq are (0.27, 0.45, 0.88) versus (0.97, 0.67, 0.31) for ~ξq. We observe also that the estimators ¯ξq and ~ξq outperform only for the parametric and nonparametric situations, respectively.
Commonly, the financial area is the natural field of expectile regression. In this area, expectile regression is used as an alternative model of risk to the expected-shortfall and the value-at-risk (VaR). It is more informative than the mentioned risk tools, such as the VaR function. This statement results from the fact that the expectile model is the tail expectation, whereas the VaR function is the quantile model of the tail probabilities. The tail expectation function covers the frequencies as well as the values, whereas the tail probability is based only on the frequency. The novelty of the spatial expectile is the possibility to fit the financial risk of the co-movements of various investments in different sectors or stock markets. In a sense, the spatio-functional correlation of the financial data regroups the time conventional correlation as well as the pairwise relationships between the stock markets through known financial metrics. It is worth noting the fact that the spatial linkage in financial data is unrelated to the geographic localization of the stock markets. Such spatial financial distances are usually deafened by some spatial matrix weighting. Thus, we seek in this computational part to inspect the behavior of the spatio-functional conditional expectile concerning some common spatial matrix weighting. For this aim, we consider the Euro Stoxx-50 index data. Such spatio-functional observations are available at https://fred.stlouisfed.org/series (accessed on 14 March 2023). We proceed with the difference logarithmic of this data for the closed prices r(⋅) of the period between 22 February 2022 to 23 February 2023. Specifically, we build a functional variable from
Z(t)=−100log(r(t)r(t−1)) |
as a continuous process. The real interest variable Y is Z (of the-last-day of month), and the functional insert variable X(⋅) illustrates the values of Z (for one month). In Figure 6, we plot the observed functional regressors.
Of course, the first step of spatial modeling is spatial detrending, which is necessary to use the stationarity assumption. To do that, we generate a stationary process (˜Yl,˜Xl)l from the initial spatial observation (Yl,Xi)l as
{˜Xl=Xl−m1(l),˜Yl=Yl−m2(l). |
Next, we compute the conditional expectile estimator ˆξq from the statistics (ˆXi,ˆYi)i instead of from the initial observations (Xi,Yi)i. Thus, the used observations are obtained by estimating the real functions m1(⋅) and m2(⋅) by
ˆm1(i0)=∑i∈InWi0iXi |
and
ˆm2(j0)=∑j∈InW′j0jYj, |
where Wij and W′ij are given spatial weighting matrices. As mentioned above, we evaluate the impact of this step in the spatio-financial risk by comparing three common weighting matrices:
∙ The first one is
W1ij={1, if i and j are in the same sector,0, if not. |
∙ The second one
W2ij={1,if i and j are in the same sector and the same country,0.5,if i and j are in the same sector,0,if not. |
∙ The third one
W3ij={1,if i=j,0,if not. |
The last matrix allows us to examine also the behavior of the spatio-functional expetile regression without the detrending part. Now, to run our spatio-functional model, we retrain the same schemes as those exercised in the artificial example to designate the parameters of ˆξq. Specifically, we use the quadratic function supported on (0,1) as the kernel and we select the single index and the smoothing parameter by the rule SCVLSE as
(θopt,hn)=argminargminh∈Hn,θinΘn∑k∈Inρq(Yk−ˆξ−kq(Xk)), |
where Θn and Hn are defined in the same manner as in the previous section. However, we use the metric of principal component analysis to specify the ball in Hn. For both subsets Θn and Hn, the principal component analysis is performed over 5-eigenfunctions associated to the 5 largest eigenvalues. The feasibility of the proposed spatial risk analysis is checked for q=0.95 and q=0.05 by dividing the data several times (exactly 55 times). The observations are divided at random into two parts: the training sample (220 observations) and the testing sample (150 observations). Finally, we assess the feasibility of the detrending steps by measuring
Eror=1150150∑i=1ρ0.01(Yi−ˆξ0.01(Xi)), |
where
ρq(s)=|q−1{s<0}|s2. |
The values Eror for the 55 random splitting operations of the sample are plotted in Figure 7.
It is obvious that there is a remarkable difference between the stationary case and the nonstationary situation. In a sense, when we execute with the initial data without detrending (matrix 3), the estimation quality is poor compared to the detrending step (using matrices 1 and 2). In particular, the detrending step permits us to increase the effectiveness of the estimator ˆξq by reducing the risk measure Eror(⋅). The scatter plot in Figure 7 demonstrates that the median of Eror-values remarkably varies between the three matrices. It is around 0.3 for the matrix W1ij, around 0.2 for W2ij and 0.6 in the nonstationary case associated with the matrix W3ij.
In this paper, we have demonstrated the BCC consistency of spatio-functional expectile regression under the FSIR structure. Such a result is established over some general postulates, allowing us to explore the nonparametric nature of the expectile operator, the functional path of the financial time series data, and the spatial correlation of the observed data. On the other hand, since the degree of correlation greatly impacts the convergence speed of the estimator, we have modeled this feature using various rules. First, we have evaluated the effect of spatial dependency on the choice of the single index model. For this purpose, we have employed vicinity-set techniques and the spatial weighting matrix. We observed that the two approaches fit the spatial dependency correctly and are a good tool for controlling the spatial covariation of the financial data. Indeed, since financial transactions are performed via the internet, the spatial dependency between them is not based only on the location of the financial institutions. Thus, the vicinity-set techniques (defined by the trace-variogram function) and the spatial weighting matrix algorithm allows for the integration of all the different elements affecting the spatial correlation. Moreover, we show that the insertion of the expectile operator in financial risk management is carried out via two principal steps: detrending and determination of the estimator. Such a strategy increases the efficiency of our algorithm in practice. Additionally, to this theoretical and practical development, the present contribution opens very interesting tracks for future research. For example, it will be a priority in the future to study the asymptotic property of the parametric estimation of the spatio-functional expectile operator. Such a prospect is motivated by the expectile operator's ability to behave as linear, nonparametric, or semiparametric forms [62,63,64]. Thus, the ideas of [30] can be extended here. Second, the asymptotic normality of the estimators is also crucial in mathematical statistics. It allows us to determine the confidence interval with a given confidence level. Extending our results to other functional time series cases (ergodic, long memory, associated process) would be interesting. However, it would require nontrivial mathematics that is well beyond the scope of this paper.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The authors would like to thank the associate-editor and the four anonymous reviewer for their valuable comments and suggestions which improved substantially the quality of an earlier version of this paper. The authors thank and extend their appreciation to the funders of this project: 1) Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2024R358), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia; 2) The Deanship of Scientific Research at King Khalid University through the Research Groups Program under grant number R.G.P. 1/366/44.
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
This section is devoted to the proof of our main result. The previously presented notation continues to be used in the following.
Proof of Theorem 1. The function ˆG is non-decreasing, and therefore its derivative is non-negative at ξq(θ,z). So, for any ϵ>0, we have
∑nPr(|^ξq(θ,z)−ξq(θ,z)|>ϵ)≤∑nPr(|ˆG(θ,z,ξq(θ,z)−ϵ)−G(θ,z,ξq(θ,z)−ϵ)|≥Cϵ)+∑nPr(|ˆG(θ,z,ξq(θ,z)+ϵ)−G(θ,z,ξq(θ,z)+ϵ)|≥Cϵ). |
Next, it suffices to use the Proposition 3.1, and for
t=ξq(θ,z)±ϵ, |
we obtain
∑nPr(supt∈[ξq(θ,z)−δ,ξq(θ,z)+δ]|ˆG(θ,z,t)−G(θ,z,t)|>Cϵ)<∞. | (A.1) |
So, in order to show the consistency of ^ξq(θ,z), we have to prove the uniform BCC of ˆG1,2(θ,z,t), as n→∞. This requirement is a consequence of Proposition 1. Hence, it suffices to prove this proposition first.
Proof of Proposition 3.1. For this proposition, we decompose ˆG(θ,z,t)−G(θ,z,t) as follows:
ˆG(θ,z,t)−G(θ,z,t)=ˆG1(θ,z,t)ˆG2(θ,z,t)−G1(θ,z,t)G2(θ,z,t)=1ˆG2(θ,z,t)[ˆG1(θ,z,t)−G1(θ,z,t)]+G(θ,z,t)ˆG2(θ,z,t)[G2(θ,z,t)−ˆG2(θ,z,t)]. | (A.2) |
Thus, the proposition, as well as the theorem, are consequences of the following intermediate results.
Lemma 1. Under the conditions of Proposition 3.1, we have
supt∈[ξq(θ,z)−δ,ξq(θ,z)+δ]|ˆG1(θ,z,t)−IE[ˆG1(θ,z,t)]|=Oa.co.(√logˆnˆnϕ(θ,z,hn))asn→∞ | (A.3) |
and
supt∈[ξq(θ,z)−δ,ξq(θ,z)+δ]|ˆG2(θ,z,t)−IE[ˆG2(θ,z,t)]|=Oa.co.(√logˆnˆnϕ(θ,z,hn))asn→∞. |
Proof of Lemma 1. Since the proof of the two terms is the same, we focus on the first term only. To evaluate the maximum of this dispersion term on a compact interval, we recover [ξq(θ,z)−δ,ξq(θ,z)+δ], by finite compact intervals [tj−ℓn,tj+ℓn], with
ℓn=ˆn−1/2ςl |
and
dn=O(ˆn1/2ς1). |
Let
Gn={tj−ℓn,tj+ℓn,1≤j≤dn}, | (A.4) |
the subset of the covering interval's extremities. Now, from the monotonicity of ˆG1(θ,z,⋅) and IE[ˆG1(θ,z,⋅)], we write, for 1≤j≤dn,
ˆG1(θ,z,tj−ℓn)≤supt∈(tj−ℓn,tj+ℓn)ˆG1(θ,z,t)≤ˆG1(θ,z,tj+ℓn),IE[ˆG1(θ,z,tj−ℓn)]≤supt∈(tj−ℓn,tj+ℓn)IE[ˆG1(θ,z,t)]≤IE[ˆG1(θ,z,tj+ℓn)]. |
It follows that
supt∈[ξq(θ,z)−δ,ξq(θ,z)+δ]|IE[ˆG1(θ,z,t)]−ˆG1(θ,z,t)|≤max1≤j≤dnmaxz∈{tj−ℓn,tj+ℓn}|IE[ˆG1(θ,z,t)]−ˆG1(θ,z,z)|+2ς1C2ℓς1n. |
Since
ℓς1n=o((√logˆnˆnϕ(θ,z,hn))1/2), |
we prove only the fact that
max1≤j≤dnmaxz∈{tj−ℓn,tj+ℓn}|IE[ˆG1(θ,z,t)]−ˆG1(θ,z,z)|=O(√logˆnˆnϕ(θ,z,hn)),a.co. | (A.5) |
To do that we write
Pr(max1≤j≤dnmaxz∈{tj−ℓn,tj+ℓn}|IE[ˆG1(θ,z,t)]−ˆG1(θ,z,z)|>η√logˆnˆnϕ(θ,z,hn))≤ 2dnmax1≤j≤dnmaxz∈{tj−ℓn,tj+ℓn}Pr(|IE[ˆG1(θ,z,t)]−ˆG1(θ,z,z)|>η√logˆnˆnϕ(θ,z,hn)). |
So, all that remains is to evaluate for all z=tj∓ℓn,1≤j≤dn,
Pr(|IE[ˆG1(θ,z,t)]−ˆG1(θ,z,z)|>η√logˆnˆnϕ(θ,z,hn)). | (A.6) |
To simplify the notation, we let
Y−i=1(Yi−t)≤0(Yi−t), |
which is not necessarily bounded. Thus, to evaluate (A.6), we use the truncation method. Indeed, we consider
^G∗1(θ,z,t)=1ˆnIE[K1(θ,z)]∑i∈InKi(θ,z)Y∗i, |
where
Y∗i=Y−i1(|Yi|<γn) |
with
Ki(θ,z)=K(h−1nNθ(z−Xi)). |
So, (A.3) is a consequence of the following results:
|IE[^G∗1(θ,z,t)]−IE[ˆG1(θ,z,t)]|=O(√logˆnˆnϕ(θ,z,hn)), | (A.7) |
|^G∗1(θ,z,t)−ˆG1(θ,z,t)|=Oa.co.(√logˆnˆnϕ(θ,z,hn)) | (A.8) |
and
|^G∗1(θ,z,t)−IE[^G∗1(θ,z,t)]|=Oa.co.(√logˆnˆnϕ(θ,z,hn)). | (A.9) |
For statement (A.7): For this equation, we use Holder's inequality (for α=p2 and β such that 1α+1β=1), which allows us to write that
|IE[|^G∗1(θ,z,t)]−IE[|ˆG1(θ,z,t)]|≤C1IE[K1(θ,z)]IE[|Y−|1{Y≥γn}K1(θ,z)]≤γ−αnIE[K1(θ,z)]IE1/α[|Y2α|]IE1/β[Kβ1(θ,z)]≤γ−αnIE[K1(θ,z)]IE1/α[|Yp|]IE1/β[Kβ1(θ,z)]. |
It follows that
|IE[|^G∗1(θ,z,t)]−IE[|ˆG1(θ,z,t)]|≤Cγ−αnϕ(1−β)/βx(θ,z,hn). |
Finally, (H5) allows us to conclude the statement (A.7).
For statement (A.8): The demonstration of this statement is based on the Markov inequality. For all ϵ>0, we have
Pr(|ˆG1(θ,z,t)−^G∗1(θ,z,t)|>ϵ)≤∑i∈InPr(Yi>γn)≤ˆnPr(Y>γn)≤ˆnγ−pnIE[Yp]. |
Choose
ϵ=ϵ0(√logˆnˆnϕ(θ,z,hn)) |
to deduce that
∑nPr(|ˆG1(θ,z,t)−^G∗1(θ,z,t)|>ϵ0(√logˆnˆnϕ(θ,z,hn)))≤∑nˆnγ−pn<∞. |
For statement (A.9): This is based on the blocks of spatial decomposition insights [30]. Specifically, we write
Λi=Ki(θ,z)Y∗i−IE[K1(θ,z)Y∗i], |
which allows us to write
^G∗1(θ,z,t)−IE[^G∗1(θ,z,t)]=1ˆnIE[K1(θ,z)]∑i∈InΛi(θ,z). |
The latter can be decomposed as
^G∗1(θ,z,t)−IE[^G∗1(θ,z,t)]=1ˆnIE[K1(θ,z)]2N∑i=1T(n,i), | (A.10) |
with, for all i∈(1,N) and ri=2nip−1n,
T(n,i)=∑l∈JM(i,n,j), |
where
J={0,⋯,r1−1}×⋯×{0,⋯,rN−1}, |
with
M(1,n,l)=2lkpn+pn∑ik=2lkpn+1k=1,⋯,NΛi(θ,z),M(2,n,l)=2lkpn+pn∑ik=2lkpn+1k=1,⋯,N−1(lN+1)pn∑iN=2lNpn+pn+1Λi(θ,z),M(3,n,l)=2lkpn+pn∑ik=2lkpn+1k=1,⋯,N−22(lN−1+1)pn∑iN−1=2lN−1pn+pn+12lNpn+pn∑iN=2lNpn+1Λi(θ,z),M(4,n,l)=2lkpn∑ik=2lkpn+1k=1,⋯,N−22(lN−1+1)pn∑iN−1=2lN−1pn+pn+12(lN+1)pn∑iN=2lNpn+pn+1Λi(θ,z),⋯ |
where the last one is
M(2N,n,l)=2(lk+1)pn∑ik=2lkpn+pn+1k=1,⋯,NΛi(θ,z). |
We note that, if the ni is not exactly equal 2ripn, we regroup the remaining terms in T(n,2N+1). Clearly, the first term T(n,1) is the leading one and the quantity ^G∗1(θ,z,t) is a finite sum. Thus, from (A.10), for η>0, we get
Pr(|IE[^G∗1(θ,z,t)]−^G∗1(θ,z,t)|≥η)≤2Nmaxi=1,⋯2NPr(T(n,i)≥ηˆnIE[K1(θ,z)]). |
Therefore, (A.8) is a consequence of
Pr(T(n,i)≥ηˆnIE[K1(θ,z)]) |
for all i=1,⋯,2N.
We now treat the leading term T(n,1). Of course, the other cases are proved by the same treatment. The proof of this last case is based on the application of Lemma 2 in [46]. Indeed, we recount the
M=N∏k=1rk=2−Nˆnp−Nn≤ˆnp−Nn |
variables in another arbitrary way, that is Z1(θ,z),⋯,ZM(θ,z). In a sense, for each Zj(θ,z) there exists j in J such that
Zj(θ,z)=∑i∈I(1,n,j)Λi(θ,z), |
where
I(1,n,l)={i:2lkpn+1≤ik≤2lkpn+pn,k=1,⋯,N}. |
Clearly, the sets I(1,n,l) are distanced by pn and they contain pNn sites. In addition, we have
K(h−1nNθ(z−Xi))Y∗i≤Cγn. |
Then, from Lemma 2 in [46], we extract independent random variables Z∗1(θ,z),⋯,Z∗M(θ,z) having the same distribution as
Zl=1,⋯,M(θ,z), |
such that
r∑j=1IE|Zj(θ,z)−Z∗j(θ,z)|≤2CγnMpNnψ(pNn(M−1),pNn)φ(pn). | (A.11) |
So, we write
Pr(T(n,i)≥ηˆnIE[K1(θ,z)])≤N1(n)+N2(n), |
where
N1(n)=Pr(|M∑j=1Z∗j|≥MηˆnIE[K1(θ,z)]2M), |
N2(n)=Pr(M∑j=1|Zj−Z∗j|≥ηˆnIE[K1(θ,z)]2). |
We start by evaluating the term N1. As an independent array, we use the Bernstein's inequality. The latter is based on the variance quantity
Var[Z∗1(θ,z)]=Var[∑i∈I(1,n,1)Λi(θ,z)]. |
For this purpose, we use the fact that
IE[Ypi|Xi]<∞, |
for p>2, to prove that
Var[Λi(θ,z)]≤CIE[K2iY∗2i]≤CIE[K2iY2i]≤CIE[K2iIE[Y2i|Xi]]≤CIE[K2i]≤Cϕ(θ,z,h). |
Therefore,
∑i∈I(1,n,1)Var[Λi(θ,z)]=O(pNnϕ(θ,z,h)). |
Next, we use the second part of (H2) to prove that, for all i≠j,
Cov(Λi(θ,z),Λj(θ,z))≤CIE[KiKj|YiYj|] ≤CIE[KiKjIE[|YiYj||XiXj]]≤CIE[KiKj]≤Cϕ(a+1)/ax(θ,z,h). |
On the other hand, since
IE[Ypi|Xi]<∞, |
then, for all i≠j, Hölder's inequality allows us to write
Cov(Λi(θ,z),Λj(θ,z))≤‖Λi(θ,z)‖2pφ1−2/p(‖j−i‖)≤Cϕ2/px(θ,z,h)φ1−2/p(‖i−j‖)). |
Hence,
∑i≠j∈I(1,n,1)|Cov(Λi(θ,z),Λj(θ,z))|≤∑{i,j∈I(1,n,1)‖i−j‖≤un}|Cov(Λi(θ,z),Λj(θ,z))|+∑{i,j∈I(1,n,1)‖i−j‖>un}|Cov(Λi(θ,z),Λj(θ,z))|≤CpNnϕ(θ,z,h)(uNnϕ(θ,z,h)1/a+u−Nanϕ2/p−1x(θ,z,h)∑i:‖i‖≥un‖i‖Naφ1−2/p(‖i‖)). |
It suffices to choose
un=ϕ(θ,z,h)2/Np(a+1)−1/Na, |
and we obtain
∑i≠j∈I(1,n,1)|Cov(Λi(θ,z),Λj(θ,z))|≤CpNnϕ(θ,z,h). |
So, we infer that
Var[∑i∈I(1,n,1)Λi(θ,z)]=O(pNnϕ(θ,z,h)). |
Now, we replace in the inequality of N1, that is,
N1(n)≤2exp(−(ηˆnIE[K1(θ,z)])2MVar[Z∗1]+CηγnpNnˆnIE[K1(θ,z)]). | (A.12) |
This gives
N1(n)≤exp(−Cη0logˆn). |
Finally, a good choice of η0 allows us to write that
∑nN1(n)<∞. |
For the term N2(n), we use the Markov inequality and (A.11) to get that
N2(n)≤2MγnpNn(ηˆnIE[K1(θ,z)])−1ψ(pNn(M−1),pNn)φ(pn). |
Now, since
{{\mathrm {I\!E}}}\left[K_{\bf 1}(\theta,\mathfrak{z})\right]\leq C \phi(\theta,\mathfrak{z},h_{\bf n}),\; \; { }\widehat{\mathbf{n}} = 2^NMp_{\mathbf{n}}^N |
and
\psi(p_{\mathbf{n}}^N(M-1),p_{\mathbf{n}}^N)\leq p_{\mathbf{n}}^N, |
by choosing
\eta = { } \eta_0\sqrt{\frac{\log \widehat{\mathbf{n}}}{\widehat{\mathbf{n}}\, \phi(\theta,\mathfrak{z},h)}}, |
we readily obtain
\mathfrak{N}_2(\mathbf{n})\leq \widehat{\mathbf{n}} \gamma_{{\bf n}} p_{\mathbf{n}}^N\left( \log \widehat{\mathbf{n}}\right)^{-1/2}\left(\widehat{\mathbf{n}}\phi(\theta,\mathfrak{z},h)\right)^{-1/2}\varphi(p_{\mathbf{n}}). |
It suffices to choose
p_{\mathbf{n}} = C\left({ }\frac{\widehat{\mathbf{n}}\phi(\theta,\mathfrak{z},h)}{\log \widehat{\mathbf{n}}\gamma_{{\bf n}}^2}\right)^{1/2N} |
to get
\begin{equation*} \mathfrak{N}_2 (\mathbf{n})\leq \widehat{\mathbf{n}}\, \varphi(p_{\mathbf{n}}). \end{equation*} |
From (H5), we conclude that
\sum\limits_{\mathbf{n}} \mathfrak{N}_2(\mathbf{n}) < \infty. |
In a similar way, as {\bf n} \rightarrow \infty , we get
\sup\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]} \left| \widehat{G}_{2}(\theta,\mathfrak{z},t)-{{\mathrm {I\!E}}}\left[\widehat{G}_{2}(\theta,\mathfrak{z},t)\right] \right| = O_{a.co.}\Bigg( \sqrt{\frac{\log \widehat{\bf n}}{\widehat{\bf n}\phi(\theta,\mathfrak{z},h_{\bf n})}}\Bigg), |
which completes the demonstration of this lemma.
Lemma 2. Assume that the conditions (H1), (H3), and (H4) are fulfilled. We have, as {\bf n} \rightarrow \infty ,
\begin{equation} \sup\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]} \left|G_{1}(\theta,\mathfrak{z},t)-{{\mathrm {I\!E}}}\left[\widehat{G}_{1}(\theta,\mathfrak{z},t)\right]\right| = O\Bigg(h_{\bf n}^{k_{ 1}}\Bigg) \end{equation} | (A.13) |
and
\sup\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]} \left|G_{2}(\theta,\mathfrak{z},t)-{{\mathrm {I\!E}}}\left[\widehat{G}_{2}(\theta,\mathfrak{z},t)\right]\right| = O\Bigg(h_{\bf n}^{k_{2}}\Bigg). |
Proof of Lemma 2. Similar to the previous lemma, as the proof of two terms is the same, we focus on the term \widehat{G}_{1} . For this goal, we use the stationarity of the observations and write
\begin{eqnarray*} { {{\mathrm {I\!E}}}\left[\widehat{G}_{1}(\theta,\mathfrak{z},t)\right]-G_{1}(\theta,\mathfrak{z},t) } \\& = & \frac{ 1}{ {{\mathrm {I\!E}}}\left[K_{\bf 1}(\theta,\mathfrak{z})\right]} \left\{{{\mathrm {I\!E}}}\left[K_{\bf 1}(\theta,\mathfrak{z})\left(G_{1}(\theta, {{\cal X}_{\bf 1}},t)-G_{1}(\theta,\mathfrak{z},t)\right)\right]\right\}. \end{eqnarray*} |
Then, from (H4), we get
\begin{eqnarray*} {{{\mathrm {I\!E}}}\left[K_{\bf 1}(\theta,\mathfrak{z})\left(G^{1}(\theta,;{\cal X}_{\bf 1},t)-G_{1}(\theta,\mathfrak{z},t)\right)\right]}\\ & = & {{\mathrm {I\!E}}}\left[K_{\bf 1}(\theta,\mathfrak{z})\left(G_{1}(\theta,{{\cal X}_{\bf 1}},t)-G_{1}(\theta,\mathfrak{z},t)\right) \mathbb{1}_{B(\mathfrak{z},h_{n})}({{\cal X}_{\bf 1}})\right]. \end{eqnarray*} |
Therefore, we infer
\begin{eqnarray*} { \Big|{{\mathrm {I\!E}}}\left[\widehat{G}_{1}(\theta,\mathfrak{z},t)\right]-G_{1}(\theta,\mathfrak{z},t)\Big|}\\ & = &\frac{ 1}{ {{\mathrm {I\!E}}}\left[K_{\bf 1}(\theta,\mathfrak{z})\right]} \Big|{{\mathrm {I\!E}}}\left[K_{\bf 1}(\theta,\mathfrak{z})\left(G^{1}(\theta,{\cal X}_{\bf 1},t)-G_{1}(\theta,\mathfrak{z},t)\right) \mathbb{1}_{B(\mathfrak{z},h_{n})}({\cal X}_{\bf 1})\right]\Big|. \end{eqnarray*} |
It follows that
\Big|G_{1}(\theta,\mathfrak{z},t)-{{\mathrm {I\!E}}}\left[\widehat{G}_{1}(\theta,\mathfrak{z},t)\right]\Big| \leq C h_{\bf n}^{k_{ 1}}. |
The last inequality is uniform in t , then write
\sup\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]} \left|G_{1}(\theta,\mathfrak{z},t)-{{\mathrm {I\!E}}}\left[\widehat{G}_{1}(\theta,\mathfrak{z},t)\right]\right| = O\Bigg(h_{\bf n}^{k_{ 1}}\Bigg). |
By the same analytical arguments, we obtain
\sup\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]} \left|G_{2}(\theta,\mathfrak{z},t)-{{\mathrm {I\!E}}}\left[\widehat{G}_{2}(\theta,\mathfrak{z},t)\right]\right| = O\Bigg(h_{\bf n}^{k_{2}}\Bigg). |
Hence, the proof is complete.
Lemma 3. Under the conditions of Proposition 3.1, as {\bf n} \rightarrow \infty , we have
\begin{equation} \sum\limits_{\bf n}Pr\Bigg(\inf\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]}\Big| \widehat{G}_{2}(\theta,\mathfrak{z},t)\Big| \leq \epsilon\Bigg) < \infty. \end{equation} | (A.14) |
Proof of Lemma 3. Note that
\begin{eqnarray*} { \inf\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]}\Big| \widehat{G}_{2}(\theta,\mathfrak{z},t)\Big| \leq \frac{1}{2}\inf\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]}G_2(\theta,\mathfrak{z},t)}\\ &\Longrightarrow& \sup\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]}\Big| \widehat{G}_{2}(\theta,\mathfrak{z},t)-G_2(\theta,\mathfrak{z},t)\Big| > \frac{1}{2}\inf\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]}G_2(\theta,\mathfrak{z},t). \end{eqnarray*} |
This statement means
\begin{eqnarray*} { Pr\left( \inf\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]}\Big| \widehat{G}_{2}(\theta,\mathfrak{z},t)\Big| \leq \frac{1}{2}G_2(\theta,\mathfrak{z},t)\right)} \\ &\leq& Pr\left(\sup\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]}\Big| \widehat{G}_{2}(\theta,\mathfrak{z},t)-G_{2}(\theta,\mathfrak{z},t)\Big| > \frac{1}{2}\inf\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]}G_2(\theta,\mathfrak{z},t)\right). \end{eqnarray*} |
Finally, we combine the Lemmas 1 and 2 and choose
\epsilon = \frac{1}{2}\inf\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]}G_2(t;x) > 0 |
to conclude that
\begin{equation} \begin{array}{ccc} \sum\limits_{\bf n}Pr\left(\inf\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta,\, \xi_q(\theta,\mathfrak{z})+\delta]}\Big| \widehat{G}_{2}(\theta,\mathfrak{z},t)\Big| \leq \frac{1}{2}\inf\limits_{t\in [\xi_q(\theta,\mathfrak{z})-\delta, \xi_q(\theta,\mathfrak{z})+\delta ]}G_2(\theta,\mathfrak{z},t\right)& < & \infty. \end{array} \end{equation} | (A.15) |
Hence, the proof is complete.
[1] | L. Preziosi (ed.), Cancer Modelling and Simulation, 1^{st} edition, CRC Press, 2003. |
[2] |
H. Byrne, T. Alarcon, M. Owen, S. Webb, P. Maini, Modelling aspects of cancer dynamics: A review, Philosoph. Transact. Ser. A Math. Phys. Eng. Sci., 364 (2006), 1563–1578. https://doi.org/10.1098/rsta.2006.1786 doi: 10.1098/rsta.2006.1786
![]() |
[3] |
K. Brindle, New approaches for imaging tumour responses to treatment, Nat. Rev. Cancer, 8 (2008), 94–104. https://doi.org/10.1038/nrc2289 doi: 10.1038/nrc2289
![]() |
[4] |
D. A. Hormuth, A. M. Jarrett, E. A. B. F. Lima, M. T. McKenna, D. T. Fuentes, T. E. Yankeelov Mechanism-based modeling of tumor growth and treatment response constrained by multiparametric imaging data, JCO Clin. Cancer Inform., 3 (2019), 1–10. https://doi.org/10.1200/CCI.18.00055 doi: 10.1200/CCI.18.00055
![]() |
[5] |
T. E. Yankeelov, R. G. Abramson, C. C. Quarles, Quantitative multimodality imaging in cancer research and therapy, Nat. Rev. Clin. Oncol., 11 (2014), 670–680. https://doi.org/10.1038/nrclinonc.2014.134 doi: 10.1038/nrclinonc.2014.134
![]() |
[6] |
C. S. Szot, C. F. Buchanan, J. W. Freeman, M. N. Rylander, 3D in vitro bioengineered tumors based on collagen I hydrogels, Biomaterials, 32 (2011), 7905–7912. https://doi.org/10.1016/j.biomaterials.2011.04.001 doi: 10.1016/j.biomaterials.2011.04.001
![]() |
[7] |
E. A. B. F. Lima, N. Ghousifam, A. Ozkan, J. T. Oden, A. Shahmoradi, M. N. Rylander, et al., Calibration of multi-parameter models of avascular tumor growth using time resolved microscopy data, Sci. Rep., 8 (2018), 1–14. https://doi.org/10.1038/s41598-018-32347-9 doi: 10.1038/s41598-018-32347-9
![]() |
[8] | C. P. Robert, G. Casella, G. Casella, Monte Carlo statistical methods, 2^{nd} edition, Springer-Verlag, New York, 2004. https://doi.org/10.1007/978-1-4757-4145-2 |
[9] |
A. Beskos, A. Jasra, E. A. Muzaffer, A. M. Stuart, Sequential Monte Carlo methods for Bayesian elliptic inverse problems, Stat. Comput., 25 (2015), 727–734. https://doi.org/10.1007/s11222-015-9556-7 doi: 10.1007/s11222-015-9556-7
![]() |
[10] |
Y. Zhou, A. M. Johansen, J. A. Aston, Toward automatic model comparison: An adaptive sequential Monte Carlo approach, J. Comput. Graph. Statist., 25 (2016), 701–726. https://doi.org/10.1080/10618600.2015.1060885 doi: 10.1080/10618600.2015.1060885
![]() |
[11] |
R. M. Neal, Annealed importance sampling, Stat. Comput., 11 (2001), 125–139. https://doi.org/10.1023/A:1008923215028 doi: 10.1023/A:1008923215028
![]() |
[12] |
N. Chopin, A sequential particle filter method for static models, Biometrika, 89 (2002), 539–552. https://doi.org/10.1093/biomet/89.3.539 doi: 10.1093/biomet/89.3.539
![]() |
[13] |
P. Del Moral, A. Doucet, A.Jasra, Sequential monte carlo samplers, J. Royal Stat. Soc. Ser. B Stat. Methodol., 68 (2006), 411–436. https://doi.org/10.1111/j.1467-9868.2006.00553.x doi: 10.1111/j.1467-9868.2006.00553.x
![]() |
[14] |
F. Liang, W. H. Wong, Real-parameter evolutionary Monte Carlo with applications to Bayesian mixture models, J. Am. Stat. Assoc., 96 (2001), 653–666. https://doi.org/10.1198/016214501753168325 doi: 10.1198/016214501753168325
![]() |
[15] |
A. Jasra, D. A. Stephens, C. C. Holmes, On population-based simulation for static inference, Stat. Comput., 17 (2007), 263–279. https://doi.org/10.1007/s11222-007-9028-9 doi: 10.1007/s11222-007-9028-9
![]() |
[16] |
G. Evensen, The ensemble Kalman filter: Theoretical formulation and practical implementation, Ocean Dynam., 53 (2003), 343–364. https://doi.org/10.1007/s10236-003-0036-9 doi: 10.1007/s10236-003-0036-9
![]() |
[17] |
C. Schillings, A. M. Stuart, Analysis of the ensemble Kalman filter for inverse problems, SIAM J. Numer. Analys., 55 (2017), 1264–1290. https://doi.org/10.1137/16M105959X doi: 10.1137/16M105959X
![]() |
[18] |
D. Blömker, C. Schillings, P. Wacker, S. Weissmann, Well posedness and convergence analysis of the ensemble Kalman inversion, Inverse Probl., 35 (2019), 085004. https://doi.org/10.1088/1361-6420/ab149c doi: 10.1088/1361-6420/ab149c
![]() |
[19] |
F. J. Richards, A flexible growth function for empirical use, J. Exp. Bot., 10 (1959), 290–301. https://doi.org/10.1093/jxb/10.2.290 doi: 10.1093/jxb/10.2.290
![]() |
[20] |
M. Stoker, H. Rubin, Density dependent inhibition of cell growth in culture, Nature, 215 (1967), 171-–172. https://doi.org/10.1038/215171a0 doi: 10.1038/215171a0
![]() |
[21] |
G. Aubert, J. F. Aujol, A variational approach to removing multiplicative noise, SIAM J. Appl. Math., 68 (2008), 925–946. https://doi.org/10.1137/060671814 doi: 10.1137/060671814
![]() |
[22] |
Y. M. Huang, Yu-Mei, M. K. Ng, and Y. W. Wen, A new total variation method for multiplicative noise removal, SIAM J. Imag. Sci., 2 (2009), 20–40. https://doi.org/10.1137/080712593 doi: 10.1137/080712593
![]() |
[23] |
G. Steidl, T. Teuber, Removing multiplicative noise by Douglas-Rachford splitting methods, J. Math. Imag. Vis., 36 (2010), 168–184. https://doi.org/10.1007/s10851-009-0179-5 doi: 10.1007/s10851-009-0179-5
![]() |
[24] |
R. Šášik, E. Calvo, J. Corbeil, Statistical analysis of high-density oligonucleotide arrays: a multiplicative noise model, Bioinformatics, 18 (2002), 1633–1640. https://doi.org/10.1093/bioinformatics/18.12.1633 doi: 10.1093/bioinformatics/18.12.1633
![]() |
[25] |
A. M. Stuart, Inverse problems: A Bayesian perspective, Acta Numer., 19 (2010), 451–559. https://doi.org/10.1017/S0962492910000061 doi: 10.1017/S0962492910000061
![]() |
[26] |
N. Kantas, A. Beskos, A. Jasra, Sequential Monte Carlo methods for high-dimensional inverse problems: A case study for the Navier–Stokes equations, SIAM/ASA J. Uncertain. Quan., 2 (2014), 464–489. https://doi.org/10.1137/130930364 doi: 10.1137/130930364
![]() |
[27] | P. Del Moral, Feynman-Kac Formulae, Springer-Verlag, New York, 2004. https://doi.org/10.1007/978-1-4684-9393-1 |
[28] | M. Bulté, J. Latz, E. Ullmann, A practical example for the non-linear Bayesian filtering of model parameters, in Quantification of Uncertainty: Improving Efficiency and Technology, Springer, Cham, (2020), 241–272. https://doi.org/10.1007/978-3-030-48721-8_11 |
[29] | H. Haario, E. Saksman, J. Tamminen, An adaptive Metropolis algorithm, Bernoulli, 7 (2001), 223–242. |
[30] | Build models that drive breakthroughs, Technology Networks, 2018. Available from: https://www.technologynetworks.com/cancer-research/ebooks/build-models-that-drive-breakthroughs-311843. |
[31] |
K. Wrzesinski, A. Rogowska-Wrzesinska, R. Kanlaya, K. Borkowski, V. Schwämmle, J. Dai, et. al., The cultural divide: Exponential growth in classical 2D and metabolic equilibrium in 3D environments, PloS One, 9 (2014), e106973. https://doi.org/10.1371/journal.pone.0118050 doi: 10.1371/journal.pone.0118050
![]() |
[32] |
S. Ferson, W. L. Oberkampf, L. Ginzburg, Model validation and predictive capability for the thermal challenge problem, Computer Methods in Applied Mechanics and Engineering, 197 (2008), 2408–2430. https://doi.org/10.1016/j.cma.2004.04.030 doi: 10.1016/j.cma.2004.04.030
![]() |
[33] | R. E. Kass, A. E. Raftery, Bayes factors, J. Am. Statist. Assoc., 90 (1995), 773–795. https://doi.org/10.1080/01621459.1995.10476572 |
[34] | O. J. Schmitz, Ecology and ecosystem conservation. Island Press, 2013. |
[35] | A. Özkan, D. L. Stolley, E. N. Cressman, M. McMillin, S. DeMorrow, T. E. Yankeelov, et. al., Tumor microenvironment alters chemoresistance of hepatocellular carcinoma through CYP3A4 metabolic activity, Front. Oncol., 11 (2021). https://doi.org/10.3389/fonc.2021.662135 |
[36] |
A. Ozkan, N. Ghousifam, P. J. Hoopes, T. E. Yankeelov, M. N. Rylander, In vitro vascularized liver and tumor tissue microenvironments on a chip for dynamic determination of nanoparticle transport and toxicity, Biotechnol. Bioeng., 116 (2019), 1201–1219. https://doi.org/10.1002/bit.26919 doi: 10.1002/bit.26919
![]() |
[37] | A. Özkan, D. L. Stolley, E. N. Cressman, M. McMillin, T. E. Yankeelov, M. N. Rylander, CYP3A4 mediates chemoresistance controlled by cirrhosis and inflammation captured in vascularized hepatocellular carcinoma-on-a-chip, Submitted to Small, (2022). |
[38] |
P. Hartman, A lemma in the theory of structural stability of differential equations, Proceed. Am. Math. Soc., 11 (1960), 610–620. https://doi.org/10.2307/2034720 doi: 10.2307/2034720
![]() |
[39] | D. M. Grobman, Homeomorphism of systems of differential equations (in Russian), Doklady Akademii Nauk SSSR, 128 (1959), 880–881. |
1. | Oussama Bouanani, Salim Bouzebda, Limit theorems for local polynomial estimation of regression for functional dependent data, 2024, 9, 2473-6988, 23651, 10.3934/math.20241150 | |
2. | Salim Bouzebda, Limit Theorems in the Nonparametric Conditional Single-Index U-Processes for Locally Stationary Functional Random Fields under Stochastic Sampling Design, 2024, 12, 2227-7390, 1996, 10.3390/math12131996 | |
3. | Salim Bouzebda, Amel Nezzal, Issam Elhattab, Limit theorems for nonparametric conditional U-statistics smoothed by asymmetric kernels, 2024, 9, 2473-6988, 26195, 10.3934/math.20241280 | |
4. | Breix Michael Agua, Salim Bouzebda, Single index regression for locally stationary functional time series, 2024, 9, 2473-6988, 36202, 10.3934/math.20241719 | |
5. | Fatimah A. Almulhim, Mohammed B. Alamari, Salim Bouzebda, Zoulikha Kaid, Ali Laksaci, Robust Estimation of L1-Modal Regression Under Functional Single-Index Models for Practical Applications, 2025, 13, 2227-7390, 602, 10.3390/math13040602 | |
6. | Changfu Yang, Wenxin Zhou, Wenjun Xiong, Junjian Zhang, Juan Ding, Single-index logistic model for high-dimensional group testing data, 2025, 10, 2473-6988, 3523, 10.3934/math.2025163 | |
7. | Said Attaoui, Oum Elkheir Benouda, Salim Bouzebda, Ali Laksaci, Limit Theorems for Kernel Regression Estimator for Quasi-Associated Functional Censored Time Series Within Single Index Structure, 2025, 13, 2227-7390, 886, 10.3390/math13050886 | |
8. | Ibrahim M. Almanjahie, Hanan Abood, Salim Bouzebda, Fatimah Alshahrani, Ali Laksaci, Nonparametric expectile shortfall regression for functional data, 2025, 58, 2391-4661, 10.1515/dema-2025-0125 |
rule | \mbox{value of } \psi | q=0.01 | q=0.05 | q=0.1 |
SCVLSE-rule | 0.5 | 0.34 | 0.37 | 0.31 |
2 | 0.19 | 0.15 | 0.22 | |
SCVML-rule | 0.5 | 0.28 | 0.32 | 0.34 |
2 | 0.11 | 0.09 | 0.08 | |
Arbitrary -rule | 0.5 | 0.77 | 0.73 | 0.84 |
2 | 0.61 | 0.68 | 0.62 |