Next Article in Journal
Novel Algorithms with Inertial Techniques for Solving Constrained Convex Minimization Problems and Applications to Image Inpainting
Previous Article in Journal
Hybrid Learning Moth Search Algorithm for Solving Multidimensional Knapsack Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Adaptive-to-Model Test for Parametric Functional Single-Index Model

1
Faculty of Science, Beijing University of Technology, Beijing 100124, China
2
School of Mathematics and Statistics, Guangxi Normal University, Guilin 541004, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(8), 1812; https://doi.org/10.3390/math11081812
Submission received: 6 March 2023 / Revised: 2 April 2023 / Accepted: 6 April 2023 / Published: 11 April 2023

Abstract

:
Model checking methods based on non-parametric estimation are widely used because of their tractable limiting null distributions and being sensitive to high-frequency oscillation alternative models. However, this kind of test suffers from the curse of dimensionality, resulting in slow convergence, especially for functional data with infinite dimensional features. In this paper, we propose an adaptive-to-model test for a parametric functional single-index model by using the orthogonality of residual and its conditional expectation. The test achieves model adaptation by sufficient dimension reduction which utilizes functional sliced inverse regression. This test procedure can be easily extended to other non-parametric test methods. Under certain conditions, we prove the asymptotic properties of the test statistic under the null hypothesis, fixed alternative hypothesis and local alternative hypothesis. Simulations show that our test has better performance than the method that does not use functional sufficient dimension reduction. An analysis of COVID-19 data verifies our conclusion.

1. Introduction

Consider the functional single-index regression model
Y = g ( β 0 , X ) + ε ,
where Y is a scalar response, X is a predictor function in a functional space H with an inner product denoted by · , · , β 0 H is an unknown functional index, β 0 = 1 , ε is a scalar random noise and E ( ε | X ) = 0 , g ( · ) is a unknown square integrable continuous function. If the model (1) is correctly specified, Refs. [1,2,3,4] proposed consistent estimates of β 0 and g ( · ) in different ways. However, compared to non-parametric regression models without specifying the model form of E ( Y | X ) , parametric functional single-index models are estimated more accurately, especially for functional data with infinite dimension. On the other hand, a wrongly specified model structure would lead to misleading statistical analysis. Therefore, we develop a goodness-of-fit test for model (1) with a given link function g ( · ) . In some cases, the presence of outliers may affect the quality of model fitting, and it is necessary to conduct an outlier test to determine whether outliers exist in the data and then perform the goodness-of-fit test based on the data with outliers removed. The method for the outlier test can be referred to [5], while the primary focus of this paper is on the goodness-of-fit test of model (1).
The special case of the model (1) in which g is linear has been extensively studied. Among them, the research on hypothesis testing has made great progress. For example, testing the nullity of the slope parameter was studied in [6,7,8], and testing whether the conditional expectation of the response given the covariate is almost surely zero or not was considered in [9,10,11]. The model checking method relies on linear model assumptions and is not flexible enough in practical applications. Thus, the proposed model-adaptive approach in this paper overcomes these limitations and offers a more flexible and reliable solution for model checking in parametric functional single-index models. In terms of model testing, the two widely used tests are global smoothing methods based on empirical processes (see [12,13]) and local smoothing methods based on non-parametric estimation (see [9,10,11]). Similar to the conclusion summarized by [14] on vector space, each of these two kinds of tests has advantages and disadvantages. For example, as summarized in [15], the global smoothing methods converge to the limiting distribution at the rate of order O ( n 1 / 2 ) , which is the faster rate that the existing methods can reach. However, they are insensitive to high-frequency oscillation alternative models and often require the use of resampling techniques to obtain critical values. In contrast, as summarized in [14,16], the local smoothing methods have high power in detecting high-frequency oscillation alternative models, but these methods converge to the limiting distribution at the rate of order O ( n 1 / 2 h p / 4 ) , where p = d i m ( X ) . Thus, the local smoothing methods suffer from the curse of dimensionality due to the use of multivariate non-parametric estimation, especially for functional data with infinite dimensional features. To overcome this shortcoming, most of the local smoothing methods mentioned above are based on the semi-metric of the function space to compute univariate non-parametric kernel functions (see [11]), which requires functional principal component analysis for dimensionality reduction. The other method is based on projection technology (see [9,10,17]), which is more complicated for choosing the projection directions. In this paper, we employ sufficient dimension reduction to construct a model adaptive test that achieves a convergence rate of O ( n 1 / 2 h 1 / 4 ) , which is faster compared to that of the locally smooth method when the dimensionality is reduced to 1 under the null hypothesis.
Our aim is to test the parametric model (1) against a parsimonious alternative model class through functional sufficient dimension reduction (FSDR) [18]. We consider the alternative model:
Y = G ( B , X ) + ε ,
where B = ( β 1 , , β K ) is composed of K linearly independent functions in H and column vectors B , X = ( β 1 , X , , β K , X ) is a column vector, G ( · ) is a function from R K R and ε is a scalar error independent of X. The FSDR method is based on the above functional multiple index model. Since the FSDR method can effectively use the information between the functional predictor and the response variable, it is widely used in regression models. In the FSDR literature, the subspace spanned by β 1 , , β K is called the functional sufficient dimension reduction subspace, which is first estimated by functional sliced inverse regression (FSIR) in [19]. Then, functional inverse regression [20], functional K-means inverse regression [21] and functional sliced average variance estimation [18,22] were created. For the test problem studied in this paper, when K > 1 , the model (2) is a functional multiple-index model; when K = 1 and B = β 0 / β 0 , model (2) reduces to the functional single-index model in (1). Therefore, we need construct a test statistic which can adapt to the related model structure through estimating B and its dimension K. In this way, we can achieve that under a null hypothesis, the estimate of B is consistent to β 1 , and under an alternative hypothesis, it is consistent to a multiple-index functional ( β 1 , , β K ) . To achieve this goal, the most critical step is to determine K. We use the modified BIC criterion introduced in [14] to determine K.
This paper is organized as follows. In Section 2, we construct a test statistic based on the orthogonality of residual and its conditional expectation, which is the extension of the model testing method for vector data in [23] to a functional model. In Section 3, we describe in detail the procedure of identifying and estimating B and its dimension K using the FSIR method. The asymptotic properties under null, alternative and local alternative hypotheses are given in Section 4. Simulation results and a real data analysis are reported in Section 5. The corresponding conclusions and discussions are given in Section 6. All proofs are given in Appendix A.

2. The Test Statistic

To verify the validity of model (1), we should develop an appropriate model checking method. Therefore, for a given g ( · ) , we consider the following null hypothesis
H 0 : β 0 H , s u c h t h a t E [ Y | X ] = g β 0 , X ,
and for any K linearly independent function set B = ( β 1 , , β K ) , under model (2), the alternative hypothesis is
H 1 : sup β H P E [ Y | X ] = E [ Y | B , X ] = g β , X < 1 .
Without loss of generality, we take H = L 2 ( [ 0 , 1 ] ) . It is worth noting that under H 0 , K = 1 and B = c 1 β 0 for some constant c 1 ; under H 1 , K 1 . Similar to the interpretation of vector spaces introduced in Stute and Zhu [24], the test statistic based on E [ Y g ( β 0 , X ) | β 0 , X ] is directional, while the test statistic based on E [ Y g ( β 0 , X ) | B , X ] is omnibus. For example, we consider the alternative model: Y = β 0 , X + sin ( β 1 , X ) + ε , where β 0 and β 1 are orthogonal functions. When E [ sin ( β 1 , X ) ] = 0 , then E [ Y g ( β 0 , X ) | β 0 , X ] = 0 providing that β 0 , X is independent of β 1 , X and X is a Gaussian process with E ( X ) = 0 ; similar examples can be found in [14]. Thus, the test based on E [ Y g ( β 0 , X ) | β 0 , X ] cannot detect the above alternative. Therefore, we can construct a test statistic based on E [ Y g ( β 0 , X ) | B , X ] .
Define ε = Y g ( X , β ) , E ( ε | B , X ) = 0 almost surely under H 0 and P ( E ( ε | B , X ) = G ( B , X ) g ( β 0 , X ) 0 ) > 0 under H 1 . Thus, it is natural to construct a test statistic based on the orthogonality between model’s residuals ε and a non-parametric estimation of E ( ε | B , X ) . Thus, we construct a test statistic based on E { ε E ( ε | B , X ) } = E { E 2 ( ε | B , X ) } .
Given the sample { X i , Y i } i = 1 n and B ^ with its structural dimension K ^ (the estimation procedure for B ^ and K ^ will be introduced in Section 3), the kernel regression estimator of E ( ε i | B , X i ) is
E ^ ( ε i | B ^ , X i ) = j i n ε j K h ( B ^ , X i X j ) j i n K h ( B ^ , X i X j ) ,
where K h ( · ) = K ( · / h ) / h K ^ is a K ^ -dimensional kernel function, and h is a bandwidth.
When g ( · ) is known, similar to the estimation methods in [25] (Refs. [26,27,28]’s methods are also available), it is easy to obtain a consistent estimator of β , which is denoted as β ^ . Define ε ^ i Y i g ( β ^ , X i ) . Next, to simplify the formula, let ε ^ = ( ε ^ 1 , , ε ^ n ) , U ^ be the n × n diagonal matrix with ( i i ) th element ε ^ i , w ^ i j K h ( B ^ , X i X j ) / l i n K h ( B ^ , X i X l ) and w ^ i i = 0 , W ^ n = { w ^ i j } be the n × n matrix with ( i , j ) th element w ^ i j , where B ^ is an FSDR estimate by [18] in following Section 3. Then, by Proposition 1 and Corollary 1 in [23], the finite sample corrected and standardized quadratic statistic is
T n F S I R = ε ^ W ^ n s ε ^ 2 s ( U ^ W ^ n s U ^ ) + K ^ 2 s ( W ^ n s ) ,
where W ^ n s = ( W ^ n + W ^ n ) / 2 which satisfies the symmetry condition required in Corollary 1 of [23], s ( A ) = ( i , j a i j 2 ) 1 / 2 is the Frobenius norm of matrix A = ( a i j ) . We call K ^ 2 s ( W ^ n s ) the finite sample correction term (FSC). When T n F S I R is large enough, the null hypothesis is rejected.
Remark 1. 
This paper mainly extends the problem considered in [14] to the functional data model. In addition to this, there is a clear difference in the construction of the test statistic. For example, [14] constructed the test statistic based on E { ε E ( ε | B X ) ω ( B X ) } for a positive weight ω ( B X ) . When ω ( B X ) was the denominator of the non-parametric estimate of E ( ε | B X ) , they obtained a simple statistic. Thus, T n F S I R is different from the test of [14] even without the standardization and the bias correction.

3. Identification and Estimation of B and K

Functional sufficient dimension reduction (FSDR) is concerned with the situations where the distribution of Y given X depends on X through a set of linear combinations of X; that is, there is a minimal set of members B = ( β 1 , , β K ) in H , such that Y X | B , X . In the applications where the conditional mean is of main interest, we can formulate FSDR through conditional mean as follows
Y E ( Y | X ) | B , X ,
where ⊥ stands for the conditional independence of Y and E ( Y | X ) given B , X . This is a fundamental assumption for functional sufficient dimension reduction [19,22,29]. It is worth noting that Y E ( Y | X ) | B , X is equivalent to E ( Y | X ) = E ( Y | B , X ) by Proposition 8.1 in [29]. For any K × K orthogonal matrix D , the above relation is unchanged if we replace B by D B , since E ( Y | B , X ) = E ( Y | D B , X ) . In fact, the identifiable parameter is the space spanned by the columns of B .
We denote the column space of B by s p a n ( B ) . S E [ Y | X ] is defined as the central mean subspace and is the intersection of all subspace s p a n ( B ) such that Y E ( Y | X ) | B , X . S Y | X is defined as the central subspace and is the intersection of all subspace s p a n ( B ) of minimal dimension such that Y X | B , X . Most of the FSDR methods are looking for central subspace S Y | X , such as the FSIR method in [19], FSAVE method in [22], and FSDR via distance covariance in [30]. As proposed in [29] Theorem 8.1, S E ( Y | X ) is a subspace of S Y | X . To ensure that S Y | X = S E ( Y | X ) for ε in model (2), we assume that ε = G 1 ( B , X ) ε ˜ , where G 1 ( · ) is an unknown smooth function and ε ˜ X . ε X is a special case. Next, we uniformly use S Y | X to represent the FSDR space spanned by B . Here, we only introduce a popular FSDR method, functional sliced inverse regression (FSIR) which refers to [18], as an example. Of course, other FSDR methods are also applicable.

3.1. A Brief Review on FSIR

Most of the FSDR methods are based on the functional multiple index model (2). As introduced in following Proposition 1, with the FSDR technique, we can achieve with probability tending to 1 K ^ = 1 under H 0 ; K ^ 1 under H 1 . In the following, we will introduce the specific dimensionality reduction technique of FSIR.
Given ( X , Y ) , by Karhunen–Loéve expansion of X, we have
X = j = 1 ξ j ϕ j ,
where E ξ j 2 = γ j , γ i and ϕ j are the eigenvalues and eigenfunctions, respectively, of the covariance operator of X, that is Γ = V a r ( X ) = E ( X X ) , where for x , y L 2 ( [ 0 , 1 ] ) , x y : L 2 ( [ 0 , 1 ] ) L 2 ( [ 0 , 1 ] ) is defined by ( x y ) ( z ) = x , z y . Define Γ S Y | X as the space spanned by Γ β 1 , , Γ β K . In addition, we write Γ V a r ( X | Y ) Γ S Y | X when ( Γ V a r ( X | Y ) ) β Γ S Y | X for all β L 2 ( [ 0 , 1 ] ) . To ensure that the estimation of FSIR belongs to the space of Γ S Y | X , we need the following Lemma (refer to Theorem 1 in [18]).
Lemma 1. 
(a) Suppose for all b L 2 ( [ 0 , 1 ] ) , the conditional expectation is a linear condition E { b , X | B , X } = c 0 + i = 1 K c i β i , X for constants c 0 , , c K , then E ( X | Y ) Γ S Y | X ; (b) If in addition V a r ( X | B , X ) is nonrandom, Γ V a r ( X | Y ) Γ S Y | X .
The conditions in (a) and (b) hold when X is a Gaussian process, although Gaussianity is not necessary. Furthermore, there is no requirement for specific distribution forms in the following theoretical study. In the simulation of Section 4, we will use Gaussian and other continuous distributions.
Since for functional data, the operator Γ 1 is not bounded, we replace Γ 1 by ( Γ + ρ I ) 1 , where ρ 0 as n . Then, the FSDR space, which is obtained by FSIR, is the space spanned by the top eigenvectors of
( Γ + ρ I ) 1 V a r ( E [ X | Y ] ) .
Given independent and identically distributed samples { ( X i , Y i ) } i = 1 n , referring to [29]’s steps in Section 3 on vector data, we summarize the estimation procedure of functional data as the following algorithm (Algorithm 1).
Algorithm 1: Estimation of B
step 1. Compute the sample mean and sample variance functions X ¯ ( t ) = 1 n i = 1 n X i ( t ) and Γ ^ n = 1 n i = 1 n [ X i ( t ) X ¯ ( t ) ] [ X i ( t ) X ¯ ( t ) ] , where t , t [ 0 , 1 ] , then, define Z i ( t ) = X i ( t ) X ¯ ( t ) , i = 1 , , n .
step 2. Discretize Y: let ( I s ˙ ) s ˙ = 1 S ˙ be slice intervals in [ min ( Y ) , max ( Y ) ] , define p s ˙ = P ( Y I s ˙ ) and e s ˙ = E ( Z | Y I s ˙ ) , v s ˙ = V a r ( Z | Y I s ˙ ) , then, p ^ s ˙ = E n [ I ( Y I s ˙ ) ] , e ^ s ˙ = i = 1 n Z i ( t ) I ( Y i I s ˙ ) / n s ˙ , where E n ( · ) is the empirical distribution function, n s is the number of elements that satisfy Y i I s ˙ .
step 3. Approximate V a r [ E ( Z | Y ) ] by Γ ^ S ˙ F S I R = s ˙ = 1 S ˙ p ^ s ˙ e ^ s ˙ e ^ s ˙ / S ˙ .
step 4. Let β ^ 1 , , β ^ r and λ ^ 1 , , λ ^ r be the top r eigenvectors and eigenvalues of ( Γ ^ n + ρ I ) 1 Γ ^ S ˙ F S I R , where r K and λ ^ 1 λ ^ 2 λ ^ r 0 . The sufficient predictors are
β ^ 1 , Z , , β ^ r , Z .
Remark 2. 
The way of discretization in step 2 may affect the outcome of Algorithm 1. A lot of discussions on the discretization were carried out, and some methods were recommended; see [18,31,32], among others.

3.2. Estimating the Structural Dimension K ^

There are several methods to estimate the structural dimension K of S Y | X for vector data. Ref. [33] introduced a Bayesian information criterion (BIC), Refs. [34,35] showed a modified BIC method, and Refs. [36,37] used a bootstrap method. Through simulation, we find that [34]’s method verifies the conclusion of Proposition 1 better. That is,
D ( k ) = n 2 l = 1 k { log ( λ ^ l + 1 ) λ ^ l } l = 1 r { log ( λ ^ l + 1 ) λ ^ l } 2 C n k ( k + 1 ) 2 r .
where λ ^ 1 λ ^ 2 λ ^ r 0 are eigenvalues of ( Γ ^ n + ρ I ) 1 Γ ^ S ˙ F S I R and r K .
We maximize D ( k ) over k = 1 , , r , and denote the maximizer as K ^ . Then, if  C n / n 0 and C n , the estimate K ^ by (6) satisfies: P ( K ^ = 1 ) 1 under H 0 ; P ( K ^ = K 1 ) 1 under H 1 ; refer to Theorem 4 in [34]. The choice of C n = n can ensure the consistency of the estimators under hypotheses H 0 and H 1 .

4. Asymptotic Properties

Define semi-metric d β ( x 1 , x 2 ) = | β , x 1 x 2 | . The following assumptions are used to prove the asymptotic properties of the test statistic.
Assumption 1. 
E X 4 < , Γ has eigenvalues λ 1 > λ 2 > > 0 , and  E ξ j 4 C λ j 2 .
Assumption 2. 
For some constants a , b that satisfy a > 1 and a 1 / 2 < b , j a C λ j , | ξ k j * | C j b , for some C > 0 , where ξ k j * = β k , ϕ j , j = 1 , 2 , are the generalized Fourier coefficients of β k , 1 k K .
Assumption 3. 
The regularization parameter is chosen as ρ n a / ( a + 2 b ) .
Assumption 4. 
h = h n is a sequence of positive numbers satisfying C 2 n τ 2 < h < C 1 n τ 1 with 0 < τ 1 < τ 2 < 1 .
Assumption 5. 
K is an asymmetrical strictly dereasing kernel in [ 0 , 1 ] such that c 1 , c 2 > 0 , c 1 I [ 0 , 1 ] ( t ) < K ( t ) < c 2 I [ 0 , 1 ] ( t ) for t R .
Assumption 6. 
(1) C is a compact subset of H such that X C , a.s.; (2) β Θ n such that c a r d Θ n = n α with α > 0 , we have P ( d β ( X 1 , X 2 ) < h | X 1 ) = C X 1 , β h + o ( h ) , a.s., where X 1 and X 2 are independent copies of X; (3) 0 < inf β Θ n C X 1 , β sup β Θ n C X 1 , β < , a.s.
Assumption 7. 
Let g β ( · ) = E [ Y | β , · ] . g β satisfies, for  β , the Holder-type condition: C > 0 , b 1 > 0 , ( x 1 , x 2 ) C 2 , | g β ( x 1 ) g β ( x 2 ) | C d β ( x 1 , x 2 ) b 1 .
Assumption 8. 
For any positive integer c 3 , E ( | Y | c 3 | X ) C c 3 , X < , a.s. and ϱ > 0 , E ( Y 2 | X = x ) = σ ( x ) ϱ , with  σ ( · ) continuous on C .
Assumption 9. 
The density f ( B , X ) of B , X on its support C 1 exists and has two bounded derivatives and satisfies 0 < inf B , X C 1 f ( B , X ) sup B , X C 1 f ( B , X ) < .
Remark 3. 
Assumptions 1–3 are used to obtain β ^ j β j 2 = O p ( n ( 2 b 1 ) / ( a + 2 b ) ) O p ( r n 2 ) , j = 1 , , K and r n = n 2 b 1 2 ( a + 2 b ) , which was shown in [18]. By choosing the regularization parameter as ρ n a / ( a + 2 b ) in Assumption 3, we aim to balance the bias and variance of the estimator and avoid overfitting. The choice of ρ in this form is motivated by the theoretical analysis of the convergence rates of the FSIR methods. Assumptions 4–8 mainly come from [1], which are necessary assumptions to obtain the following Theorem 1. Assumption 9 is needed for the asymptotic normality of our statistic.
For a given g ( · ) , ε ^ i = Y i g ( β ^ , X i ) . Similar to the proof of [25], under Assumptions 1–3 and H 0 , we have β ^ β 2 = O p ( r n 2 ) for least squares estimate (or penalized likelihood estimation in [38,39], but other assumptions may be needed). Therefore, the limiting behavior of the test statistic can be obtained.
Theorem 1. 
Given Assumptions 1–9, if  n h 1 / 2 r n 2 0 and r n / h 0 as n , then, under  H 0 , we have
T n F S I R d N ( 0 , 1 )
Proof. 
See Appendix A (Part I).    □
Remark 4. 
It should be noted that although theoretically, we find that F S C p 0 as n , through simulation, we found that FSC can improve the power to some degree for finite samples.
We now investigate the power performance of the test statistic T n under the fixed alternatives model (2). Suppose that under H 1 , there exists a unique function θ L 2 ( [ 0 , 1 ] ) satisfying θ = a r g m i n β E ( Y g ( β , X ) ) 2 ; then, ε ^ i = Y g ( θ , X i ) { g ( β ^ , X i ) g ( θ , X i ) } .
Theorem 2. 
Under the alternative model (2) and Assumptions 1–9, it holds that β ^ θ = O p ( r n ) and
T n F S I R .
Proof. 
See Appendix A (Part II).    □
Next, we consider the following sequence of the local alternative model as
H 1 n : Y = g ( β 0 , X ) + C 1 n G ( B , X ) + ε ,
where C 1 n goes to zero, G ( · ) satisfies E [ G 2 ( B , X ) ] < and inf β L 2 ( [ 0 , 1 ] ) E [ G ( B , X ) g ( β , X ) ] 2 > 0 .
In order to obtain the main results about the power behaviors under H 1 n of (7), we need to show the asymptotic property of K ^ when C 1 n 0 .
Proposition 1. 
Under the local alternative hypothesis H 1 n in (7) and Assumptions 1–3, when C n h 1 / 2 ( C n was introduced in (6)), if  C 1 n = O ( n 1 / 2 h 1 / 4 ) , we have P ( K ^ = 1 ) 1 .
Proof. 
See Appendix A (Part III).    □
We now present the power performance of our tests under the local alternative hypothesis. Before this, the following condition needs to be met.
Similar to the parameter estimation method under H 0 , we can obtain the parameter estimate β ^ with observations { X i , Y i } i = 1 n under H 1 n . Then, we have ε ^ i = ε i + C 1 n G ( B , X ) { g ( β ^ , X i ) g ( β 0 , X i ) } .
Theorem 3. 
If Assumptions 1–9 hold, n h 1 / 2 r n 2 0 and r n / h 0 as n . Then, under  H 1 n , if  C 1 n = n 1 / 2 h 1 / 4 and β ^ β 0 = O p ( r n ) , we have
T n F S I R d N ( ω , 1 ) ,
where E [ G 2 ( B , X ) ] 2 σ 2 E [ f 1 ( B , X ) ] , σ 2 = V a r ( ε ) .
Proof. 
See Appendix A (Part III).    □

5. Simulations Results and Real Data Analysis

We now perform simulations to check the performance of the proposed test. In the non-parametric estimation, we use the Quartic kernel function K ( u ) = 15 / 16 ( 1 u 2 ) 2 I ( u [ 0 , 1 ] ) , and the bandwidth h is selected by generalized cross-validation (GCV), which is similar to the method in [40]. Similar to the kernel estimators in Euclid spaces, the choice of the kernel affects less, and other kernel functions are also applicable, such as Gaussian kernel function and Triangular kernel function. The reason for choosing the Quartic kernel function is usually because it has a better fitting effect than the Triangular kernel, and compared to Gaussian kernel function, the Quartic kernel function is relatively more flexible in bandwidth selection. In the functional regression analysis, the optimal bandwidth h can be chosen by the R function f r e g r e . n p in the f d a . u s c package. Since the conclusion in [18] shows that the results are less sensitive to the slice number S ˙ , we choose S ˙ = 3 in the simulations if not specifically stated. The regularization parameter is chosen as ρ = 0.02 , which is suggested by the results in [18]. All studies are based on 1000 repetitions, and the significance level is always taken as α = 0.05 . In the following subsections, we check the adequacy of the parametric functional single-index model by considering that g ( · ) is both linear and nonlinear.

5.1. Study 1: Linear Link Function

When g ( · ) is a linear function, we need to check the adequacy of the functional linear models (FLM). Among the existing model testing methods for FLM, [13] constructed two test statistics CvM and KS based on a random projection empirical process. We also compare the proposed test with the weighted kernel smoothing test suggested in [11], and we denote the test by S n when its critical value is directly given by the asymptotic normal distribution and  S n B when its critical value is obtained by bootstrap. In addition, we also extend the method of [11] to sufficient dimension reduction with the kernel function K ( B , X i X j / h ) and denote the test as S n FSIR (see the following Algorithm 2). Similar to S n and S n B , we can drop out the FSDR and just use a univariate kernel K ( X i X j / h ) to substitute the multivariate kernel in T n FSIR ; then, we obtain the corresponding statistics as T n and T n B , respectively. In the simulation, the KS and CvM approach can be directly implemented by the r p . f l m . t e s t function in R package f d a . u s c when g ( · ) is a linear link function, and the number of projections is taken as 3. The bootstrap resampling number is 500 in each case. The specific steps to calculate the test statistic T n FSIR and S n FSIR are as follows:   
Algorithm 2: Calculate T n FSIR and S n FSIR
step 1. To obtain ε ^ in T n FSIR and S n FSIR , we use the p f r function in the r e f u n d package of R for the FLM parameter estimation.
step 2. B ^ in the kernel function K h ( B ^ , X i X j ) of the test statistic T n FSIR is estimated by Algorithm 1.
step 3. K ^ is determined by (6). When K ^ = 1 , K h ( · ) is a one-dimensional kernel function; when K ^ > 1 , K h ( · ) is the product of K ^ one-dimensional kernel functions.
step 4. Then, we compute S n F S I R = 1 n ( n 1 ) j i n 1 h K h ( B ^ , X i X j ) ε ^ i ε ^ j and T n F S I R in (5).
We consider the following functional linear link model
Y = X , β 1 + c · m X , B + ε ,
where X is a standard Brownian motion on [ 0 , 1 ] and is observed on a grid of equally spaced 101 points, β j = sin ( j π t ) is the Fourier orthogonal basis with j = 1 , 2 , , ε N ( 0 , 0 . 1 2 ) , c is a constant, and the null hypothesis H 0 holds when c = 0 . We carry out simulations from different perspectives with two examples.
Example 1. 
In this example, our aim is to examine the performance of different test statistics under following alternative models:
  • H 11 : Y = X , β 1 + 0.2 c · cos ( 2 π X , β 1 ) + ε ;
  • H 12 : Y = X , β 1 + 0.2 c · exp ( X , β 1 ) + ε ;
  • H 13 : Y = X , β 1 + 0.2 c · X , β 1 X , β 2 + ε ;
  • H 14 : Y = X , β 1 + c exp ( X , β 2 ) + cos ( X , β 3 ) + ε ;
  • H 15 : Y = X , β 1 + c exp ( X , β 2 ) + cos ( X , β 3 ) + X , β 4 2 + ε ;
  • H 16 : Y = X , β 1 + c exp ( X , β 2 ) + cos ( X , β 3 ) + X , β 4 2 + X , β 5 3 + ε ,
where H 11 and H 12 are used to examine whether and how much power our test loses by sufficient dimension reduction. H 13 is the case of the multiplicative effect, H 14 H 16 are used to analyze the effect of dimensionality on the different test statistics.
Table A1 shows the empirical sizes and powers of different test statistics with the nominal size α = 0.05 . Through the comparative analysis in different scenarios, we conclude that: (1) when the critical value is directly given by the asymptotic normal distribution, the empirical sizes of T n FSIR are relatively closer to the pre-specified significance level α = 0.05 , while S n and T n are more conservative. (2) T n FSIR has the highest power, except for small and moderate cs in the cases of H 12 and H 13 . (3) S n , CvM and KS test statistics are mainly applicable to the cases of H 11 and H 12 , see [11,13]. In H 11 and H 12 , B = β 1 under both the null and alternative hypothesis. Therefore, FSDR is not needed in these two alternative models, but we can examine whether and how much the power would be lost if FSDR were used. Comparing the results of using FSDR and not using FSDR in Table A1 for H 11 and H 12 , it is evident that the use of FSDR would result in some loss of power in H 12 , while T n FSIR still has the highest power under H 1 n . (4) Under the multiplicative effects model H 13 , the considered test statistics remain effective, and the test statistic using FSIR is slightly more powerful for large c and n. (5) From H 14 to H 16 , the power of the considered test statistic increases to varying degrees with the increasing dimensionality. (6) From H 14 to H 16 , there is a significant improvement in the powers of almost all test statistics. This indicates that the test statistics are more sensitive to the quadratic model than the high-frequency oscillation model. The power is slightly improved after adding a cubic model from H 15 to H 16 . (7) The powers of all statistics increase significantly as the sample size increases, and the power of the FSIR method increases more rapidly with increasing c, especially for T n FSIR .
Example 2. 
To analyze the effect of different data from what we considered in Example 1, we take model Y = X , β 1 + c exp ( X , β 2 ) + ε as an example to analyze the performance of different statistics under the following data generation mechanism (DGM):
  • D G M 1 : X O r n s t e i n U h l e n b e c k p r o c e s s , β 1 = sin ( π t ) , β 2 = sin ( 2 π t ) ;
  • D G M 2 : X s t a n d a r d B r o w n i a n m o t i o n , β 1 = sin ( π t ) , β 2 = sin ( 2 π t ) ;
  • D G M 3 : X s t a n d a r d B r o w n i a n m o t i o n , β 1 = 2 ( 2 t 1 ) 2 1 , β 2 = 4 ( 2 t 1 ) 3 3 ( 2 t 1 ) ;
  • D G M 4 : X s t a n d a r d B r o w n i a n m o t i o n , β 1 = 2 ( 2 t 1 ) 2 1 , β 2 = sin ( 2 π t ) ,
where ε N ( 0 , 0 . 1 2 ) , β 1 and β 2 in DGM 1 DGM 2 and DGM 3 are the Fourier orthogonal basis and Chebyshev orthogonal basis, respectively, and DGM 4 has a non-orthogonal basis.
We present the test results for Example 2 in Table A2. In DGM 1 and DGM 2 , although X is different, the test results are similar. Comparing the results of DGM 2 DGM 4 , which has different β s, we find that the powers in the Chebyshev basis are lower than those in the Fourier basis, while the powers in the case of the non-orthogonal basis are somewhere in between.

5.2. Study 2: Non-Linear Link Function

When g ( · ) is a given nonlinear function, we only give results of T n , T n B and T n FSIR , since S n , S n B , T n , T n B , CvM and KS are studied for FLM in [11,13]. Except for step 1 in Algorithm 2, the other steps are similar. Step 1 can be implemented in a similar way to functional principal component analysis in [25] when g ( · ) is known, and ε ^ can be obtained through nonlinear least squares estimation.
We consider the nonlinear link function model:
Y = g ( X , β 1 ) + c · m X , B + ε ,
where g ( · ) is a smooth function. Similar to Study 1, we mainly analyze the performance of the test under different alternative models. The data generation mechanisms are the same as those in DGM 2 . In order to make our test more widely applicable, we also analyze the case of monotonic and non-monotonic link functions, as demonstrated by the following two examples, respectively.
Example 3. 
In this example, we consider the case where g ( · ) is monotonic. The following models with the exponential link function are implemented.
  • H 31 : Y = exp ( X , β 1 ) + c · cos ( X , β 2 ) + ε ;
  • H 32 : Y = exp ( X , β 1 ) + c · cos ( X , β 2 ) + | X , β 3 | + ε ;
  • H 33 : Y = exp ( X , β 1 ) + c · cos ( X , β 2 ) + | X , β 3 | + X , β 4 2 + ε .
Example 4. 
To examine the power of the test when g ( · ) is non-monotonic, we analyze the following models with the sin link function,
  • H 41 : Y = sin ( X , β 1 ) + c · exp ( X , β 2 ) + ε ;
  • H 42 : Y = sin ( X , β 1 ) + c · exp ( X , β 2 ) + | X , β 3 | + ε ;
  • H 43 : Y = sin ( X , β 1 ) + c · exp ( X , β 2 ) + | X , β 3 | + X , β 4 2 + ε .
We present the results for Examples 3 and 4 in Figure A1. We make three observations. First, for the testing of nonlinear models, the statistic constructed based on the FSIR method always has higher power than others. Second, the difference between monotonic and non-monotonic power is not significant, and the effect of dimensionality on power is not very large, but the differences between T n FSIR and T n B , T n in the monotonic cases are more obvious than the non-monotonic cases. For example, the power increment of T n FSIR relative to T n in H 33 is slightly larger than that in H 43 when n = 200 . Third, there is a significant increase in power with increasing sample size.
To summarize, through the above examples, including different data generation mechanisms, multiplicative effect models, and monotone/non-monotone cases under functional data, we find that the use of the FSDR method in functional data not only increases the power with the increase of sample size n and deviation parameter c but also improves the power of the test to some extent compared with the method without using FSDR [11,13] even under the multiplicative effect model. The above conclusion is largely consistent with the findings in vector data analysis [14,41,42]. In functional data analysis, similar research has been conducted, such as [43], but Ref. [43]’s method is only applicable to linear models. Our study is not only effective in linear situations but also performs well in nonlinear cases.

5.3. Analysis of the COVID-19 Data

As one of the most important aggregate indicators reflecting the overall macroeconomic operation of each country, Gross Domestic Product (GDP) is an important reference value for decision makers to make rational decisions based on the economic situation. The large shock to the world’s macroeconomy caused by the COVID-19 crisis has imposed tremendous challenges on now accurately forecasting a country’s GDP. Performing relevant data modeling is an important way to make predictions. Goodness-of-fit tests are a prerequisite for accurate modeling. Therefore, when the relationship between the number of new COVID-19 cases per day and the GDP of a country is known, the model can be tested first using the method in this paper, and then the related statistical analysis can be performed. Below, we describe the specific operations.
If X and Y have the following relationship model
E [ Y i G D P | X i ] = g X i , β ,
where Y i is the DP growth (annual %) (https://advisor.visualcapitalist.com/gdp-growth-by-country-in-2021/, accessed on 18 January 2023) of country i in 2021, X i ( t ) is the rolling mean of newly confirmed COVID-19 cases (in R package t i d y c o v i d 19 ) per day as reported by JHU CSSE (https://www.worldbank.org/en/home, accessed on 18 January 2023) in 2021, i = 1 , 2 , , 183 . Rolling mean can be implemented by the r o l l m e a n function in the R language and the width of the rolling window is 7; missing values are replaced by the nearby rolling mean. The specific trends of X and Y and the selected countries are shown in Figure A2, where the legend shows the ISO numbers for each country. In the legend of X, we have only indicated the three countries with a high number of confirmed cases. It is worth mentioning that there are more than 183 countries in the world, and we selected 183 countries for our analysis mainly because we excluded those countries that did not publish their GDP in 2021 and had zero cumulative confirmed cases in 2021.
Using the test statistic and parameter values, such as bandwidth and slice number, as mentioned in the simulation, if g is a linear function, the test results are displayed in Table A3. Table A3 shows that except for CvM and K S , all other statistics reject functional linear models, where K ^ = 3 with FSIR . These conclusions suggests that a nonlinear model should be used: for example, a functional multi-index model may be appropriate. However, CvM and KS fail to detect this nonlinear model structure.

6. Conclusions and Discussion

In this paper, we extend [14]’s method to functional data and construct a different test statistic for checking a single-index model. We theoretically prove the large-sample properties of the constructed test statistic under the null hypothesis, alternative hypothesis, and local alternative hypothesis. In addition, the simulations of Section 4 also show that the FSIR-based statistic has higher power compared to other statistic in most cases and can be used to check the adequacy of functional linear or nonlinear models.
We can extend our methodology to missing, censored or dependent data. For example, the case of missing data can be referred to the method of [41]. Of course, in addition to mean regression, we can also consider the quantile regression or composite quantile regression.

Author Contributions

L.X.: Conceptualization, Formal analysis, Methodology, Software, Writing—original draft, Data curation; T.L.: Formal analysis, Methodology, Software, Writing—review and editing; Z.Z.: Conceptualization, Formal analysis, Methodology, Investigation, Supervision, Funding acquisition, Writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research is partly supported by the National Natural Science Foundation of China (Grant No. 12271014 and No. 11971045).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

“GDP growth (annual %) Data” at https://advisor.visualcapitalist.com/gdp-growth-by-country-in-2021/ (accessed on 18 January 2023), “COVID-19 Data (in R package t i d y c o v i d 19 )” at https://www.worldbank.org/en/home (accessed on 18 January 2023).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1. Part I

Proof of Theorem 1. 
We first give several notations used in the proof. Specifically, Σ n is an n × n diagonal matrix with σ i i = V a r 1 / 2 ( ε i ) , and for any n × n matrix A n , denote the spectral radius of A n by
ρ ( A n ) sup v R n , v 0 A n v v .
When A n is symmetric and the eigenvalue satisfies | k 1 | | k 2 | | k n | , ρ ( A n ) = | k 1 | .
Next, we prove Theorem 1 in five steps.
Step 1: ε W n s ε 2 s ( Σ n W n s Σ n ) d N ( 0 , 1 ) .
Both W n and W n s are n × n matrices, with their ( i j ) th elements being w i j = K h ( B , X i X j ) / l i n K h ( B , X i X l ) and ( w i j + w j i ) / 2 , respectively.
To present the proof of Step 1, we first introduce two pivotal lemmas. The first lemma is a direct conclusion of Theorem 1.1 in [44].
Lemma A1. 
Assume that { ε } is a sequence of independent uniformly square integrable random variables with mean zero, and ε = ( ε 1 , , ε n ) . Define Σ n as an n × n diagonal matrix with ii-th element σ i i . For a real symmetric matrix A n = ( a i j ) n × n with zero diagonal, if ρ ( A n ) / s ( A n ) 0 as n , then
T = ε A n ε 2 s ( Σ n A n Σ n ) d N ( 0 , 1 ) .
The second lemma comes from [23].
Lemma A2. 
Let X n = ( X 1 , X 2 , , X n ) be an i.i.d. sample from X, and x n for a realization of X n . ε is the same as in Lemma 1. P n ( X n ) and T ( X n , ε ) are measurable functions. If
P n ( x n ) 0 T ( X n , ε ) d N ( 0 , 1 ) ,
and
P n ( X n ) p 0 ,
then T ( X n , ε ) d N ( 0 , 1 ) .
By Lemmas A1 and A2, to prove the conclusion in Step 1, we need to prove ρ ( W n s ) / s ( W n s ) p 0 as n . Firstly, since W n and W n have the same eigenvalues, it is easy to show
ρ ( W n s ) = ρ ( W n ) .
Secondly, noting that w i j 0 , it holds
s 2 ( W n s ) 1 2 s 2 ( W n ) .
Thus,
ρ ( W n s ) s ( W n s ) 2 ρ ( W n ) s ( W n ) ,
and it suffices to show that ρ ( W n ) / s ( W n ) p 0 .
Since the row sum of W n is 1, ρ ( W n ) 1 . Thus, we need only show s ( W n ) p . Define Z = B , X j as a random variable and the density function exists (Assumption 9). Since
E j = 1 n 1 K h 2 ( B , X j x ) = ( n 1 ) E K h 2 ( Z B , x ) = ( n 1 ) K 2 Z B , x h f ( Z ) d Z = ( n 1 ) K 2 ( u ) f u h + B , x h K d u = ( n 1 ) h K K 2 ( u ) f B , x + h u f ( B , x ) + o ( h ) d u = ( n 1 ) h K f ( B , x ) K 2 ( u ) d u + ( n 1 ) h K + 1 u f ( B , x ) K 2 ( u ) d u + o ( n h K + 1 ) .
Similarly, we have
E j = 1 n 1 K h ( B , X j x ) = ( n 1 ) E K h ( B , X x ) | x = ( n 1 ) h K f ( B , x ) + ( n 1 ) h K + 1 f ( B , x ) u K ( u ) d u + o ( n h K + 1 ) .
These two equations lead to
s 2 ( W n ) = i j w i j 2 = i = 1 n 1 n h K f ( B , X i ) K 2 ( u ) d u + o p ( 1 ) = 1 h K E [ f 1 ( B , X ) ] K 2 ( u ) d u + o p ( 1 ) = O p ( h K ) .
It is easy to derive s ( W n ) = O p ( h K / 2 ) . This result will be used in the following steps.
By Assumption 5: c 1 , c 2 > 0 , c 1 I [ 0 , 1 ] ( t ) < K ( t ) < c 2 I [ 0 , 1 ] ( t ) for t R , we have
E j = 1 n 1 K h 2 ( B , X j x ) ( n 1 ) h K f ( B , x ) c 1 2 + o p ( 1 ) .
Thus,
s 2 ( W n ) 1 h K E [ f 1 ( B , X ) ] c 1 2 + o p ( 1 ) .
Under Assumption 4, when K = 1 under H 0 , we have s ( W n ) .
Step 2: ε W ^ n s ε ε W n s ε 2 s ( Σ n W n s Σ n ) p 0 .
To begin, note that
ε W ^ n s ε ε W n s ε = i = 1 n j i n ( w ^ i j s w i j s ) ε i ε j = 1 2 i = 1 n j i n K h ( B ^ , X i X j ) l i n K h ( B ^ , X i X l ) K h ( B , X i X j ) l i n K h ( B , X i X l ) + K h ( B ^ , X i X j ) l j n K h ( B ^ , X j X l ) K h ( B , X i X j ) l j n K h ( B , X j X l ) ε i ε j = 1 2 i = 1 n j i n K h ( B ^ , X i X j ) l i n K h ( B ^ , X i X l ) K h ( B , X i X j ) l i n K h ( B , X i X l ) ε i ε j + i = 1 n j i n K h ( B ^ , X i X j ) l j n K h ( B ^ , X j X l ) K h ( B , X i X j ) l j n K h ( B , X j X l ) ε i ε j 1 2 ( D 1 n + D 2 n )
Then, an application of Taylor expansion of D 1 n yields
D 1 n = i = 1 n j i n K h ( B ^ , X i X j ) l i n K h ( B ^ , X i X l ) K h ( B , X i X j ) l i n K h ( B , X i X l ) ε i ε j = D 1 n + o p ( D 1 n ) ,
where
D 1 n = i = 1 n j i n K h ( B , X i X j ) l i n K h ( B , X i X l ) K h ( B , X i X j ) l i n K h ( B , X i X l ) l i n K h ( B , X i , X l ) 2 · B ^ B , X i X j h K ^ ε i ε j = i = 1 n j i n 1 h K ^ s i j B ^ B , X i X j ε i ε j ,
and superscript “′” represents the gradient vector, B ^ B , X i X j = β ^ 1 β 1 , X i X j , β ^ 2 β 2 , X i X j , , β ^ K β K , X i X j and we define s i j as the K-dimensional column vectors.
Similar to D 1 n , D 2 n = D 2 n + o p ( D 2 n ) , where
D 2 n = i = 1 n j i n K h ( B , X i X j ) l j n K h ( B , X j X l ) K h ( B , X i X j ) l j n K h ( B , X j X l ) l j n K h ( B , X j X l ) 2 · B ^ B , X i X j h K ^ ε i ε j i = 1 n j i n 1 h K ^ ( s i j ) B ^ B , X i X j ε i ε j ,
and s i j is the K-dimensional column vector corresponding to the coefficients of B ^ B , X i X j h K ^ ε i ε j .
Let the elements in the K-dimensional vector 1 / 2 ( s i j + s i j ) X i X j be a i j k , A n k is a n × n matrix with the element a i j k . Then,
ε W ^ n s ε ε W n s ε 2 s ( Σ n W n s Σ n ) = 1 2 ( D 1 n + D 2 n ) + 1 2 o p ( D 1 n ) + o p ( D 2 n ) ,
where
D 1 n + D 2 n = i = 1 n j i n 1 2 h K ^ ( s i j + s i j ) B ^ B , X i X j ε i ε j 2 s ( Σ n W n s Σ n ) i = 1 n j i n 1 2 ( s i j + s i j ) ε i ε j X i X j 2 s ( Σ n W n s Σ n ) · B ^ B h K ^ = i = 1 n j i n 1 2 ( s i j + s i j ) X i X j ε i ε j d i a g ( 2 s ( Σ n A n k Σ n ) 1 ) · d i a g ( 2 s ( Σ n A n k Σ n ) ) 2 s ( Σ n W n s Σ n ) · B ^ B h K ^ ,
where B ^ B = ( β ^ 1 β 1 , β ^ 2 β 2 , , β ^ K β K ) .
Similar to the proof in Step 1, the elements of the K-dimensional vector i = 1 n j i n 1 2 ( s i j + s i j ) X i X j ε i ε j · d i a g ( ( 2 s ( Σ n A n k Σ n ) ) 1 ) tend to a multivariate standard normal distribution. Thus, the first term of (A2) is O p ( 1 ) .
We define the k-th component of the K-dimensional vector s i j as s i j k , k = 1 , , K , where s i j comes from D 1 n . Letting the K-dimension vector R = i = 1 n j = 1 n ( s i j k ) 2 X i X j 2 , it is easy to see that 2 s ( Σ n A n k Σ n ) and R 1 / 2 have the same convergence rate. By Assumption 5, we have
R = i = 1 n j i n K h ( B , X i X j ) l i n K h ( B , X i X l ) K h ( B , X i X j ) l i n K h ( B , X i X l ) 2 X i X j 2 l i n K h ( B , X i , X l ) 4 = i = 1 n j = 1 , j i n [ K h ( B , X i X j ) ] 2 X i X j 2 l i n K h ( B , X i , X l ) 2 + i = 1 n j = 1 , j i n K h 2 ( B , X i X j ) l i n K h ( B , X i X l ) 2 X i X j 2 l i n K h ( B , X i , X l ) 4 2 i = 1 n j = 1 , j i n K h ( B , X i X j ) K h ( B , X i X j ) X i X j 2 l i K h ( B , X i X l ) l i n K h ( B , X i , X l ) 3 = O p ( h K ) .
Then, the second term of (A2) satisfies d i a g ( 2 s ( Σ n A n k Σ n ) ) / 2 s ( Σ n W n s Σ n ) = O p ( 1 ) . Under the above results, together with B ^ B = O p ( r n ) (refer to [18]) and r n / h 0 , under H 0 , we have
ε W ^ n s ε ε W n s ε 2 s ( Σ n W n s Σ n ) = O p ( r n h ) = o p ( 1 ) .
Step 3: ε ^ W ^ n s ε ^ ε W ^ n s ε 2 s ( Σ n W n s Σ n ) p 0 .
Decompose the quantity into two parts:
ε ^ W ^ n s ε ^ ε W ^ n s ε 2 s ( Σ n W n s Σ n ) = ( ε ^ ε ) W ^ n s ( ε ^ ε ) 2 s ( Σ n W n s Σ n ) + 2 ( ε ^ ε ) W ^ n s ε 2 s ( Σ n W n s Σ n ) .
First, we consider the numerator of the first term on the right side of Equation (A3). It is easy to see that ( ε ^ ε ) W ^ n s ( ε ^ ε ) = ( ε ^ ε ) W n s ( ε ^ ε ) ( 1 + o p ( 1 ) ) ; then,
( ε ^ ε ) W n s ( ε ^ ε ) = 1 2 i = 1 n j i n K h ( B , X i X j ) l i n K h ( B , X i X l ) + K h ( B , X i X j ) l j n K h ( B , X j X l ) ( ε ^ i ε i ) ( ε ^ j ε j ) = 1 2 i = 1 n j i n K h ( B , X i X j ) l i n K h ( B , X i X l ) ( ε ^ i ε i ) ( ε ^ j ε j ) + 1 2 i = 1 n j i n K h ( B , X i X j ) l j n K h ( B , X j X l ) · ( ε ^ i ε i ) ( ε ^ j ε j ) 1 2 D 31 n + 1 2 D 32 n .
Then, since ε ^ ε = g ( β ^ , X ) g ( β , X ) , the Taylor expansion for D 31 n with respect to β leads to
D 31 n = i = 1 n j i n K h ( B , X i X j ) l i n K h ( B , X i X l ) ( ε ^ i ε i ) ( ε ^ j ε j ) = i = 1 n j i n w i j g ( X i , β ) g ( X j , β ) X i , β ^ β X j , β ^ β + o ( β ^ β ) = O p ( n r n 2 ) ,
where β ^ β = O p ( r n ) and K = 1 under H 0 .
Similar to D 31 n , we can prove D 32 n = O P ( n r n 2 ) . Then, since s ( W n s ) = O p ( h 1 / 2 ) (see Step 1), r n / h 0 and n h 1 / 2 r n 2 0 , we have
( ε ^ ε ) W ^ n s ( ε ^ ε ) 2 s ( Σ n W n s Σ n ) p 0 .
Next, we consider the numerator of the second term on the right side of Equation (A3), that is,
( ε ^ ε ) W ^ n s ε = 1 2 i = 1 n j i n K h ( B ^ , X i X j ) l i n K h ( B ^ , X i X l ) + K h ( B ^ , X i X j ) l j n K h ( B ^ , X j X l ) ( ε ^ i ε i ) ε j = 1 2 i = 1 n j i n K h ( B ^ , X i X j ) l i n K h ( B ^ , X i X l ) ( ε ^ i ε i ) ε j + 1 2 i = 1 n j i n K h ( B ^ , X i X j ) l j n K h ( B ^ , X j X l ) ( ε ^ i ε i ) ε j 1 2 D 33 n + 1 2 D 34 n .
Similar to D 31 n , by Taylor expansion of D 33 n with respect to B and β , we have
D 33 n = D 33 n + o p ( D 33 n ) ,
where
D 33 n = i = 1 n j i n w i j ( ε ^ i ε i ) ε j + i = 1 n j i n 1 h s i j X i X j , B ^ B ( ε ^ i ε i ) ε j = i = 1 n j i n w i j g ( X i , β ) ε j X i , β ^ β + i = 1 n j i n 1 h s i j X i X j , B ^ B g ( X i , β ) ε j X i , β ^ β = D 331 n + D 332 n .
Similar to the Lemma 5 in Ait-Saïdi et al. [1], under H 0 with K = 1 , we have
j w i j ε j = C n h + o p ( 1 n h ) + o p ( log ( n ) n h ) + o p ( h b 1 ) .
Since ε ^ ε = O p ( r n ) , B ^ B = O p ( r n ) , β ^ β = O p ( r n ) and s ( W n ) = O p ( h 1 / 2 ) under H 0 , we have
D 331 n 2 s ( Σ n W n s Σ n ) = O p ( r n ) .
Combining the proofs of D 31 n and D 331 n , we obtain
D 332 n = i = 1 n j i n 1 h K h ( B , X i X j ) l i n K h ( B , X i X l ) ε j X i X j , B ^ B g ( X i , β ) X i , β ^ β i = 1 n j i n 1 h l i n K h ( B , X i X l ) l i n K h ( B , X i X l ) w i j ε j X i X j , B ^ B g ( X i , β ) X i , β ^ β = i = 1 n j i n 1 h K h ( B , X i X j ) K h ( B , X i X j ) w i j ε j X i X j , B ^ B g ( X i , β ) X i , β ^ β i = 1 n j i n 1 h l i n K h ( B , X i X l ) l i n K h ( B , X i X l ) w i j ε j X i X j , B ^ B g ( X i , β ) X i , β ^ β = D 3321 n + D 3322 n .
Combining Lemma A1 (i) in [11], and similar to (A6), we have
D 3321 n 2 s ( Σ n W n s Σ n ) = O p ( r n 2 h 2 ) = o p ( 1 ) , D 3322 n 2 s ( Σ n W n s Σ n ) = O p ( r n 2 h ) = o p ( 1 ) .
Then, D 33 n = o ( 1 ) . Similar to D 33 n , it is easy to obtain D 34 n = o p ( 1 ) ; then, we have
2 ( ε ^ ε ) W ^ n s ε 2 s ( Σ n W n s Σ n ) p 0 .
Step 4: s ( U ^ W ^ n s U ^ ) s ( Σ n W n s Σ n ) p 1 .
We prove this in two steps. That is,
i j ( w i j s ) 2 ε i 2 ε j 2 s ( Σ n W n s Σ n ) s 2 ( Σ n W n s Σ n ) p 0 ,
and
i j ( w ^ i j s ) s ε ^ i 2 ε ^ j 2 i j ( w i j s ) s ε i 2 ε j 2 s 2 ( Σ n W n s Σ n ) p 0 ,
From Step 3, we find that we only need to prove that the rate of W n (or W ^ n ) in the formula is the same as the rate of W n s (or W ^ n s ). Then, referring to Lemma 3 and Lemma A.3 in [23], we can easily obtain the conclusion of (A7). Next, combining Step 2 and Step 3, we can obtain (A8).
Step 5: K ^ s ( W ^ n s ) p 0 .
By Taylor expansion of W n s in (A4), and B ^ B = O p ( r n ) , we know that s ( W ^ n s ) = O p ( h K / 2 ) with K = 1 under H 0 . Since P ( K ^ = 1 ) p 1 under H 0 , then
K ^ s ( W ^ n s ) p 0 .
Theorem 1 is proved. □

Appendix A.2. Part II

Proof of Theorem 2. 
Define ε ˜ i = Y i g ( X i , θ ) and ε ˜ = ( ε ˜ 1 , ε ˜ 2 , , ε ˜ n ) . Then, combining the results in Step 2, we have
ε ^ W ^ n s ε ^ 2 s ( W n s ) = ε ˜ W n s ε ˜ 2 s ( W n s ) + ( ε ^ ε ˜ ) W n s ( ε ^ ε ˜ ) 2 s ( W n s ) + 2 ( ε ^ ε ˜ ) W n s ε ˜ 2 s ( W n s ) + o ( 1 ) .
It is easy to show that the first part dominates the third part and goes to infinity in probability and the second part is bounded in probability. The details are similar to the proof of Theorem 1 and are therefore omitted. □

Appendix A.3. Part III

Proof of Proposition 1. 
For the FSIR method, we first define
M s ˙ = ( Γ + ρ I ) 1 V a r [ E ( X | Y I s ˙ ) ] = ( Γ + ρ I ) 1 E [ X I ( Y I s ˙ ) ] p s ˙ E [ X I ( Y I s ˙ ) ] 1 p s ˙ E [ X I ( Y I s ˙ ) ] p s ˙ E [ X I ( Y I s ˙ ) ] 1 p s ˙ p s ˙ ( 1 p s ˙ ) = ( Γ + ρ I ) 1 E [ X I ( Y I s ˙ ) ] E ( X ) E [ I ( Y I s ˙ ) ] p s ˙ ( 1 p s ˙ ) E [ X I ( Y I s ˙ ) ] E ( X ) E [ I ( Y I s ˙ ) ] p s ˙ ( 1 p s ˙ ) p s ˙ ( 1 p s ˙ ) = ( Γ + ρ I ) 1 E { ( X E ( X ) ) I ( Y I s ) } E { ( X E ( X ) ) I ( Y I s ) } / ( p s ˙ ( 1 p s ˙ ) ) ,
where s ˙ is the number of slicing space, which is described in detail in Algorithm 1 of Section 3.1. By Lemma 1 in [14], ρ ( · ) , such that s p a n [ E { M s ˙ } ] = s p a n [ E { M s ˙ ρ ( s ˙ ) } ] , then M s ˙ can also be taken to be
M s ˙ = ( Γ + ρ I ) 1 E { ( X E ( X ) ) I ( Y I s ) } E { ( X E ( X ) ) I ( Y I s ) } .
the sample estimate of M s ˙ is
M n s ˙ = ( Γ ^ + ρ I ) 1 1 n i = 1 n ( X i X ¯ ) I ( Y i I s ) 1 n i = 1 n ( X i X ¯ ) I ( Y i I s ) .
Referring to Theorem 4 in [34], we find that to prove P ( D ( 1 ) > D ( k ) ) 1 under H 1 n for k > 1 , we only need to prove
P n l = 2 k { log ( λ ^ l + 1 ) λ ^ l } 2 l = 1 r { log ( λ ^ l + 1 ) λ ^ l } < C n k ( k + 1 ) 2 r 1 .
The biggest difference between the above formula under H 0 and H 1 n is that the convergence rate about λ ^ is different. To obtain the rate of convergence about λ ^ under H 1 n , we first have to prove M n s ˙ M s ˙ = O p ( C 1 n ) ( C 1 n is defined in model (7)).
Let Y and Y n be the responser under H 0 and H 1 n , respectively. Then, Y n = Y + C 1 n G ( B , X ) . Define I s ˙ = ( a s ˙ , b s ˙ ] , if C 1 n = n 1 / 2 h 1 / 4 , then
1 n i = 1 n X i I ( Y i n I s ˙ ) E [ X I ( Y I s ˙ ) ] = 1 n i = 1 n X i I ( Y i n I s ˙ ) E [ X I ( Y n I s ˙ ) ] + E [ X I ( Y n I s ˙ ) ] E [ X I ( Y I s ˙ ) ] = O p ( n 1 / 2 ) + E X P ( Y n I s ˙ | X ) P ( Y I s ˙ | X ) = E X F Y | X ( b s ˙ C 1 n G ( B , X ) ) F Y | X ( b s ˙ ) F Y | X ( a s ˙ C 1 n G ( B , X ) ) + F Y | X ( a s ˙ ) + O p ( n 1 / 2 ) = E X C 1 n G ( B , X ) { f Y | X ( a s ˙ ) f Y | X ( b s ˙ ) } + O p ( C 1 n ) + O p ( n 1 / 2 ) = O p ( C 1 n ) .
Similarly, M n s ˙ M s ˙ = O p ( C 1 n ) . Then, by Theorem 3 in [45], we have a λ ^ i λ i = O p ( C 1 n ) . Then, l = 2 k { log ( λ ^ l + 1 ) λ ^ l } = O p ( C 1 n 2 ) . To obtain the result of (A9), we need to make n C 1 n 2 / C n 0 hold. When C 1 n = n 1 / 2 h 1 / 4 , since C n h 1 / 2 , we have
n C 1 n 2 C n = 1 C n h 0 .
Therefore, P ( D ( 1 ) > D ( k ) ) 1 under H 1 n . □
Proof of Theorem 3. 
Define δ i = ε i ^ ε i , i = 1 , , n , δ = ( δ 1 , , δ n ) . Similar to the proof of Theorem 1, we have,
T n = ε W n s ε 2 σ s ( W n s ) + 2 δ W n s ε 2 σ s ( W n s ) + δ W n s δ 2 σ s ( W n s ) + o ( 1 ) = D 41 n + D 42 n + D 43 n + o ( 1 ) .
Similar to the proof in Theorem 1, it is easy to obtain D 41 n 2 s ( Σ n W n s Σ n ) d N ( 0 , 1 ) . Similar to (A5), we can prove that D 42 n 2 s ( Σ n W n s Σ n ) = o ( 1 ) .
Lastly, we consider the term D 43 n . Define Z = B , X j as a random variable and the density function exists (Assumption 9). From step 1 in Theorem 1 and Proposition 1, we know that
s 2 ( W n ) = 1 h E [ f 1 ( B , X ) ] K 2 ( u ) d u + o p ( 1 ) .
Similar to the first part of (A3), we have
δ W n δ = i = 1 n j i n w i j [ g ( β ^ , X i ) g ( β 0 , X i ) + C 1 n G ( B , X i ) ] [ g ( β ^ , X j ) g ( β 0 , X j ) + C 1 n G ( B , X j ) ] = i = 1 n j i n w i j [ g ( β ^ , X i ) g ( β 0 , X i ) ] [ g ( β ^ , X j ) g ( β 0 , X j ) ] + 2 C 1 n i = 1 n j i n w i j [ g ( β ^ , X i ) g ( β 0 , X i ) ] G ( B , X j ) + C 1 n 2 i = 1 n j i n w i j G ( B , X i ) G ( B , X j ) = i = 1 n j i n w i j g ( β , X i ) g ( β , X j ) β ^ β 0 , X i β ^ β 0 , X j [ 1 + o p ( 1 ) ] + 2 C 1 n i = 1 n j i n w i j g ( β , X i ) β ^ β 0 , X i G ( B , X j ) + C 1 n 2 i = 1 n j i n w i j G ( B , X i ) G ( B , X j ) Q 1 n + C 1 n Q 2 n + C 1 n 2 Q 3 n .
Similar to Theorem 1, when C 1 n = n 1 / 2 h 1 / 4 , we have
D 43 n = δ W n s δ 2 σ s ( W n s ) = Q 1 n 2 σ s ( W n s ) + C 1 n Q 2 n 2 σ s ( W n s ) + C 1 n 2 Q 3 n 2 σ s ( W n s ) = O p ( n r n 2 h 1 / 2 ) + O p ( C 1 n n r n h 1 / 2 ) + n h 1 / 2 C 1 n 2 E [ G 2 ( B , X ) ] 2 σ 2 E [ f 1 ( B , X ) ] = o p ( 1 ) + o p ( 1 ) + E [ G 2 ( B , X ) ] 2 σ 2 E [ f 1 ( B , X ) ] ,
which leads to
D 43 n p E [ G 2 ( B , X ) ] 2 σ 2 E [ f 1 ( B , X ) ] ω .
This completes the proof of Theorem 3. □

Appendix B

Table A1. Empirical sizes and powers of test statistics under five alternative models in Example 1.
Table A1. Empirical sizes and powers of test statistics under five alternative models in Example 1.
ModelsSamplec S n S n B S n FSIR T n T n B T n FSIR CvM KS
H 11 n = 100 00.0320.0470.0350.0360.0520.0420.0540.059
20.1750.2460.3900.1860.2580.4240.1360.125
40.6990.7710.9410.7490.8140.9570.3270.270
60.9810.9900.9970.9890.9950.9990.5590.471
n = 200 00.0380.0450.0400.0390.0510.0480.0520.047
20.3970.4920.7400.4030.5310.7400.2290.196
40.9810.9880.9990.9890.9910.9990.5820.493
61.0001.0001.0001.0001.0001.0000.8120.723
H 12 n = 100 00.0320.0470.0350.0360.0520.0420.0540.059
20.1700.2170.3900.2040.2500.5570.4040.299
40.7620.7420.5900.8200.7900.6500.7610.673
60.9780.9570.8770.9830.9680.8930.8930.839
n = 200 00.0380.0450.0400.0390.0510.0480.0520.047
20.4260.4900.7980.4910.5690.9350.6910.556
40.9890.9900.9310.9970.9960.9510.9720.938
61.0001.0000.9971.0001.0000.9970.9930.989
H 13 n = 100 00.0320.0470.0350.0360.0520.0420.0540.059
20.0900.1410.1230.1150.1530.1980.2570.181
40.4530.4940.5790.5280.5500.6460.6340.522
60.8190.7980.8440.8770.8480.8790.7880.710
n = 200 00.0380.0450.0400.0390.0510.0480.0520.047
20.1870.2650.2630.2700.3450.4340.5030.368
40.8480.8770.9170.9280.9320.9440.8830.818
60.9920.9930.9950.9990.9991.0000.9540.926
H 14 n = 100 00.0320.0470.0350.0360.0520.0420.0540.059
20.1170.1660.1510.1460.2010.2330.1950.158
40.5230.5470.7580.6130.6240.7990.3240.279
60.8530.8430.9500.8720.8830.9540.3710.326
n = 200 00.0380.0450.0400.0390.0510.0480.0520.047
20.2910.3750.3520.3990.4860.5210.2850.235
40.8350.8480.9770.8760.8760.9850.4600.414
61.0001.0001.0001.0001.0000.9990.5090.467
H 15 n = 100 00.0320.0470.0350.0360.0520.0420.0540.059
20.2440.3140.3480.2830.3500.4440.3490.278
40.8260.8360.9370.8570.8580.9560.5410.496
60.8810.8760.9960.8830.8670.9950.5700.520
n = 200 00.0380.0450.0400.0390.0510.0480.0520.047
20.5810.6610.7040.6790.7450.8320.5170.454
40.8960.8981.0000.8980.8990.9990.7170.671
61.0001.0001.0001.0001.0001.0000.7260.687
H 16 n = 100 00.0320.0470.0350.0360.0520.0420.0540.059
20.2420.3150.3430.2810.3500.4790.3440.296
40.7300.7240.9460.7550.7430.9580.5300.469
60.8900.8800.9950.8900.8820.9960.5790.528
n = 200 00.0380.0450.0400.0390.0510.0480.0520.047
20.5820.6610.7370.6760.7340.8570.5290.466
40.8980.8971.0000.9010.9021.0000.7170.670
61.0001.0001.0001.0001.0001.0000.7400.703
Table A2. Empirical sizes and powers of test statistics under four DGM in Example 2.
Table A2. Empirical sizes and powers of test statistics under four DGM in Example 2.
DGMSamplec S n S n B S n FSIR T n T n B T n FSIR CvM KS
D G M 1 n = 100 00.0370.0450.0430.0390.0440.0460.0530.049
20.1680.2420.2580.1850.2540.3700.2620.192
40.6670.7030.9070.7160.7430.9470.4170.369
60.8480.8430.9900.8650.8510.9940.5440.486
n = 200 00.0410.0540.0380.0470.0410.0530.0570.059
20.4010.4900.5600.4420.5620.7500.3960.319
40.8850.8841.0000.8960.8951.0000.6490.601
61.0001.0001.0001.0001.0001.0000.6980.646
D G M 2 n = 100 00.0390.0490.0520.0380.0480.0560.0500.054
20.2380.2980.3590.2620.3100.4910.3420.266
40.8060.8140.9660.8380.8310.9850.5630.485
60.8750.8700.9990.8840.8740.9990.6660.586
n = 200 00.0330.0610.0570.0390.0460.0520.0620.049
20.5340.6140.7300.6090.6920.8670.5380.443
40.8970.8991.0000.8970.8981.0000.7250.661
61.0001.0001.0001.0001.0001.0000.8050.773
D G M 3 n = 100 00.0300.0580.0440.0400.0630.0520.0610.066
20.0590.0970.0600.0590.1000.0580.1130.100
40.1750.2510.2840.1850.2520.3920.2110.177
60.4300.5060.7940.4790.5420.8820.3370.274
n = 200 00.0360.0470.0450.0380.0460.0480.0540.053
20.1000.1650.1490.1180.1920.1660.1860.165
40.4680.5580.6830.5300.6200.8180.3520.292
60.7780.7860.9920.8120.8330.9970.5010.424
D G M 4 n = 100 00.0380.0540.0420.0390.0530.0550.0530.061
20.2250.3020.2660.2420.3190.2450.3150.284
40.7810.7850.9530.8190.8100.9750.5640.506
60.8830.8720.9970.8850.8750.9980.6460.591
n = 200 00.0390.0420.0410.0430.0450.0510.0610.071
20.5270.6230.5860.6190.6840.6200.5130.456
40.9960.9971.0000.9970.9971.0000.7470.686
61.0001.0001.0001.0001.0001.0000.8160.776
Table A3. The p-value of different test statistics in real data analysis.
Table A3. The p-value of different test statistics in real data analysis.
Statistic S n S n B S n FSIR T n T n B T n FSIR CvM KS
p-value0.0000.0090.0000.0000.0130.0000.5280.374
Figure A1. Empirical sizes and powers of T n , T n B and T n FSIR against the value of c at 5 % significance level in Examples 3 and 4.
Figure A1. Empirical sizes and powers of T n , T n B and T n FSIR against the value of c at 5 % significance level in Examples 3 and 4.
Mathematics 11 01812 g0a1
Figure A2. COVID-19 data in real data analysis.
Figure A2. COVID-19 data in real data analysis.
Mathematics 11 01812 g0a2

References

  1. Ait-Saïdi, A.; Ferraty, F.; Kassa, R.; Vieu, P. Cross-validated estimations in the single-functional index model. Statistics 2008, 42, 475–494. [Google Scholar] [CrossRef]
  2. Ma, S. Estimation and inference in functional single-index models. Ann. Inst. Stat. Math. 2016, 68, 181–208. [Google Scholar] [CrossRef]
  3. Wang, G.; Feng, X.N.; Chen, M. Functional partial linear single-index model. Scand. J. Stat. 2016, 43, 261–274. [Google Scholar] [CrossRef]
  4. Jiang, F.; Baek, S.; Cao, J.; Ma, Y. A Functional single-index model. Stat. Sin. 2020, 30, 303–324. [Google Scholar] [CrossRef]
  5. Jäntschi, L. A test detecting the outliers for continuous distributions based on the cumulative distribution function of the data being tested. Symmetry 2019, 11, 835. [Google Scholar] [CrossRef] [Green Version]
  6. Cardot, H.; Ferraty, F.; Mas, A.; Sarda, P. Testing hypotheses in the functional linear model. Scand. J. Stat. 2003, 30, 241–255. [Google Scholar] [CrossRef]
  7. Kokoszka, P.; Maslova, I.; Sojka, J.; Zhu, L. Testing for lack of dependence in the functional linear model. Can. J. Stat. 2008, 36, 207–222. [Google Scholar] [CrossRef]
  8. Hilgert, N.; Mas, A.; Verzelen, N. Minimax adaptive tests for the functional linear model. Ann. Stat. 2013, 41, 838–869. [Google Scholar] [CrossRef]
  9. Patilea, V.; Sánchez-Sellero, C.; Saumard, M. Projection-based nonparametric goodness-of-fit testing with functional covariates. arXiv 2012, arXiv:1205.5578. [Google Scholar]
  10. Patilea, V.; Sánchez-Sellero, C.; Saumard, M. Testing the predictor effect on a functional response. J. Am. Stat. Assoc. 2016, 111, 1684–1695. [Google Scholar] [CrossRef]
  11. Shi, G.; Du, J.; Sun, Z.; Zhang, Z. Checking the adequacy of functional linear quantile regression model. J. Stat. Plan. Inference 2021, 210, 64–75. [Google Scholar] [CrossRef]
  12. García-Portugués, E.; González-Manteiga, W.; Febrero-Bande, M. A goodness-of-fit test for the functional linear model with scalar response. J. Comput. Graph. Stat. 2014, 23, 761–778. [Google Scholar] [CrossRef] [Green Version]
  13. Cuesta-Albertos, J.A.; García-Portugués, E.; Febrero-Bande, M.; González-Manteiga, W. Goodness-of-fit tests for the functional linear model based on randomly projected empirical processes. Ann. Stat. 2019, 47, 439–467. [Google Scholar] [CrossRef] [Green Version]
  14. Guo, X.; Wang, T.; Zhu, L. Model checking for parametric single-index models: A dimension reduction model-adaptive approach. J. R. Stat. Soc. Ser. B Stat. Methodol. 2016, 78, 1013–1035. [Google Scholar] [CrossRef]
  15. Stute, W.; Thies, S.; Zhu, L.X. Model checks for regression: An innovation process approach. Ann. Stat. 1998, 26, 1916–1934. [Google Scholar] [CrossRef]
  16. Hardle, W.; Mammen, E. Comparing nonparametric versus parametric regression fits. Ann. Stat. 1993, 21, 1926–1947. [Google Scholar] [CrossRef]
  17. Patilea, V.; Sánchez-Sellero, C. Testing for lack-of-fit in functional regression models against general alternatives. J. Stat. Plan. Inference 2020, 209, 229–251. [Google Scholar] [CrossRef]
  18. Lian, H. Functional sufficient dimension reduction: Convergence rates and multiple functional case. J. Stat. Plan. Inference 2015, 167, 58–68. [Google Scholar] [CrossRef]
  19. Ferré, L.; Yao, A.F. Functional sliced inverse regression analysis. Statistics 2003, 37, 475–488. [Google Scholar] [CrossRef]
  20. Ferré, L.; Yao, A.F. Smoothed functional inverse regression. Stat. Sin. 2005, 15, 665–683. [Google Scholar]
  21. Wang, G.; Lin, N.; Zhang, B. Functional k-means inverse regression. Comput. Stat. Data Anal. 2014, 70, 172–182. [Google Scholar] [CrossRef]
  22. Lian, H.; Li, G. Series expansion for functional sufficient dimension reduction. J. Multivar. Anal. 2014, 124, 150–165. [Google Scholar] [CrossRef]
  23. Ellison, G.; Ellison, S.F. A simple framework for nonparametric specification testing. J. Econom. 2000, 96, 1–23. [Google Scholar] [CrossRef] [Green Version]
  24. Stute, W.; Zhu, L.X. Model checks for generalized linear models. Scand. J. Stat. 2002, 29, 535–545. [Google Scholar] [CrossRef]
  25. Hall, P.; Horowitz, J.L. Methodology and convergence rates for functional linear regression. Ann. Stat. 2007, 35, 70–91. [Google Scholar] [CrossRef] [Green Version]
  26. Yuan, M.; Cai, T.T. A reproducing kernel Hilbert space approach to functional linear regression. Ann. Stat. 2010, 38, 3412–3444. [Google Scholar] [CrossRef]
  27. Crambes, C.; Kneip, A.; Sarda, P. Smoothing splines estimators for functional linear regression. Ann. Stat. 2009, 37, 35–72. [Google Scholar] [CrossRef]
  28. Zhao, Y.; Ogden, R.T.; Reiss, P.T. Wavelet-based LASSO in functional linear regression. J. Comput. Graph. Stat. 2012, 21, 600–617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Li, B. Sufficient Dimension Reduction: Methods and Applications with R; Chapman and Hall/CRC: Boca Raton, FL, USA, 2018. [Google Scholar]
  30. Xu, J.; Cheng, H.; Cui, W.; Li, Y. Sufficient dimension reduction via distance covariance for functional and longitudinal data. arXiv 2022, arXiv:2202.13579. [Google Scholar]
  31. Li, K.C. Sliced inverse regression for dimension reduction. J. Am. Stat. Assoc. 1991, 86, 316–327. [Google Scholar] [CrossRef]
  32. Lin, Q.; Zhao, Z.; Liu, J.S. On consistency and sparsity for sliced inverse regression in high dimensions. Ann. Stat. 2018, 46, 580–610. [Google Scholar] [CrossRef] [Green Version]
  33. Zhu, L.; Miao, B.; Peng, H. On sliced inverse regression with high-dimensional covariates. J. Am. Stat. Assoc. 2006, 101, 630–643. [Google Scholar] [CrossRef]
  34. Zhu, L.; Wang, T.; Zhu, L.; Ferré, L. Sufficient dimension reduction through discretization-expectation estimation. Biometrika 2010, 97, 295–304. [Google Scholar] [CrossRef]
  35. Zhu, L.P.; Zhu, L.X.; Feng, Z.H. Dimension reduction in regressions through cumulative slicing estimation. J. Am. Stat. Assoc. 2010, 105, 1455–1466. [Google Scholar] [CrossRef]
  36. Sheng, W.; Yin, X. Sufficient dimension reduction via distance covariance. J. Comput. Graph. Stat. 2016, 25, 91–104. [Google Scholar] [CrossRef]
  37. Zhang, J.; Chen, X. Robust sufficient dimension reduction via ball covariance. Comput. Stat. Data Anal. 2019, 140, 144–154. [Google Scholar] [CrossRef]
  38. Du, P.; Wang, X. Penalized likelihood functional regression. Stat. Sin. 2014, 24, 1017–1041. [Google Scholar] [CrossRef] [Green Version]
  39. Cardot, H.; Sarda, P. Estimation in generalized linear models for functional data via penalized likelihood. J. Multivar. Anal. 2005, 92, 24–41. [Google Scholar] [CrossRef] [Green Version]
  40. Zheng, J.X. A consistent nonparametric test of parametric regression models under conditional quantile restrictions. Econom. Theory 1998, 14, 123–138. [Google Scholar] [CrossRef]
  41. Niu, C.; Zhu, L. An adaptive-to-model test for parametric single-index models with missing responses. Electron. J. Stat. 2017, 11, 1491–1526. [Google Scholar] [CrossRef]
  42. Zhu, X.; Guo, X.; Zhu, L. An adaptive-to-model test for partially parametric single-index models. Stat. Comput. 2017, 27, 1193–1204. [Google Scholar] [CrossRef] [Green Version]
  43. Shi, E.; Liu, Y.; Sun, K.; Li, L.; Kong, L. An adaptive model checking test for functional linear model. arXiv 2022, arXiv:2204.01831. [Google Scholar]
  44. Mikosch, T. Functional limit theorems for random quadratic forms. Stoch. Process. Their Appl. 1991, 37, 81–98. [Google Scholar] [CrossRef] [Green Version]
  45. Li, B.; Wen, S.; Zhu, L. On a projective resampling method for dimension reduction with multivariate responses. J. Am. Stat. Assoc. 2008, 103, 1177–1186. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xia, L.; Lai, T.; Zhang, Z. An Adaptive-to-Model Test for Parametric Functional Single-Index Model. Mathematics 2023, 11, 1812. https://doi.org/10.3390/math11081812

AMA Style

Xia L, Lai T, Zhang Z. An Adaptive-to-Model Test for Parametric Functional Single-Index Model. Mathematics. 2023; 11(8):1812. https://doi.org/10.3390/math11081812

Chicago/Turabian Style

Xia, Lili, Tingyu Lai, and Zhongzhan Zhang. 2023. "An Adaptive-to-Model Test for Parametric Functional Single-Index Model" Mathematics 11, no. 8: 1812. https://doi.org/10.3390/math11081812

APA Style

Xia, L., Lai, T., & Zhang, Z. (2023). An Adaptive-to-Model Test for Parametric Functional Single-Index Model. Mathematics, 11(8), 1812. https://doi.org/10.3390/math11081812

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop