Functional-coefficient Spatial Durbin Models with Nonparametric Spatial Weights: an Application to Economic Growth

This paper considers a functional-coefficient spatial Durbin model with nonparametric spatial weights. Applying the series approximation method, we estimate the unknown functional coefficients and spatial weighting functions via a nonparametric two-stage least squares (or 2SLS) estimation method. To further improve estimation accuracy, we also construct a second-step estimator of the unknown functional coefficients by a local linear regression approach. Some Monte Carlo simulation results are reported to assess the finite sample performance of our proposed estimators. We then apply the proposed model to reexamine national economic growth by augmenting the conventional Solow economic growth convergence model with unknown spatial interactive structures of the national economy, as well as country-specific Solow parameters, where the spatial weighting functions and Solow parameters are allowed to be a function of geographical distance and the countries' openness to trade, respectively.


Introduction
Ever since the seminal work of [1], there has been a significant amount of empirical work studying variation in economic growth rates across countries.Particularly, more and more economists have paid attention to the impacts of economic interaction and spillover effects on the regional and national economy in the past two decades; see, e.g., [2] for taxation and the global allocation of capital, [3] for cross-border foreign direct investment decisions, [4,5] for economic growth models with worldwide interactions, [6] for country interactions in discretionary fiscal policy and [7] for an overview of empirical studies of strategic interaction among governments over environmental standards and public expenditures.In the meanwhile, econometric theory in parametric spatial regression models has been introduced and well developed to analyse the spatial and economic externalities; for a detailed survey on parametric spatial econometric models, see [8][9][10][11].This paper joins the others to examine the impact of cross-country economic externalities on national growth through a Solow growth model augmented with economic externalities.
The role of spatial dependence in regional economic growth has received substantial attention in the empirical growth literature in the recent decade; see, e.g., [12] for a survey on economic growth and space.It has been recognized that a nation's per capita GDP growth rate is affected not only by its own values of determinants, such as savings, population growth rate and initial level of income, but also by its neighbouring nations' per capita GDP growth rates and the values of these determinants.For example, Ertur and Koch [4] developed a theoretical growth model with spatial externality resulting from technological interdependence among economies and proposed a spatially-augmented Solow economic growth model yielding a conditional convergence equation with heterogeneous Solow parameters.Note that heterogeneous Solow parameters are also supported by similar studies with no spatial interactions; see, e.g., [13][14][15].
Fitting a spatial Durbin model (SDM) using data from 91 non-oil regions/and countries for the period from 1960 to 1995, Ertur and Koch [4] found positive and significant spatial dependence across these economies together with predicted signs for all coefficients.However, their study suffers from two potential problems.First, the parametric SDM requires researchers to pre-determine the non-stochastic spatial dependence structure among economies before estimating parameters appearing in the model, and the misspecified spatial interactive relations can incur inconsistent estimation and misleading inference.Second, the subsampling method may not be the best way of studying heterogeneous Solow parameters.In this paper, working on the same dataset used in [4], we therefore aim to re-examine the spatial spillover effects of economic growth, while estimating in a nonparametric way the true spatial dependence structure among economies and allowing the Solow parameters to vary with respect to the trade openness of an economy.
Specifically, we propose a functional-coefficient spatial Durbin model with nonparametric spatial weights and estimate the unknown spatial weights and coefficient curves via a series approximation approach by a nonparametric two-stage least squares method.Based on the first-step consistent estimator, we then construct a second-step estimator for the unknown functional coefficients, which is oracle efficient in the sense that the limit distribution of the second-step estimator is the same regardless of whether the spatial weights are known.Moreover, we give our inference on spatial dependence through average direct and indirect impact values with standard errors calculated from the bootstrap method.
The remainder of this paper is organized as follows.Section 2 introduces our proposed semiparametric spatial Durbin model.Section 3 presents our estimation methodology.Section 4 reports results from a small Monte Carlo simulation to examine the finite sample performance of our proposed estimators.Section 5 gives our empirical results.Section 6 concludes.

Model
In the empirical economic growth literature, DeLong and Summers [16], to the best of our knowledge, is the first study to investigate spatial correlation taking geographical distance into account.Using a sample of 61 countries for the period from 1960 to 1985, they find no significant spatial correlation in their sample.Moreno and Trehan [17], on the other hand, augments [1]'s model with a spatial interactive term and find highly significant spillover effects between geographical neighbours, and they argue that using a border dummy variable instead of a spatial lag term neglects the influence of neighbour countries that do not have a common border with the country of interest; relevant literature includes [18][19][20].Moreover, Ertur, le Gallo, and Baumont [21] provide strong evidence of spatial dependence in economic convergence processes among European regional economies.Using 155 European regions over the period 1988-2000, Basile [22] also finds some evidence of spatial spillovers across countries.
The question of how to measure spatial interactive relations between any pair of spatial units is answered by defining a neighbourhood set for each spatial unit according to some selected relevant variables.For example, Cliff and Ord [8] specify spatial weights for spatial unit i as the ratio of the length of the common border between units i and j to the geographical distance between them, , with some parameters α > 0 and β > 0. A common approach in practice is to use only distance-based weights with a decay parameter α; i.e., the spatial weight from unit j on unit i is defined as g ij = g d ij , α , where g (•) is a known function and α is a parameter to be estimated.The popularly-used distance function includes an inverse power function g ij = 1/d α ij or a negative exponential function g ij = exp(−αd ij ) (e.g., [23]) for some α > 0.Moreover, if a cut-off distance is not used, then a non-sparse spatial weight matrix is constructed.This implies that every region is a neighbour of other regions, but the spatial weights depreciate as the distance between two regions increases.For recent studies on spatial weights, see [24,25].Moreover, the term "neighbours" can also refer to contiguities defined by economic distances (e.g., [26][27][28]) and social networks ( [29]).
As the functional form g (•) is unknown in practice, to avoid misspecifying the wrong spatial weighting function, one can alternatively estimate the unknown spatial weighting function g (•) from the data via a nonparametric series estimation method.Compared to the parametric spatial modelling approach, the nonparametric approach enables researchers to impose less restrictive assumptions on the spatial weight function; see [27,30] for details.Alternatively, Ahrens and Bhattacharjee [31] proposed to estimate the unknown spatial weights via the LASSO estimation method when the unknown spatial weights matrix is sufficiently sparse.
Ertur and Koch [4] derive a theoretical Solow economic growth model augmented with global technological interdependence.They then approximate their theoretical model by a parametric spatial Durbin model via a linearization procedure and calculate the spatial weights from both the inverse power function and the exponential function of geographic distances with α = 2 for the sake of robustness, as the true spatial weighting function is unknown.As the linearization and selected parametric spatial weights may result in a model misspecification problem, in this paper, in order to better approximate [4]'s theoretical model, we therefore propose a semiparametric growth model that extends [4]'s parametric model by allowing nonparametric spatial weights, as well as varying Solow coefficients.Specifically, our proposed functional-coefficient spatial Durbin model with nonparametric spatial weights is given by: where Y i is a scalar dependent variable, X i = X 1i , . . ., X pi T is a p × 1 vector, D i is a continuous scalar random variable and Z is a non-stochastic spatial covariate with Z ii = 0 and Z ij > 0 for i = j.Moreover, g(•), m t (•) for t = 1, ..., p and θ(•) are all unknown measurable smooth functions with g(0) = 0 and m t (0) = 0, t = 1, .., p.In Model (1), X i has a functional coefficient depending on D i , and the unknown spatial weights are a function of non-stochastic geographic distance Z ij .The first two terms in the right-hand side of Equation (1), ∑ j =i g(Z ij )Y j and ∑ j =i m t (Z ij )X tj for each t = 1, ..., p, are called the spatial lag of the dependent variable and the spatially-lagged exogenous variables, respectively.In the spatial econometrics literature, this kind of model specification is referred to as the spatial Durbin model, where (X i , D i , u i ) are independently distributed across i while only the dependent variable Y i is dependently distributed across spatial units.The detailed data information on X i , D i , Y i and Z ij is delayed to Section 5.If the model includes only the spatial lag of the dependent variable, it is called a pure spatial autoregressive (SAR) model.Basile et al. [32] introduced the spatial autoregressive semiparametric geoadditive models to account for spatial dependence, spatial unobserved heterogeneity and unknown functional curves of regressors simultaneously, where the spatial autoregression is represented by a pre-determined spatial lag term of the dependent variable.Let G n and M n,t be n × n unknown spatial weight matrices with its (i, j) th element being g ij = g(Z ij ) and m t,ij = m t (Z ij ), respectively, for t = 1, ..., p, i = 1, ..., n, j = 1, ..., n.We then obtain a reduced form of Model (1) written in matrix form: where Y = [Y 1 , . . . ,Y n ] T , X t = [X t1 , • • • , X tn ] T and U = [u 1 , . . . ,u n ] T are all n × 1 vectors and mtk {X, θ(D)} is an n × 1 vector with the i th element equal to X T i θ(D i ).Furthermore, I n is an n × n identity matrix. Let j=1 a ij and A n 1 = max 1≤j≤n ∑ n i=1 a ij be the respective row and column norm of A n .Furthermore, C > 0 is a finite positive number that takes different values at different appearances.Below, we impose some regularity conditions on Model (1).
Assumption A1: (i) {Y i } is generated from Model (1), and {(X i , D i )} is independently distributed with finite second moments; (ii) g (•), m t (•) and θ (•) are all uniformly bounded up to their respective p th -order derivatives for some p > 2; (iii) {u i } is an independent sequence with zero mean, Assumption A2: (i) There exist a positive integer N and a constant c G ∈ (0, 1), such that for all Assumption A1 (i) states that the explanatory variables (X i , D i ) are independent, while the dependent variable Y i exhibits spatial dependence; and Assumption A1 (iii) allows the error term, u i , to be independent with heteroskedasticity, and the bounded higher order moment is required for deriving the limiting normal distribution of the proposed estimator.By [33] n , which implies that {Y i } is spatially stationary.In addition, we have max 1≤i,j≤n g ij ≤ ρ (G n ) < 1 by Properties 4.66 and 4.67 in [33] is continuously differentiable and uniformly bounded under Assumptions A1 and A2.In addition, Assumption A2 (ii) is a regularity condition (see, e.g., Assumption 1 in [34]), and it holds if the spatial weight function, g (z), decreases to zero for large z and: where the indicator function I (A) = 1 if event A holds, and zero otherwise.Note that a parametric SDM is given by n, where W n and P n,t are spatial weight matrices with their respective (i, j) th element equal to w ij and p t,ij .Therefore, if the parametric spatial Durbin model holds true, the spatial weight matrices G n and M n,t in Model (1) are equivalent to ρW n and δP n,t in the spatial Durbin model, respectively.From an estimation and econometric modelling viewpoint, the normalization of spatial weight matrices in the spatial Durbin model is used to identify the spatial multiplier parameters (ρ, δ), but this is not necessary in our proposed Model (1).Therefore, allowing nonparametric spatial weights saves us from applying an ad hoc spatial weight matrix normalization procedure as in the parametric SDM.

Estimation Methodology
If the dependent variable exhibits spatial autocorrelation, it must be accounted for by incorporating the spatially-lagged dependent variable into the model.If this variable is not included in the model, there would be an omitted variable type specification error due to the fact that unobserved factors may have a direct effect on the response variable.Moreover, the presence of the spatial lag-dependent variable in the model results in a simultaneity bias problem.This can be seen explicitly from the reduced form Model (2).From Model (2), we see a non-zero matrix, so that the spatial lag of the dependent variable, G n Y, is correlated with the error term, U, which, therefore, results in an endogeneity problem in Model (1).Therefore, the ordinary least squares (or OLS) estimator would be biased and inconsistent.Moreover, the term (I n − G n ) −1 in Equation ( 2) explains that region i is affected not only by its own determinants, but also by its neighbouring regions' values.This has been called a global interaction effect in [4] (p. 1044) and [24] (p.15).Another source of the spatial endogeneity problem is due to the endogeneity of spatial covariate in the model.Recently, Kelejian and Piras [35] and Sun [30] estimated the spatial panel data model and the SAR model with an endogenous spatial weight matrix in a nonparametric way, respectively.Moreover, Qu and Lee [36] proposed estimators for the parametric SAR model with an endogenous spatial covariate.This paper only deals with the endogeneity of the spatial lag of the dependent variable, as our spatial weight matrix is non-stochastic.
The endogeneity problem can be addressed by using the maximum likelihood estimation (or MLE) method, as well as the instrumental variable (or IV) approach.Ord [37] was the first to examine the MLE of SAR models.He proposed to use the eigenvalues of the spatial weights matrix to alleviate the computational complexity of the MLE method in large sample sizes.Lee [38] derived the large sample properties of the quasi-MLE without a normality assumption on error terms, while Bao and Ullah [39] obtained the second order bias of the maximum likelihood estimator for spatial autoregressive models.As the (quasi-) maximum likelihood estimator can be computationally difficult in moderate or large-sized samples, Kelejian and Prucha [40] proposed a two-stage least squares (or 2SLS) estimator for a SAR model with spatial autoregressive errors, while Lee [41] proposed an asymptotically-optimal 2SLS estimator.As the spatial weights in Model (1) are unknown, the 2SLS estimation methods derived in [40,41] are not feasible; we therefore use a series approximation method to recover the unknown spatial weight function and estimate all unknown functions via a nonparametric 2SLS (or NP2SLS) estimation method.For an overview of the sieve estimation method, see [42].
Assumption A3: (i) There exists a positive constant sequence {v n }, such that: (ii) There exist L n × 1 vectors α, γ t , and β t for t ∈ {1, ..., p}, such that: and: sup It is not necessary to know the exact order of υ n in Assumption A3 (i), and the consistency of our proposed estimator does not require υ n ≡ 1, as assumed in [27]'s Assumption (vi).From approximation theory in mathematics, Assumption A1 (ii) is a necessary condition for Assumption A3 (ii).However, ( 8) and ( 9) also require spatial units expanding sparsely as more spatial units are included, for example when (3) holds true; and the consistency of our proposed estimator relies on increasing domain asymptotic theory.Moreover, we use Laguerre polynomial series to approximate the unknown functions, as it is one of the common choices for series expansions when a function has a domain over [0, ∞) ( [42] (p.5574)).In addition, L n acts as a smoothing parameter that increases slowly with the sample size.In other words, it is required to have L n → ∞ and L n /n → 0 as n → ∞.An introduction of series estimation methods in a nonparametric framework can be found in [43] (Chapter 15).Now, we approximate Model (1) by: To derive our first-step estimator, we rewrite Model (10) in matrix form as follows: where we denote the i th row vector of an n × [(2p + 1)L n ] matrix V n by: T .The specification of the instrumental variable matrix is of great importance to obtain a consistent estimator.Since the number of endogenous variables increases with the number of approximating functions, L n , it is intuitively appealing to instrument the endogenous variables, .., p} for p > 1 as in our empirical application; see, e.g., [44].Since X t 1 j X t 2 i and D j are exogenous and relevant in predicting Y j , we would expect the proposed instrumental variables to serve as valid instruments for ∑ j =i Φ L n (Z ij )Y j .Therefore, we define the i th row vector of an n × [(2p + 2)L n ] instrumental matrix Q n as: We then can estimate ξ from (11) by the 2SLS estimation method.Note that we do not pursue optimal instrument variables in this paper due to the complexity of this approach in the semiparametric setup and the fact that the oracle efficiency of the second-step estimator of θ (•) does not rely on the use of optimal instruments in the first-step estimation.
To ensure the existence of our 2SLS estimator, we assume that the exogenous regressors matrix X n , the instrumental variables matrix Q n and V V n all have full column rank.Moreover, for the relevance of the instruments, we assume that E[Q T n V n ] has a full column rank.Otherwise, we can remove linearly-dependent terms as long as the number of instruments in Q n is more than the number of endogenous variables L n plus the number of exogenous regressors 2p. Lee [45] (p.493) argues that the 2SLS estimator would be inconsistent if (X i , D i ) are both irrelevant in predicting {Y i }.Therefore, throughout this paper, we assume that X and D contain relevant variables in predicting {Y i } and θ (•) takes non-zero values over any non-empty interval, so that there is no need to use the quadratic moments as additional orthogonal relations, as suggested in [45].Our empirical application in this paper satisfies this assumption by both economic theory and empirical findings observed from the economic growth literature.
To construct a consistent estimator for g(•), m t (•) and θ t (•), t = 1, ..., p, we consider the following nonparametric 2SLS objective function: The nonparametric 2SLS estimator of ξ solves (12) and is given by: and hence, the corresponding nonparametric 2SLS estimators 1 of unknown functions are given by: Next, we propose a second-step estimator for the functional coefficients, θ(d), using the local linear regression approach.We would expect the local linear estimate of θ(d), θ(d), to have an improvement over the first-step estimator, θ(d).Sun [30] considered a semiparametric spatial autoregressive model that has a mathematical representation of Model (1) with M n,t = 0 for all t = 1, . . ., p and has recently shown that the local linear estimator of θ(•) can be oracle efficient under some regularity conditions in the sense that its limiting distribution does not depend on whether or not the spatial weights are known.This is a general result from non-/semi-parametric additive models.As the unknown functions g Z ij , m 1 Z ij , . . ., m p Z ij and θ (D i ) enter Model (1) additively, we expect that θ(d) is oracle efficient, as well.We assume that as n → ∞, h → 0, nh → ∞ and nh 5 → c ∈ (0, ∞), where h is the bandwidth, which controls the size of the local neighbourhood around an interior point d.Moreover, let K(•) be a kernel function, which assigns more weights to the data closer to point d, satisfying: (i) K(a)da = 1; (ii) K(a) = K(−a); and (iii) a 2 K(a)da > 0.
The estimation procedure for θ(d) is given as follows: (i) We replace g(z) and m t (z) in (1) by ĝ(z) and mt (z), respectively, and treat Ŷi The asymptotic properties of the series estimator without an endogenous variable follows from [46] and with some extensions from [47].Moreover, we expect the asymptotic properties of our proposed estimators closely follow [30].
(ii) Applying the first-order Taylor series expansion of θ(D) around d, θ(D) ≈ θ(d) + θ (d)(D − d), we calculate the local linear estimator from a minimization of a kernel-weighted objective function: where θ(d) estimates θ(d) and θ (d) estimates θ (d), the first order derivative of θ (d).
For a complete treatment of local linear estimator, see [48].As the mathematical proofs of the consistency of the first-step estimator of ĝ(z), m1 (z), . . ., mp (z), θ(d) T T and the limiting result of the second-step estimator θ(d) closely follow those given in [30], the proofs are omitted from the paper.

Monte Carlo Simulations
In this section, we present the results from a very small Monte Carlo simulation study to assess the finite-sample properties of our estimators and more simulation results can be obtained from the authors upon request.We generate the data from the following regression model: where we randomly draw u i ∼ i.i.d.N(0, 0.5), D i ∼ i.i.d.U[0, 1] and X i = 0.5D i +η i with η i ∼ i.i.d.N(0, 1) independent of {u i }.For the exogenous variable, Z, we first randomly generate n observations from the U[0, R n ] distribution with R n = 0.001n 1.6 , by which we control the sparseness of spatial units.Then, we calculate Z ij as the absolute distance between observations i and j.
The specification of spatial weight functions requires that g(•) and m(•) are both decreasing and non-negative functions.We therefore set g(z) = m(z) = 0.01exp(−z/0.01)for z > 0 with g(0) = 0 and m(0) = 0.The random variables, u i , D i and Z i , are all mutually independent.We consider a sample size n ∈ {100, 200, 400}.The number of replications is 1000 for each n in the Monte Carlo experiments.Moreover, we set L n = 1, 2, 3 for each sample size, respectively.In the second-step estimation of coefficient functions, we select the bandwidth via a cross-validation method and use the Gaussian kernel function.To measure the performance of the estimators, we compute the root mean squared errors (or RMSEs) for each simulation.In Table 1, we report the averages of the RMSEs computed over 1000 repetitions, where ĝ(•), m(•) and θ(•) denote the NP2SLS estimators of g(•), m(•) and θ(•), respectively; θ(•) is the second-step estimator of θ(•), and θ(•) is the local linear estimator of θ(•), while g(•) and m(•) are known.Furthermore, we also estimate the average direct impact (ADI) and the average indirect impact (AII) and report their corresponding RMSEs in Table 1.Specifically, we first obtain the reduced form model from (17): where Y, X, θ (D) and U are all n × 1 vectors, and " " denotes the Hadamard multiplication.Then, the expected marginal effect of X is given by the following n × n matrix: from which we obtain ADI= n −1 tr {S (G n , M n , D)} and AII= n −1 i n S (G n , M n , D) i n −ADI; see LeSage and Pace [10], where i n is the n × 1 vector of ones and diag{θ (D)} is a diagonal matrix.
Replacing the two unknown spatial weight matrices and θ (D) by their estimates, we obtain the estimates for ADI and AII.
From Table 1, we observe that there is a decrease in the RMSEs for all three estimators as the sample size increases in each design.Moreover, the second-step estimator always performs better than the nonparametric 2SLS estimator.The relative ratios of the RMSEs of the second-step estimator and θ(•) generally reduce as the sample size increases.Therefore, our simulation results support the consistency of our proposed estimators.As for the ADI and AII, we also see an overall decreasing pattern in the RMSEs as the sample size increases, where the AII is less accurately estimated than the ADI, as the former is calculated from n (n − 1) terms and the latter is calculated from n elements only.

Empirical Application
Monte Carlo simulations results given in Section 4 support the consistency of our proposed estimation method.We are now in a position to re-investigate cross-country growth patterns.We want to evaluate the impact of a country's initial income, savings rate, population growth rate and openness, as well as neighbour countries' economic growth spillovers on a country's economic growth rate.We follow [4] in using a sample of 91 countries listed in [1], which is the Heston-Summers data taken from Penn World Table 6.1.Consider the following conditional convergence Solow growth model 2 : for i = 1, 2, ..., 91, where Y i is the i th country's average growth rate of real GDP per capita between 1960 and 1995, y60 i is the i th country's initial real GDP per capita in 1960, s i is the i th country's average saving rate, n i is the i th country's average growth rate of working-age population (ages between 15 and 64), lopen i is a scalar development index of a country defined as the logarithm of the i th country's average ratio of total imports plus exports over its real GDP over the period from 1960 to 1995, and Z ij is the great-circle distance between i th and j th countries' capitals 3 .We approximate g(•), m(•) and θ i (•), i = 1, 2, 3, 4, using the Laguerre polynomials with L n = 2.Moreover, a cross-validation selected bandwidth, h opt , is calculated as 0.5285.We obtained ρ( Ĝn ) = 0.041, which suggests a spatial stationarity in the data.The distribution of the estimated residuals is approximately normal as the q-q plot of estimated residuals is close to linear.The coefficient estimates, θ1 (•), θ2 (•), θ3 (•) and θ4 (•) are presented in Figure 1, where the solid lines with circles display the second-step estimates.The dashed lines represent the estimates of the spatially-augmented Solow growth model of [4] using the inverse power spatial weight function, which we include as a baseline.We interpret Figure 1 as follows.First, in Figure 1b, we see that there is a negative relation between the initial level of income and the economic growth rate, except for Mauritius, Hong Kong, Zambia, Cameroon and Singapore, which confirms a conditional β-convergence hypothesis.Moreover, we observe that θ2 (•) is increasing in openness, which, 2 Since the sample size is less than 100, we include only one spatially-lagged exogenous variable, M n ln(y60), to have better finite sample estimation accuracy.Moreover, the reason behind this choice is that spatial lag effects from the savings rate and the population growth rate were not found significant in Ertur and Koch [4]'s Table IV  We follow [4] in calculating the variable Z ij : where radius is taken as the Earth's radius and lat i and long i are the latitude and longitude for country i, respectively.however, results in a gradually declining degree of convergence.In addition, we see that the nonparametric model reveals slightly weaker conditional economic growth convergence as compared to the parametric model.Second, in Figure 1c, we see that θ3 (•) exhibits a positive, but not a monotonic, relation between the real investment rate and the real GDP per capita growth rate.Our estimate of the coefficient of the investment rate fluctuates as the trade openness of countries increases.For the economies with a trade openness higher than 15% of GDP, our result indicates that the nonparametric model sees stronger positive impact of the investment rate on the real GDP per capita growth rate than the parametric model does.Third, in Figure 1d, it is observed that the population growth rate has a negative impact on the real GDP per capita growth rate.For the countries whose trade openness ranges between 29% and 65%, our estimates for the coefficient of the population growth rate are relatively flat.Moreover, we note that the magnitude of the negative effect of the population growth rate is getting larger as the trade openness of countries increases, especially when the trade openness is over 65% of GDP.Overall, Figure 1 can be interpreted as the fact that an open economy suffers from higher negative impact of the population growth rate, but at the same time takes the advantage of high initial real GDP per capita.
Next, due to the cross-country interactions through spatial weights, the functional coefficient estimates have a different interpretation than the one obtained from the non-spatial model.In order to correctly interpret these estimates, we rewrite the estimated model in a reduced from as follows: Y = (I n − Ĝn ) −1 [ Mn ln(y60) + θ1 (lopen) + θ2 (lopen) ln(y60) + θ3 (lopen) ln (s) + θ4 (lopen) ln(n + 0.05) + Ũ], (20) where Ũ is the n × 1 vector of residuals.Then, from Equation ( 20), the marginal effects are given by the following n × n matrices: where we define diag{a} as an n × n diagonal matrix with the elements of an n × 1 vector, a, on the main diagonal.Following LeSage and Pace [10], we label the diagonal elements of each matrix given above as the direct impacts and off-diagonal elements as the indirect impacts.
In Table 2, we report the estimated average direct impact (ADI) and average indirect impact (AII) of the explanatory variables, where the latter can be easily defined as the difference between average total impact 4 and the average direct impact.Average direct and indirect impacts from the parametric model of [4] are denoted as ADI EK and AII EK , respectively.The interpretation of Table 2 is as follows.Firstly, we observe that a 1% increase in the real initial GDP per capita of an economy, holding other factors fixed, results in a decrease by 0.5% in its own real GDP per capita growth rate.However, this change increases the rest of the economies' economic growth rates by 0.01% on average due to the spatial dependence.From another point of view, a 1% increase in all of the regions'/nations' initial real GDP per capita speeds up this economy's real GDP per capita growth by 0.01%.This result indicates a positive spillover effect of the initial level of income.Secondly, a 1% increase in this economy's real investment rate increases its own real GDP growth rate by 2.08% on average.However, this change slows down the rest of the nations' real GDP growth by 0.08% on average.Thirdly, we see that a 1% increase in the population growth rate of this economy retards its own economic growth by 3.8% on average, but helps to improve the rate of economic growth of the rest of the countries by 0.16% on average.
Fourthly, when comparing our results to [4]'s results, we find that both the nonparametric and the parametric model give almost the same average direct effects of the initial per capita income, investment rate and population growth rate on the economic growth rate and that both models result in the same signs in the average direct and indirect effects.However, the AII values from the nonparametric model are much smaller than the results from the parametric model in absolute value, especially for the initial per capita income and the population growth rate.This is not surprising, as the parametric model assumes that all of the spatial weights take non-negative values, while our nonparametric spatial weights are estimated from the data without such a restriction.Although it is popular practice to assume non-negative spatial weights, this is an assumption imposed without support from econometric or economic theory.For example, trade treaties and monetary policies are both double-edged swords that may bring opposite impacts to different national economies, and non-negative spatial weights may not be able to capture the opposite interactions among different economies.As both [4]'s parametric SDM and our proposed semiparametric SDM approximate the unknown true relationships in their own best capacity, however, our model imposes less restrictions than the parametric SDM and is believed to bring a better fit to the data and more reliable inference.Although the numbers are different, both models give the same sign in estimated direct and indirect 4 As it is stated by LeSage and Pace [10] (p.37), average total impact can be expressed in two different ways, however, which give the same numerical results.The first viewpoint states an influence from a change in the initial real GDP per capita of an economy on all of the regions, while the second viewpoint states an impact of changes in the initial real GDP per capita of the entire economy on a region/nation.The confidence intervals are reported in Table 2.We see that the ADI values are statistically significantly different from zero at the 5% significance level.Moreover, we find that there is no significant effect on average from neighbouring countries per capita initial income, savings rate and population growth rate on economic growth rate of the country of interest.The same inference is obtained for the parametric model, except that the parametric model sees a significant average indirect impact of initial per capita income.Note that insignificant AIIs do not imply that the indirect impact from economy i on economy j is insignificant for all (i, j).
Figure 2 presents estimated spatial weighting functions.We plot both estimated spatial weighting functions, ĝ (•) and m (•), for the geographic distances ranging from zero to 20 in 100 km, as the estimated spatial weights in the absolute value have an average of 5.069 × 10 −7 and 3.893 × 10 −9 , respectively, when z > 20.In Figure 2a,b, firstly, we see negative spatial weights, which greatly contradicts traditional parametric spatial regression models, which often assume non-negative spatial weights.Negative spatial interactions are indeed common in practice, especially in social networks; see [52,53] for strategic interactions within the monetary policy committee of the Bank of England.Secondly, both spatial weight functions are not strictly monotonic and exhibit convexity among nations that are not very far apart, concavity among nations with moderately far distances and a zero spatial weight function among nations that are far away.Moreover, the spatial weight functions take bigger absolute values among nations with smaller distance apart and smaller absolute values among far-away nations, which implies a relatively larger economic interaction among nearby nations than among far-away nations.Lastly, in Figure 2b, we observe positive estimated spatial weights when the distance ranges between 0.229 and 1.894, which correspond to 22.9 and 189.4 km, respectively.For the distances greater than 189.4 km, negative spatial weights are getting closer to zero as the distance between two countries increases.As the turning point 1.894 is really small and occurs for countries with a small area, our results imply that spatial interaction is very strong and different between two nearby small countries with small areas than between two countries with longer distances when at least one country has a large area.

Conclusions
We employ a spatial Durbin model combined with the nonparametric spatial weighting functions, as well as the unknown functional coefficients to estimate the augmented Solow growth model with a sample of 91 countries over the period 1960 to 1995.We find a negative spatial lag effect of neighbouring country's GDP per capita growth rate and initial GDP per capita on the economic growth rate of country i.These effects are declining in magnitude as the geographical distance between the two countries increases.Finally, allowing coefficients as a function of trade openness of a country enables us to see the true country-specific effect of each determinant of economic growth.Moreover, we find significant average direct impact from each production factors.However, our findings show that the average indirect impact of these variables is insignificant at the 5% significance level.