Next Article in Journal
Positivity-Preserving Rational Cubic Fractal Interpolation Function Together with Its Zipper Form
Next Article in Special Issue
Spatio-Functional Nadaraya–Watson Estimator of the Expectile Shortfall Regression
Previous Article in Journal
Extending L-Topologies to Bipolar L-Fuzzy Topologies
Previous Article in Special Issue
On Modulus Statistical Convergence in Partial Metric Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generalized Partially Functional Linear Model with Interaction between Functional Predictors

1
School of Science, North China University of Technology, Beijing 100144, China
2
Department of Statistics, University of Leeds, Leeds LS2 9JT, UK
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(9), 583; https://doi.org/10.3390/axioms13090583
Submission received: 10 July 2024 / Revised: 15 August 2024 / Accepted: 19 August 2024 / Published: 27 August 2024
(This article belongs to the Special Issue Advances in Functional and Topological Data Analysis)

Abstract

:
This paper proposes a generalized partially functional linear model with interaction terms. It is suitable for cases where the response variable is scalar, and the predictor variables include a mix of functional and scalar types, while considering the correlations among functional predictor variables. The model uses principal component analysis for dimensionality reduction, employs maximum likelihood estimation to obtain parameter values, proves the asymptotic properties of the estimates, and validates the model’s accuracy through data simulation experiments. Finally, the proposed model was applied to investigate the influence of air quality, climate factors, and medical and social indicators, along with their interactions, on cancer incidence, which is a binary response.

1. Introduction

With the advent of the big data era, more and more functional data, providing information about objects varying over a continuum, are collected.
Currently, functional data analysis is being applied in various fields such as medicine, environmental science, and economics and is receiving increasing attention. For details on functional data analysis, see monographs by Ramsay and Silverman [1], Horváth and Kokoszka [2], and Hsing and Eubank [3].
Several variants of functional linear regression models have been proposed to investigate the influence of functional and/or scalar predictors on functional or scalar response and, therefore, to make predictions. Cardot [4], Tony [5], and others have utilized spline methods for estimation and prediction in functional linear regression models. In 2007, Cardot et al. [6] extended the population least squares method to functional linear models, proposing smooth spline estimates for model function coefficients and providing asymptotic results for this estimation. In 2012, Delaigle and Hall [7] utilized partial least squares to demonstrate the consistency and convergence of functional linear models. Tony and Ming [8] studied the estimation and prediction issues of functional linear regression models within the framework of reproducing kernel Hilbert spaces. Nevertheless, these models cannot deal with general responses such as binary and Poisson.
An important tool for functional data analysis is the functional linear regression model, while the generalized functional regression model is an extension of the functional linear regression model. As research progressed, the generalized functional regression model was introduced to handle more complex response variables. This model was first introduced by Nelder and Wedderburn [9] in 1972, and it investigates the relationship between continuous and discrete response variables and the predictor variables through a link function. In 2002, James [10] proposed generalized linear models with functional predictors and applied it to standard missing data problems. In 2005, Müller and Stadtmüller [11] proposed a generalized functional linear regression model where the response variable is a discrete scalar and the predictor is functional. In 2011, Goldsmith et al. [12] developed fast fitting methods for generalized functional linear models that can be applied to various functional data designs, including functions measured with and without error and sparsely or densely sampled. In 2021, Xiao et al. [13] proposed a generalized partially functional linear regression model where the response variable is general and the predictors are scalar and functional. However, none of these models incorporate the interaction of functional predictors.
To better address complex data that include both functional predictors and scalar predictors, scholars have improved the functional linear model and proposed a functional regression model with mixed predictors. In 2016, Kong et al. [14] explored the estimation and variable selection problems in cases where the parametric part is high-dimensional and the functional predictors are multidimensional. Yao [15] and Ma et al. [16] further built upon the work of Kong [14], conducting more in-depth research and proving the large sample properties of the estimators. In 2020, Xu et al. [17] studied the estimation and hypothesis testing issues for models with multiple functional predictors and demonstrated the corresponding large sample properties.
In many practical applications, we need to consider the interactions between variables, and failure to consider the interaction term may lead to the problem of missing variables in the model, thus introducing inaccurate predictions and inappropriate interpretations. By introducing interaction terms, the inaccuracy can be reduced and the model can be made more reliable, thereby improving the prediction by the model and providing more reliable decision support. Indeed, functional linear regressions models with interaction between functional predictors have been proposed recently; several examples follow. In 2016, Usset et al. [18] proposed a functional regression model with a scalar response and multiple functional predictors with two-way interactions in addition to their main effects. In 2019, Luo and Qin [19] proposed function-on-function regression models with interaction and quadratic effects, together with an efficient estimation method that has a minimum prediction error. In 2013, Yang et al. [20] introduced a class of nonlinear multivariate time-frequency functional models that can identify important features in each signal as well as the interaction of signals. Some models considered the interaction of two different time points in the functional data. In 2020, Matsui [21] proposed a functional quadratic model which took the interaction between two different time points of the functional data into consideration. In 2020, Sun and Wang [22] also considered a quadratic regression model where the predictor and the response are both functional; it estimated predictions for the coefficient functions, and unknown responses and asymptotics were demonstrated. Nonetheless, these models cannot be applied to general scalar responses. As far as we know, only Fuchs et al. [23] in 2015 considered general scalar response with functional predictors to include linear functional interaction terms. However, one drawback of the method of Fuchs et al. [23] is that scalar predictors are not included, and a second drawback is that the asymptotic properties of estimated regression coefficients were not established.
A practical motivation of this paper is the investigation of the influence of air qualities, climate factors, medical and social indicators, and their interactions on cancer incidence, which is a binary response. Cancer is one of the leading causes of death in humans; therefore, it is crucial to analyze the factors related to cancer incidence. Studying cancer incidence can help improve public health and quality of life, reduce social medical costs, and promote human health and socio-economic development. In 2022, Qiu et al. [24] pointed out that cancer incidence in China is much higher than those in the United States and the United Kingdom due to the fact that China faces problems such as a large population, uneven development in various regions, and a relative lag in cancer control strategies. In 2014, Qin et al. [25] indicated that long-term exposure to air pollutants or short-term exposure to some high concentrations of air pollutants such as PM2.5 may be associated with some increased incidence rates of overall cancer, especially prostate cancer and female breast cancer. In 2022, Wu et al. [26] found that areas with high green coverage have a lower risk of cancer. In 2023, Cao et al. [27] analyzed the relationship between per capita GDP and cancer incidence in 55 regions of China, showing that regions with high GDP have high cancer incidence. In 2017, Xu et al. [28] conducted a statistical analysis of the current situation of PM2.5 in Changzhou in China and considered an interaction between PM2.5 and relative humidity during the same period, indicating a certain degree of interaction between the two. In 2022, Yang et al. [29] used the generalized linear model to study the effects of PM2.5 and relative humidity on visibility and found a significant interaction between PM2.5 and relative humidity.
Therefore, we collected data on average daily PM2.5 concentration (from 1 January 2015 to 31 December 2020), average daily humidity (from 1 January 2015 to 31 December 2020), per capita GDP, green coverage rate in built-up areas, the proportion of medical personnel (PMP) (which is the ratio of the number of licensed (assistant) doctors to the population in the locality), and the binary cancer incidence in 49 cities in China from http://www.cnemc.cn/, http://www.stats.gov.cn/sj/ndsj/ and http://www.chinancpcn.org.cn/home. Our aim was to investigate the influence of PM2.5 concentration, air humidity, per capita GDP, green coverage, and PMP on cancer incidence, with the focus not only on the main effects but also on the interaction between PM2.5 concentration and air humidity to, therefore, make predictions.
Existing models with interaction terms between functional predictors and general scalar responses cannot deal with multiple functional and scalar predictors, which is the case in our motivated datasets. Moreover, the asymptotic properties of estimators have not been addressed in existing models. Therefore, in Section 2, we fully consider the combined influence of functional predictors, scalar predictors, and interactions between functional predictors on general scalar response by proposing a generalized partially functional linear model with interaction terms. In Section 3, the asymptotic properties of our proposed estimators are established. Extensive simulation studies are given in Section 4. Section 5 is reserved for the real data analysis.

2. Model and Estimation

2.1. Model Introduction

Suppose we have n subjects, and the data we observe for the i-th subject are { ( X i 1 ( t 1 ) , t 1 T 1 ) , ( X i 2 ( t 2 ) , t 2 T 2 ) , Z i , Y i } , i = 1 , , n . For j = 1 , 2 , the functional predictor X i j ( t j ) is a random curve, which is observed for subject i and X i j ( t j ) L 2 ( T j ) , where T j is a bounded interval of R . Notice that, for the sake of simplicity in notations, we only consider the case with two functional predictors, and the case with multiple functional predictors can be easily similarly established. The scalar predictor vector Z = ( Z 1 , Z 2 , , Z q ) T is a q-dimensional random vector. The response Y i is a real-valued random variable that may be continuous or discrete (e.g., binary, count, etc.).
We assume that there is a known link function g ( · ) , which is a monotone and twice continuously differentiable function with bounded derivatives that is, thus, invertible.
We introduce the following generalized partially functional linear model with interaction between the functional predictors:
Y i = g α + T 1 X i 1 t 1 β 1 t 1 d t 1 + T 2 X i 2 t 2 β 2 t 2 d t 2 + T 1 × T 2 X i 1 t 1 X i 2 t 2 β t 1 , t 2 d t 1 d t 2 + Z i T γ + ε i ,
where α R is the intercept, β 1 ( t 1 ) , β 2 ( t 2 ) , and β ( t 1 , t 2 ) are the regression coefficient functions corresponding to the two functional predictors and the interaction term, respectively, Z i represents the i-th random variable in Z, and γ = ( γ 1 , γ 2 , , γ q ) T is the regression coefficient corresponding to the multiple scalar predictors Z. It is assumed that ε i has mean 0 and variance σ 2 and that ε i is independent with ε j if i j .
Define the linear operator :
= α + T 1 X 1 t 1 β 1 t 1 d t 1 + T 2 X 2 t 2 β 2 t 2 d t 2 + T 1 × T 2 X 1 t 1 X 2 t 2 β t 1 , t 2 d t 1 d t 2 + Z T γ .
We specify
E ( Y | X 1 ( · ) , X 2 ( · ) , Z ) = η = g ( ) , V a r ( Y | X 1 ( · ) , X 2 ( · ) , Z ) = σ 2 ( η ) .
For simplicity, we assume that the predictors X j ( t j ) and Z are both centralized, i.e., E ( X j ( t j ) ) = 0 , j = 1 , 2 and E ( Z l ) = 0 , l = 1 , , q . The Karhunen–Loève expansion is a mathematical technique for representing stochastic processes in terms of an orthonormal basis derived from the process’s covariance function. It decomposes a random process into a series of orthogonal functions, each weighted by uncorrelated random coefficients. The basis functions are the eigenfunctions of the covariance function, and the expansion efficiently captures the process’s variability, often using only a few terms. For a detailed introduction to the Karhunen–Loève expansion, please refer to Equation (2.8) in [2]. Based on Karhunen–Loève expansion, X i j ( t j ) can be expanded as
X i 1 ( t 1 ) = k = 1 χ i 1 k φ 1 k ( t 1 ) ,
X i 2 ( t 2 ) = l = 1 χ i 2 l φ 2 l ( t 2 ) ,
where χ i 1 k , χ i 2 l are the functional principal component scores, φ 1 k ( t 1 ) , φ 2 l ( t 2 ) are the functional principal component bases, and
T 1 φ 1 k 2 ( t 1 ) d t 1 = 1 , T 2 φ 2 l 2 ( t 2 ) d t 2 = 1 .
Using the functional principal component bases, the regression coefficient functions β j ( t j ) , β ( t 1 , t 2 ) are expanded as
β 1 ( t 1 ) = k = 1 b 1 k φ 1 k ( t 1 ) ,
β 2 ( t 2 ) = l = 1 b 2 l φ 2 l ( t 2 ) ,
β ( t 1 , t 2 ) = k = 1 l = 1 u k l φ 1 k ( t 1 ) φ 2 l ( t 2 ) .
Plugging the above expansions into Model (1) and truncating the predictors at p j , which increases asymptotically with n , we can obtain truncated Model (2):
Y i = g ( α + k = 1 p 1 χ i 1 k b 1 k + l = 1 p 2 χ i 2 l b 2 l + k = 1 p 1 l = 1 p 2 ρ i k l u k l + Z i T γ ) + ε i ,
where ρ i k l = χ i 1 k · χ i 2 l .

2.2. Parameter Estimation

Define the parameter vector
Ω = b 11 , , b 1 p 1 , b 21 , , b 2 p 2 , u 11 , , u 1 p 2 , u 21 , , u 2 p 2 , , u p 1 1 , , u p 1 p 2 , γ 0 , γ 1 , , γ q T ,
and define
i = α + k = 1 p 1 χ i 1 k b 1 k + l = 1 p 2 χ i 2 l b 2 l + ρ i T u + Z i T γ ,
η i = g ( i ) ,
ω i = χ i 11 , , χ i 1 p 1 , χ i 21 , , χ i 2 p 2 , ρ i 11 , , ρ i 1 p 2 , ρ i 21 , , ρ i 2 p 2 , , ρ i p 1 1 , , ρ i p 1 p 2 , z i 0 , z i 1 , , z i q T ,
where
b j = ( b j 1 , , b j p j ) T , j = 1 , 2 ,
u = u 11 , , u 1 p 2 , u 21 , , u 2 p 2 , , u p 1 1 , , u p 1 p 2 T ,
γ = ( γ 0 , γ 1 , , γ q ) T ,
ρ i = ( ρ i 11 , , ρ i 1 p 2 , ρ i 21 , , ρ i 2 p 2 , , ρ i p 1 1 , , ρ i p 1 p 2 ) T ,
z i 0 = 1 , and γ 0 = α .
The maximum likelihood estimate Ω ^ of Ω can be obtained by solving Equation (3):
U ( Ω ) = i = 1 n ( Y i g ( i ) ) g ( i ) σ 2 ( η i ) ω i = 0 ,
Ω ^ = b ^ 11 , , b ^ 1 p 1 , b ^ 21 , , b ^ 2 p 2 , u ^ 11 , , u ^ 1 p 2 , u ^ 21 , , u ^ 2 p 2 , , u ^ p 1 1 , , u ^ p 1 p 2 , γ ^ 0 , γ ^ 1 , , γ ^ q T ,
where
b ^ j = ( b ^ j 1 , , b ^ j p j ) T , j = 1 , 2 ,
u ^ = u ^ 11 , , u ^ 1 p 2 , u ^ 21 , , u ^ 2 p 2 , , u ^ p 1 1 , , u ^ p 1 p 2 T ,
α ^ = γ ^ 0 , and γ ^ = ( γ ^ 0 , γ ^ 1 , , γ ^ q ) T are the estimates of b j , u , α , γ , respectively.
Introducing the following matrices:
V = d i a g ( σ 2 ( η 1 ) , , σ 2 ( η n ) ) ,
W = d i a g ( g ( 1 ) , g ( 2 ) , , g ( n ) ) ,
A 0 = A n , q + 1 = g ( i ) z i m σ ( η i ) 1 i n , 0 m q ,
A j = A n , p j = g ( i ) χ i j r σ ( η i ) 1 i n , 0 r p j , j { 1 , 2 } ,
A 12 = A n , t = g ( i ) ρ i t σ ( η i ) 1 i n , 1 t p 1 p 2 ,
A = A n , q + 1 + p 1 + p 2 + p 1 p 2 = diag A 1 , A 2 , A 12 , A 0 ,
and the vectors Y = ( Y 1 , , Y n ) T ,   η = ( η 1 , , η n ) T , then Equation (3) can be written as
A T V 1 2 ( Y η ) = 0 .
The estimation of Ω is usually solved iteratively using a weighted least squares method. By Taylor expansion, we have
g 1 ( Y ) = g 1 ( η ) + [ g 1 ( η ) ] ( Y η ) = + W 1 ( Y η ) ;
thus, there is
A T H ( g 1 ( Y ) ) = 0 ,
where H = V 1 2 W .
Simplify to obtain estimates of b j , γ , u :
b ˜ j = ( A j T A j ) 1 A j T H g 1 ( Y ) ,
γ ˜ = A 0 T A 0 1 A 0 T H g 1 ( Y ) ,
u ˜ = A 12 T A 12 1 A 12 T H g 1 ( Y ) .
Repeat the above process until convergence; then, the estimate of Ω is obtained:
Ω ^ = b ^ 11 , , b ^ 1 p 1 , b ^ 21 , , b ^ 2 p 2 , u ^ 11 , , u ^ 1 p 2 , u ^ 21 , , u ^ 2 p 2 , , u ^ p 1 1 , , u ^ p 1 p 2 , γ ^ 0 , γ ^ 1 , , γ ^ q T .

3. Asymptotic Properties

Considering the truncated Model (2), we have the metric
d G 2 ( β ^ j , β j ) = ( b ^ j b j ) T Γ ˜ j ( b ^ j b j ) + k 1 , k 2 = p j + 1 λ j , k 1 k 2 b ¯ j 2 , j = 1 , 2 ,
where
b j = ( b j 1 , b j 2 , , b j p j ) ,
Γ ˜ j = λ j , k 1 k 2 1 k 1 , k 2 p j
is a symmetric positive definite matrix and
λ j , k 1 k 2 = E g ( ) 2 σ 2 ( η ) χ j k 1 χ j k 2 1 k 1 , k 2 p j
is an eigenvalue of the generalized self-covariance operator A G j with kernel
G j ( s , t ) = E [ g ( ) 2 σ 2 ( η ) X j ( s ) X j ( t ) ] ,
and we have
Γ ˜ j 1 = ( ξ j , k 1 k 2 ) 1 k 1 , k 2 p j
and
b ¯ j = ( b j ( p j + 1 ) , b j ( p j + 2 ) , ) T .
Combined with Corollary 4.1 in Müller (2005) [11], we have
k 1 , k 2 = p j + 1 λ j , k 1 k 2 b ¯ j 2 = o p j n .
We specify
f 2 = S f ( s ) 2 d s , f L 2 ( S ) ,
g 2 = S T g ( s , t ) 2 d s d t , g L 2 ( S × T )
and
( f g ) ( x , y ) = f ( x ) g ( y ) , x X , y Y ,
where X , Y are the domains of f , g , respectively.
Define C X j as the covariance function of a random function X j , for j = 1 , 2 . By Mercer’s theorem,
C X 1 ( t 11 , t 12 ) = k 1 λ k φ 1 k ( t 11 ) φ 1 k ( t 12 ) ,
C X 2 ( t 21 , t 22 ) = l 1 σ l φ 2 l ( t 21 ) φ 2 l ( t 22 ) ,
where t 11 , t 12 T 1 ; t 21 , t 22 T 2 ; λ k and φ 1 k , k = 1 , 2 , , are the non-negative eigenvalues and the corresponding eigenfunctions of the covariance function C X 1 ( t 11 , t 12 ) ; and σ l and φ 2 l , l = 1 , 2 , , are the non-negative eigenvalues and the corresponding eigenfunctions of the covariance function C X 2 ( t 21 , t 22 ) .
In order to derive the asymptotic nature of the regression coefficients, we have made the following assumptions in addition to the basic conditions in Section 2:
(i)
The connected function g ( · ) is monotonically invertible and has bounded second-order derivatives, the derivative of the variance function σ 2 ( · ) is continuously bounded, and there exists an σ ( · ) > Δ > 0 ;
(ii)
The scalar predictor variable Z and the functional predictor variable X j ( t j ) are independent of each other;
(iii)
When n , p j satisfies p j and p j n 1 4 0 ;
(iv)
E [ T j { X j ( t j ) } 4 d t j ] < ;
(v)
Define μ X 1 , k = min 1 k p 1 ( λ p 1 λ k + 1 ) , μ X 2 , l = min 1 l p 2 ( σ l σ l + 1 ) , and μ X 1 , k > 0 , μ X 2 , l > 0 ;
(vi)
Define d n = C ^ X C X , K ˜ n = min k 1 : λ k 2 d n 1 , L ˜ n = min l 1 : σ l 2 d n 1 ; d n 0 , K ˜ n , and L ˜ n when n .
Lemma 1.
If the above basic conditions and assumptions hold, while p 1 K ˜ n , p 2 L ˜ n and
k = 1 p 1 l = 1 p 2 1 n ( 1 μ X 1 , k 2 + 1 μ X 2 , l 2 ) 0 ,
we have β ^ β 2 = O p t = 1 p 1 p 2 ( u t u ^ t ) 2 .
Proof. 
β β ^ 2 = k = 1 l = 1 u k l φ 1 k φ 2 l k = 1 p 1 l = 1 p 2 u ^ k l φ ^ 1 k φ ^ 2 l 2 = k = 1 p 1 l = 1 p 2 u k l ( φ 1 k φ 2 l φ ^ 1 k φ ^ 2 l ) + k = 1 p 1 l = 1 p 2 ( u k l u ^ k l ) φ ^ 1 k φ ^ 2 l + k = 1 p 1 l > p 2 u k l φ 1 k φ 2 l + k > p 1 l = 1 u k l φ 1 k φ 2 l 2 4 k = 1 p 1 l = 1 p 2 u k l ( φ 1 k φ 2 l φ ^ 1 k φ ^ 2 l ) 2 + 4 k = 1 p 1 l = 1 p 2 ( u k l u ^ k l ) 2 + 4 k = 1 p 1 l > p 2 u k l 2 + 4 k > p 1 l = 1 u k l 2 = 4 I 1 + 4 t = 1 p 1 p 2 ( u t u ^ t ) 2 + 4 R β ( p 1 , p 2 ) 2 ,
where R β ( p 1 , p 2 ) = k = 1 p 1 l > p 2 u k l 2 + k > p 1 l = 1 u k l 2 1 2 0 , β 2 < ( p 1 , p 2 ) .
From the Cauchy–Schwarz’s inequality and Yifan Sun (2020) [22] Lemmas 1 and 2, we obtain
I 1 = k = 1 p 1 l = 1 p 2 u k l ( φ 1 k φ 2 l φ ^ 1 k φ ^ 2 l ) 2 = T 1 T 2 k = 1 p 1 l = 1 p 2 u k l ( φ 1 k φ 2 l φ ^ 1 k φ ^ 2 l ) 2 d s d t T 1 T 2 k = 1 p 1 l = 1 p 2 u k l 2 k = 1 p 1 l = 1 p 2 ( φ 1 k φ 2 l φ ^ 1 k φ ^ 2 l ) 2 d s d t β 2 k = 1 p 1 l = 1 p 2 φ 1 k φ 2 l φ ^ 1 k φ ^ 2 l 2 = β 2 k = 1 p 1 l = 1 p 2 C n 1 μ X 1 , k 2 + 1 μ X 2 , l 2 .
Therefore, we have
I 1 = O p k = 1 p 1 l = 1 p 2 1 n 1 μ X 1 , k 2 + 1 μ X 2 , l 2 0 .
Thus, we conclude that
β β ^ = O p t = 1 p 1 p 2 ( u t u ^ t ) 2 .
Therefore, Lemma 1 is proven. □
Theorem 1.
If the above conditions and assumptions hold, then we have
n d G 2 ( β ^ 1 , β 1 ) p 1 2 p 1 n d G 2 ( β ^ 2 , β 2 ) p 2 2 p 2 n d 2 ( β , β ^ ) p 1 p 2 τ 2 p 1 p 2 τ n Θ 0 ( γ 0 γ ^ 0 ) n Θ 1 ( γ 1 γ ^ 1 ) n Θ q ( γ q γ ^ q ) N ( 0 , I ) ,
where Θ m = E g ( i ) 2 σ 2 ( η i ) z i m 2 , τ = E g ( i ) 2 σ 2 ( η i ) ρ i t 2 , and I is a unit matrix of ( q + 1 + p 1 + p 2 ) × ( q + 1 + p 1 + p 2 ) .
Proof. 
A Taylor expansion-based approach is used to prove the asymptotic normality of the estimates. The Hessian of the proposed likelihood is J Ω = Δ Ω U ( Ω ) and
A T A = i = 1 n g ( i ) 2 σ 2 ( η i ) ω i ω i T .
Thus, we have
J Ω = U ( Ω ) Ω = U ( Ω ) i i Ω = i = 1 n g 2 ( i ) ω i ω i σ 2 ( g ( i ) ) + i = 1 n g ( i ) σ 2 ( η i ) g 2 ( i ) σ 2 ( η i ) σ 4 ( η i ) ( Y i g ( i ) ) ω i ω i = A T A + R .
The remainder term R can be ignored and, using Taylor expansions, we obtain a Ω ˜ that lies between Ω and Ω ^ . We have
U ( Ω ) U ( Ω ^ ) Ω Ω ^ = J Ω ˜ ;
therefore,
n ( Ω Ω ^ ) = I + M + N 1 A T A n 1 U ( Ω ) n ,
where M = A T A n 1 J Ω ˜ J Ω n , N = A T A n 1 J Ω A T A n .
From Lemma 7.1 in Müller (2005) [11], it follows that
n ( Ω Ω ^ ) A T A n 1 U ( Ω ) n .
Asymptotic convergence of the lower proof A T A n 1 U ( Ω ) n gives
A T A n 1 U ( Ω ) n = A T A n 1 A T V 1 2 ( Y η ) n = A T A n 1 A T ε ¯ n ,
where ε ¯ = ε σ ( η ) and follows a standard normal distribution.
Thus, we have
n ( Ω Ω ^ ) A T A n 1 A T ε ¯ n .
Since β j ( t j ) , β ( t 1 , t 2 ) , and γ are of different data types, n ( Ω Ω ^ ) is divided into three terms, i.e.,
n ( b j b ^ j ) A j T A j n 1 A j T ε ¯ n , j = 1 , 2 ,
n ( u u ^ ) A 12 T A 12 n 1 A 12 T ε ¯ n ,
n ( γ γ ^ ) A 0 T A 0 n 1 A 0 T ε ¯ n ,
where A j T A j are symmetric matrices and A 0 T A 0 are diagonal matrices.
First, we prove
n ( b j b ^ j ) A j T A j n 1 A j T ε ¯ n , j = 1 , 2 .
Let
X n j = Λ ˜ j 1 2 A j T ε ¯ n , Z n j = A j T A j n 1 A j T ε ¯ n ,
Ψ n j = Γ ˜ j 1 2 A j T A j n 1 Γ ˜ j 1 2 .
Thus, we have
n d G 2 ( β , β ^ ) = Z n j T Γ ˜ j Z n j = X n j T Ψ n j 2 X n j = X n j T X n j + 2 X n j T ( Ψ n j I n j ) X n j + X n j T ( Ψ n j I n j ) ( Ψ n j I n j ) X n j ,
From Lemma 7.2 in Müller (2005) [11], it follows that
n d G 2 ( β ^ j , β j ) = X n j T X n j .
Then,
X n j T X n j = 1 n k 1 = 1 p j k 2 = 1 p j ζ j , k 1 k 2 1 2 i = 1 n g ( η i ) χ i j k 2 σ ( μ i ) ε ¯ i 2 = E + F ,
where
E = 1 n i = 1 n ε ¯ i 2 k 2 , k 2 = 1 p j g ( i ) 2 σ 2 ( η i ) χ i j k 2 χ i j k 2 k 1 = 1 p j ζ j , k 1 k 2 1 2 ζ j , k 1 k 2 1 2 ,
F = 1 n i 1 i 2 = 1 n ε ¯ i 1 ε ¯ i 2 g ( i 1 ) g ( i 2 ) σ ( η i 1 ) σ ( η i 2 ) k 2 , k 2 = 1 p j χ i 1 j k 2 χ i 2 j k 2 k 1 = 1 p j ζ j , k 1 k 2 1 2 ζ j , k 1 k 2 1 2 .
Since ε ¯ follows a standard normal distribution, we have
E X n j T X n j = p j , V a r X n j T X n j = 2 p j .
Therefore,
n d G 2 ( β ^ j , β j ) p j 2 p j N ( 0 , 1 ) j = 1 , 2 .
The following proves that
n ( u u ^ ) A 12 T A 12 n 1 A 12 T ε ¯ n .
For the coefficient function of the interaction term, we have the metric d 2 ( β , β ^ ) = β β ^ 2 , so, according to Lemma 1, we have
d 2 ( β , β ^ ) = t = 1 p 1 p 2 ( u t u ^ t ) 2 = ( u u ^ ) T ( u u ^ ) .
Let
Q n t = A 12 T A 12 n 1 A 12 T ε ¯ n ,
A n t = A 12 T ε ¯ n , F n t = A 12 T A 12 n 1 .
Thus, we have
n d 2 ( β , β ^ ) = Q n t T Q n t = A n t T F n t 2 A n t = A n t T A n t + 2 A n t T ( F n t I n t ) A n t + A n t T ( F n t I n t ) ( F n t I n t ) A n t ,
From Lemma 7.2 in Müller (2005) [11], it follows that
n d 2 ( β , β ^ ) = A n t T A n t .
Then,
A n t T A n t = 1 n t = 1 p 1 p 2 i = 1 n g ( i ) ρ i t σ ( η i ) ε ¯ i 2 = 1 n i = 1 n ε ¯ i 2 t = 1 p 1 p 2 g ( i ) 2 σ 2 ( η i ) ρ i t 2 + 1 n i 1 i 2 = 1 n ε ¯ i 1 ε ¯ i 2 g ( i 1 ) σ ( η i 1 ) g ( i 2 ) σ ( η i 2 ) t = 1 p 1 p 2 ρ i 1 t ρ i 2 t .
We have
E A n k T A n k = p 1 p 2 τ , V a r A n k T A n k = 2 p 1 p 2 τ .
Thus, there is
n d 2 ( β , β ^ ) p 1 p 2 τ 2 p 1 p 2 τ N ( 0 , 1 ) .
Next, prove that
n ( γ γ ^ ) A 0 T A 0 n 1 A 0 T ε ¯ n .
Let
Z 0 = A 0 T A 0 n 1 A 0 T ε ¯ n .
Then, its matrix form is
Z 0 = n i = 1 n g ( i ) 2 σ 2 ( η i ) z i m 2 1 i = 1 n g ( i ) z i m σ ( η i ) ε ¯ i .
Therefore, we have
n γ m γ ^ m n i = 1 n g ( i ) 2 σ 2 ( η i ) z i m 2 1 i = 1 n g ( i ) z i m σ ( η i ) ε ¯ i .
Since
E n γ m γ ^ m = 0 ,
V a r n γ m γ ^ m = E g ( i ) 2 σ 2 ( η i ) z i m 2 1 = Θ m 1 .
There is
n Θ m γ m γ ^ m N ( 0 , 1 ) .
Therefore, Theorem 1 is proven. □

4. Simulation

In this simulation, we consider the case that has two functional predictors, three scalar predictors, an interaction term between the two functional predictors, and a binary response. In order to include the case in which the functional predictors do not have the same domain, we define the functional predictors X i 1 ( t 1 ) , t 1 [ 0 , 1 ] and X i 2 ( t 2 ) , t 2 [ 1 , 1 ] , i = 1 , , n , where n can be any positive integer. In the latter sample size, n takes the values of 50, 100, and 500, and, for each n, we run 100 simulations. First, we define two standard orthogonal bases φ 1 k ( t 1 ) , t 1 [ 0 , 1 ] and φ 2 l ( t 2 ) , t 2 [ 1 , 1 ] , satisfying
φ 1 k ( t 1 ) = 2 cos ( 2 k π t 1 ) , k = 1 , , 4 ,
φ 2 l ( t 2 ) = 2 sin ( 2 l π t 2 ) , l = 1 , , 5 .
Under the Gaussian assumption, we define the two randomly generated functional principal component scores χ i 1 k , χ i 2 l that satisfy
χ i 1 k N ( 0 , λ 1 k ) , k = 1 , , 4 ,
χ i 2 l N ( 0 , λ 2 l ) , l = 1 , , 5 ,
where λ 11 = 8 , λ 12 = 6 , λ 13 = 4 , λ 14 = 2 , λ 21 = 4 , λ 22 = 2 , λ 23 = 1 , λ 24 = 1 2 , λ 25 = 1 4 . Notice that the first three functional principal components explain up to 90% of the variation in the two predictors. So, we have
X i 1 ( t 1 ) = k = 1 4 χ i 1 k φ 1 k ( t 1 ) ,
X i 2 ( t 2 ) = l = 1 5 χ i 2 l φ 2 l ( t 2 ) .
Fifty images of X 1 ( t 1 ) and X 2 ( t 2 ) are shown in Figure 1.
For scalar predictors, we assume Z 1 N ( 0 , 2 ) , Z 2 N ( 0 , 5 ) , and Z 3 N ( 0 , 6 ) .
We assume that the theoretical values of the regression coefficients are
γ = ( 4 , 6 , 8 ) T ,
β 1 ( t 1 ) = k = 1 4 b 1 k φ 1 k ( t 1 ) ,
β 2 ( t 2 ) = l = 1 5 b 2 l φ 2 l ( t 2 ) ,
where b 1 k = k 5 , b 2 l = l 2 25 .
For the interaction term, its principal component score is denoted by ρ i k l and satisfies
ρ i k l = χ i 1 k χ i 2 l ,
ψ k l ( t 1 , t 2 ) = φ 1 k ( t 1 ) φ 2 l ( t 2 ) ,
β ( t 1 , t 2 ) = k = 1 4 l = 1 5 u k l φ 1 k ( t 1 ) φ 2 l ( t 2 ) ,
where u k l = k 10 2 .
The corresponding response variable is generated by
p ( X i , Z i ) = g T 1 X i 1 t 1 β 1 t 1 d t 1 + T 2 X i 2 t 2 β 2 t 2 d t 2 + T 1 × T 2 X i 1 t 1 X i 2 t 2 β t 1 , t 2 d t 1 d t 2 + Z i T γ ,
where the link function g ( x ) = exp ( x ) 1 + exp ( x ) and Y ( X , Z ) B e r n o u l l i ( p ( X , Z ) ) is a sequence of pseudo-random numbers.
The principal component analysis was performed for n = 50 , 100 , 500 , and the running results showed that the principal component scores of X 1 with 90 % cumulative contribution were 3, 3, 3 for each sample size, respectively and the principal component scores of X 2 with 90 % cumulative contribution were 2, 2, 2.
Table 1 shows how the standardized prediction error (SPE) varies with different sample sizes, and the results show that the model’s predictions become more and more accurate as the sample size increases. Here, SPE is defined by i | Y ^ i Y i | / i | Y i | .
Figure 2 shows β ^ 1 ( t 1 ) , β ^ 2 ( t 2 ) and the corresponding 95 % confidence interval bands for different sample sizes, where the red curves are the theoretical values of β 1 ( t 1 ) and β 2 ( t 2 ) and the black curves are the corresponding estimates β ^ 1 ( t 1 ) and β ^ 2 ( t 2 ) . From Figure 2, it can be seen that, as the sample size increases, the estimated value becomes closer to the theoretical value. Figure 3  β ^ ( t 1 , t 2 ) shows the visualized 3D plot with β ^ ( t 1 , t 2 ) in the middle panel and the 95 % confidence intervals for β ^ ( t 1 , t 2 ) in the left and right panels.
Table 2 shows the estimated values of γ ^ and their corresponding standard deviations for different sample sizes. It can be seen that, as n increases, the standard deviation becomes smaller and the estimated value of γ becomes closer to the theoretical value, where the theoretical values of γ are 4, 6, and 8, respectively. Table 3 shows the standard deviation and root mean square error for β ^ 1 , β ^ 2 and β ^ ( t 1 , t 2 ) for different sample sizes. Here, we use the coefficients of the basis expansion of the regression coefficient function to calculate the root mean square error. For example, the root mean square error of β ^ 1 is k = 1 4 ( b ^ 1 k b 1 k ) 2 . The results show that, as n increases, both the standard deviation and the RMS error become smaller, indicating that, as sample size increases, the prediction becomes more accurate.

5. Application

To investigate the influence of the influence of air qualities, climate factors, medical and social indicators, and their interactions on cancer incidence using the proposed model, we collected data on average daily PM2.5 concentration, average daily humidity, per capita GDP, green coverage rate in built-up areas, the proportion of medical personnel (PMP), and the incidence of cancer in 49 cities in China from http://www.cnemc.cn/, http://www.stats.gov.cn/sj/ndsj/, and http://www.chinancpcn.org.cn/home.
There are two functional predictors (average daily PM2.5 concentration and average daily humidity from 1 January 2015 to 31 December 2020), three scalar predictors (per capita GDP, greenery coverage, and PMP in 2020), and the response is the cancer incidence in 2020. The ratio of the number of new cancer cases to the total number of people in China in 2020 is 0.3156 % . The data of the cancer incidence can only contain 0 and 1, indicating high or low cancer incidence rate. When the cancer incidence of a city was less than 0.3156 % , the city was considered to have a low cancer incidence rate, denoted by 0; otherwise, the cancer incidence was high, denoted by 1. Figure 4 shows average daily PM2.5 concentration and daily relative humidity in 21 cities selected from the 49 cities.
We chose g ( x ) = exp ( x ) 1 + exp ( x ) as the link function. The model was first subjected to principal component analysis and then the number of principal components was determined based on the cumulative contribution to obtain the number of functional principal components for PM2.5 concentrations and relative humidity, which were chosen as p P M 2.5 = 7 , p Humidity = 14 in order to explain 75% of the variation.
The prediction accuracy is shown by the Generalized Cross Validation (GCV) with a value of 0.0038.
The results of the regression coefficients for the scalar predictor variable γ ^ are shown in Table 4, where we can see that the per capita GDP is positively correlated with the incidence of cancer, i.e., the higher the GDP per capita, the higher the incidence of cancer in that city, which is consistent with the findings of Cao et al. [27]. The reason for this situation is that the promotion of cancer screening, early diagnosis, and treatment in the more economically developed regions has, to some extent, facilitated the detection of the disease. The greenery coverage is negatively correlated with the cancer incidence, i.e., the higher the greenery coverage, the lower the cancer incidence, which is also consistent with the findings of Wu et al. [26]. A high green coverage rate implies better air quality, which in turn reduces the risk of cancer. Additionally, a high green coverage rate may provide more outdoor recreational spaces, promoting physical activity and exercise, contributing to maintaining good physical health, and, thus, reducing the risk of cancer. The PMP is positively correlated with the incidence of cancer. As we all know, cancer incidence is age-related, and older people are more susceptible to cancer. The higher PMP, the better the medical conditions, the longer the average life expectancy of the people, and, therefore, the higher the cancer incidence.
The regression coefficient functions β ^ 1 ( t 1 ) and β ^ 2 ( t 2 ) for the functional predictors are shown in Figure 5. From Figure 5, we can see that the effect of PM2.5 concentration on cancer incidence is generally positively correlated, i.e., the higher the PM2.5 concentration, the higher the cancer incidence. This result is consistent with Qin et al. [25] from 2014. Regarding the effect of humidity on cancer incidence, there is a more significant positive correlation between humidity and cancer incidence, i.e., the higher the humidity, the higher the cancer incidence. In high-humidity environments, there may be a higher presence of mold and fungi, and the spores and harmful substances released by these microorganisms may have negative effects on human health, increasing the risk of cancer. In high-humidity environments, pollutants in the air are more likely to adhere to suspended particles, making them more easily inhalable by humans. These pollutants include PM2.5, organic compounds, and heavy metals, which are believed to be associated with the occurrence of cancer. High humidity increases the survival time of bacteria and viruses in the air, increasing the chances of people becoming infected with diseases. Certain viruses such as hepatitis B virus and human papillomavirus (HPV) are believed to be associated with the occurrence of cancer.
The interaction surface estimate β ^ ( t 1 , t 2 ) (middle) ± two times the estimated standard errors (left and right) are given in Figure 6. Figure 7 shows the contour map of β ^ ( t 1 , t 2 ) , from which it can be seen that β ^ ( t 1 , t 2 ) decreases and then increases with t 1 when t 2 [ 0 , 1100 ] and increases and then decreases with t 1 when t 2 [ 1100 , 2192 ] . In the conditions of higher humidity, PM2.5 particles may be more prone to settling, reducing the suspended harmful particles in the air, potentially lowering the incidence of cancer. Conversely, in lower-humidity conditions, PM2.5 may be more likely to remain suspended in the air, increasing the risk of respiratory system exposure, thereby raising the incidence of cancer. Additionally, the concentrations of PM2.5 and humidity may not fluctuate synchronously throughout the day. By introducing interaction terms, the model can capture the temporal complexities, making the estimation results more in line with real-world conditions.
To verify the necessity of considering the interaction term, i.e., to demonstrate the effectiveness of our proposed method, we compare mod1 proposed in this paper with mod2, which does not include the interaction term, i.e.,
m o d 1 : Y i = g α + T 1 X i 1 t 1 β 1 t 1 d t 1 + T 2 X i 2 t 2 β 2 t 2 d t 2 + T 1 × T 2 X i 1 t 1 X i 2 t 2 β t 1 , t 2 d t 1 d t 2 + Z i T γ + ε i .
m o d 2 : Y i = g α + T 1 X i 1 t 1 β 1 t 1 d t 1 + T 2 X i 2 t 2 β 2 t 2 d t 2 + Z i T γ + ε i .
The general standards for evaluating model performance are AIC (Akaike Information Criterion), residual, R-squared, RMSE (root mean square error), and MAE (mean absolute error). The smaller values of AIC, residuals, RMSE, and MAE indicate that the model’s fitting effect and generalization ability are better. The R-squared takes a value between 0 and 1, and, the bigger the value, the better the model’s fitting effect. According to Table 5, we can see that the AIC, residuals, RMSE, and MAE values of mod1 are smaller and that R-squared is much closer to 1 compared to that of mod2, which indicates that mod1 has a better performance. Thus, including the interaction term between PM2.5 concentration and relative humidity makes the research results more meaningful.

6. Discussion

This paper proposes a generalized partially functional linear model with interaction terms. We first use principal component analysis to reduce the dimensionality of the functional data, followed by maximum likelihood estimation to obtain the estimates of the unknown parameters, then prove the asymptotic property of the estimators, and finally perform data simulations and apply our model to a real data example.
As the incidence and mortality of cancer in China are increasing year by year, it is necessary to study the influencing factors and formulate corresponding measures. The effect of PM2.5 concentration, average daily humidity, per capita GDP, the greenery coverage of built-up areas, and PMP on cancer incidence in 49 cities in China was investigated, which showed that the effect of PM2.5 concentration and relative humidity on cancer incidence was generally positively correlated. The effect of greenery coverage in built-up areas on cancer incidence is negatively correlated, while the effect of per capita GDP and the proportion of medical personnel on cancer incidence is positively correlated. The higher the economic level and the more developed the medical conditions, the longer the average life expectancy of people and, therefore, the higher the cancer incidence. Comparing this model with the model without the interaction term shows that considering the role of the interaction term leads to more accurate and meaningful predictions.
Our research lays a foundation for further study on the generalized partially functional linear model with interaction terms and of unknown link function or variance function.

Author Contributions

W.X.: methodology, validation, writing—review, supervision, funding acquisition. K.M.: methodology, software, data curation, writing—original draft. H.L.: writing—review, supervision. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by the Yujie Talent Project of North China University of Technology (Grant No. 107051360024XN153).

Data Availability Statement

The original data supporting the results of this study can be obtained from the National Meteorological Science Data Sharing Service Platform, the National Environmental Monitoring Station, and local statistical bulletins.

Acknowledgments

The authors would like to thank the referees and the editor for their useful suggestions, which helped us to improve the paper.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Ramsay, J.O.; Silverman, B.W. Functional Data Analysis, 2nd ed.; Springer: New York, NY, USA, 2005. [Google Scholar]
  2. Horváth, L.; Kokoszka, P. Inference for Functional Data with Application; Springer: New York, NY, USA, 2012. [Google Scholar]
  3. Hsing, T.; Eubank, R. Theoretical Foundations of Functional Data Analysis, with an Introduction to Linear Operators; Springer: New York, NY, USA, 2015. [Google Scholar]
  4. Cardot, H.; Ferraty, F.; Sarda, P. Functional linear model. Stat. Probab. Lett. 1999, 45, 11–22. [Google Scholar] [CrossRef]
  5. Tony, C.; Peter, H. Prediction in functional linear regression. Ann. Stat. 2006, 34, 2159–2179. [Google Scholar] [CrossRef]
  6. Cardot, H.; Crambes, C.; Kneip, A.; Sarda, P. Smoothing Splines Estimators in Functional Linear Regression with Errors in Variables. Comput. Stat. Data Anal. 2007, 51, 4832–4848. [Google Scholar] [CrossRef]
  7. Delaigle, A.; Hall, P. Methodology and theory for partial least squares applied to functional data. Ann. Stat. 2012, 40, 322–352. [Google Scholar] [CrossRef]
  8. Cai, T.; Yuan, M. Minimax and adaptive prediction for functional linear regression. J. Am. Stat. Assoc. 2012, 107, 1201–1216. [Google Scholar] [CrossRef]
  9. Nelder, J.A.; Wedderburn, R.W.M. Generalized Linear Models. J. R. Stat. Soc. 1972, 135, 370–384. [Google Scholar] [CrossRef]
  10. James, G.M. Generalized linear models with functional predictors. J. R. Stat. Soc. Ser. B 2002, 64, 411–432. [Google Scholar] [CrossRef]
  11. Müller, H.G.; Stadtmüller, U. Generalized functional linear models. Ann. Stat. 2005, 33, 774–805. [Google Scholar] [CrossRef]
  12. Goldsmith, J.; Bobb, J.; Crainiceanu, C.M.; Caffo, B.; Reich, D. Penalized Functional Regression. J. Comput. Graph. Stat. 2011, 20, 830–851. [Google Scholar] [CrossRef]
  13. Xiao, W.W.; Wang, Y.X.; Liu, H.Y. Generalized partially functional linear model. Sci. Rep. 2021, 11, 23428. [Google Scholar] [CrossRef]
  14. Kong, D.; Xue, K.; Yao, F.; Zhang, H.H. Partially functional linear regression in high dimensions. Biometrika 2016, 103, 147–159. [Google Scholar] [CrossRef]
  15. Yao, F.; Sue, C.S.; Wang, F. Regularized partially functional quantile regression. J. Multivar. Anal. 2017, 156, 39–56. [Google Scholar] [CrossRef]
  16. Ma, H.; Li, T.; Zhu, H.; Zhu, Z. Quantile regression for functional partially linear model in ultra-high dimensions. Comput. Stat. Data Anal. 2019, 129, 135–147. [Google Scholar] [CrossRef]
  17. Xu, W.; Ding, H.; Zhang, R.; Liang, H. Estimation and inference in partially functional linear regression with multiple functional covariates. J. Stat. Plan. Inference 2020, 209, 44–61. [Google Scholar] [CrossRef]
  18. Usset, J.; Staicu, A.M.; Maity, A. Interaction models for functional regression. Comput. Stat. Data Anal. 2016, 94, 317–329. [Google Scholar] [CrossRef]
  19. Luo, R.; Qi, X. Interaction Model and Model Selection for Function-on-Function Regression. J. Comput. Graph. Stat. 2019, 28, 309–322. [Google Scholar] [CrossRef]
  20. Yang, W.H.; Wikle, C.K.; Holan, S.H.; Wildhaber, M.L. Ecological Prediction with Nonlinear Multivariate Time-Frequency Functional Data Models. J. Agric. Biol. Environ. Stat. 2013, 18, 450–474. [Google Scholar] [CrossRef]
  21. Matsui, H. Quadratic regression for functional response models. Econom. Stat. 2020, 13, 125–136. [Google Scholar] [CrossRef]
  22. Sun, Y.; Wang, Q. Function-on-function quadratic regression models. Comput. Stat. Data Anal. 2020, 142, 106814. [Google Scholar] [CrossRef]
  23. Fuchs, K.; Scheipl, F.; Greven, S. Penalized scalar-on-functions regression with interaction term. Comput. Stat. Data Anal. 2015, 81, 38–51. [Google Scholar] [CrossRef]
  24. Qiu, H.; Cao, S.; Xu, R. Cancer incidence, mortality, and burden in China: A time-trend analysis and comparison with the United States and United Kingdom based on the global epidemiological data released in 2020. Cancer Commun. 2021, 41, 1037–1048. [Google Scholar] [CrossRef]
  25. Qin, X.; Wan, F.; Zhang, H.; Dai, B.; Shi, G.; Zhu, Y.; Ye, D. Relationship between air pollution PM2.5 concentration and cancer. In Proceedings of the 8th Chinese Oncology Academic Conference and the 13th Cross-Strait Oncology Academic Conference, Jinan, China, 12–13 September 2014; p. 76. [Google Scholar]
  26. Wu, X.; Feng, Y.; Chang, K.; Jia, X.; Xue, F. Analysis of the causal relationship between green coverage and the incidence of cancer. J. Shandong Univ. (Health Sci.) 2022, 60, 115–119. [Google Scholar]
  27. Cao, W.; Li, F.; Liang, Y.; Yu, D. Analysis of the relationship between the level of economic development and cancer incidence and mortality in selected regions of China. Chin. J. Dis. Control Prev. 2023, 27, 209–215. [Google Scholar] [CrossRef]
  28. Xu, J.; Kuang, H.Y.; Wang, G.Q.; Chen, M.; Lu, L. Analysis of the relationship between PM2.5 and air relative humidity. Agric. Technol. 2017, 37, 148–149. [Google Scholar]
  29. Yang, Z.; Wang, Y.K.; Xu, X.H.; Yang, J.; Ou, C.Q. Quantifying and characterizing the impacts of PM2.5 and humidity on atmospheric visibility in 182 Chinese cities: A nationwide time-series study. J. Clean. Prod. 2022, 368, 133182. [Google Scholar] [CrossRef]
Figure 1. Functional predictors X 1 ( t 1 ) and X 2 ( t 2 ) .
Figure 1. Functional predictors X 1 ( t 1 ) and X 2 ( t 2 ) .
Axioms 13 00583 g001
Figure 2. Estimated regression coefficient functions β ^ 1 ( t 1 ) , β ^ 2 ( t 2 ) (black curves) and their 95 % confidence bands (grey area) for different sample sizes, where the red curves are the theoretical regression coefficient functions β 1 ( t 1 ) , β 2 ( t 2 ) .
Figure 2. Estimated regression coefficient functions β ^ 1 ( t 1 ) , β ^ 2 ( t 2 ) (black curves) and their 95 % confidence bands (grey area) for different sample sizes, where the red curves are the theoretical regression coefficient functions β 1 ( t 1 ) , β 2 ( t 2 ) .
Axioms 13 00583 g002aAxioms 13 00583 g002b
Figure 3. β ( t 1 , t 2 ) and β ^ ( t 1 , t 2 ) are visualized in 3D.
Figure 3. β ( t 1 , t 2 ) and β ^ ( t 1 , t 2 ) are visualized in 3D.
Axioms 13 00583 g003
Figure 4. PM2.5 concentrations and average daily humidity in 21 selected cities in 2020.
Figure 4. PM2.5 concentrations and average daily humidity in 21 selected cities in 2020.
Axioms 13 00583 g004
Figure 5. Regression coefficient functions β ^ 1 ( t 1 ) , β ^ 2 ( t 2 ) and their 95 % confidence bands.
Figure 5. Regression coefficient functions β ^ 1 ( t 1 ) , β ^ 2 ( t 2 ) and their 95 % confidence bands.
Axioms 13 00583 g005
Figure 6. Visualization of β ^ ( t 1 , t 2 ) in 3D.
Figure 6. Visualization of β ^ ( t 1 , t 2 ) in 3D.
Axioms 13 00583 g006
Figure 7. Contour map of β ^ ( t 1 , t 2 ) .
Figure 7. Contour map of β ^ ( t 1 , t 2 ) .
Axioms 13 00583 g007
Table 1. Standardized prediction error for different sample sizes.
Table 1. Standardized prediction error for different sample sizes.
nSPE
500.0156
1000.0130
5000.0106
Table 2. Estimates of the regression coefficients and their standard deviations.
Table 2. Estimates of the regression coefficients and their standard deviations.
n γ ^ 1 γ ^ 2 γ ^ 3
503.912 (0.193)5.902 (0.068)8.053 (0.045)
1004.013 (0.062)5.958 (0.026)8.015 (0.018)
5003.998 (0.022)5.993 (0.006)7.996 (0.007)
Table 3. Standard deviation and root mean square error of the estimated values of the regression coefficient function.
Table 3. Standard deviation and root mean square error of the estimated values of the regression coefficient function.
nSdRMSE
β ^ 1 500.0340.026
1000.0180.006
5000.0080.001
β ^ 2 500.3250.204
1000.1260.044
5000.0430.004
β ^ ( t 1 , t 2 ) 500.1370.094
1000.0540.024
5000.0200.004
Table 4. Estimates of regression coefficients and their levels of significance.
Table 4. Estimates of regression coefficients and their levels of significance.
EstimateStd. Errort ValuePr (>|t|)
GDP1.165 × 10−62.632 × 10−74.4246.45 × 10−5 ***
Greenery coverage−2.062 × 10−11.920 × 10−1−2.6540.0378 *
PMP3.676 × 10−41.642 × 10−42.2390.0491 *
PMP—Proportion of medical personnel. * indicates that the result is significant at the 5% significance level. *** indicates that the result is significant at the 0.1% significance level.
Table 5. Results of model comparison.
Table 5. Results of model comparison.
AICR-SquaredResidualRMSEMAE
mod18.2810.92870.78160.12630.1036
mod235.5920.64653.21580.25620.1989
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

Xiao, W.; Mao, K.; Liu, H. Generalized Partially Functional Linear Model with Interaction between Functional Predictors. Axioms 2024, 13, 583. https://doi.org/10.3390/axioms13090583

AMA Style

Xiao W, Mao K, Liu H. Generalized Partially Functional Linear Model with Interaction between Functional Predictors. Axioms. 2024; 13(9):583. https://doi.org/10.3390/axioms13090583

Chicago/Turabian Style

Xiao, Weiwei, Kejing Mao, and Haiyan Liu. 2024. "Generalized Partially Functional Linear Model with Interaction between Functional Predictors" Axioms 13, no. 9: 583. https://doi.org/10.3390/axioms13090583

APA Style

Xiao, W., Mao, K., & Liu, H. (2024). Generalized Partially Functional Linear Model with Interaction between Functional Predictors. Axioms, 13(9), 583. https://doi.org/10.3390/axioms13090583

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