Comparing estimation methods. econometrics
|
|
- Albert Carroll
- 5 years ago
- Views:
Transcription
1 for spatial econometrics Recent Advances in Spatial Econometrics (in honor of James LeSage), ERSA 2012 Roger Bivand Gianfranco Piras NHH Norwegian School of Economics Regional Research Institute at West Virginia University Thursday, 23 August 2012
2 Outline Background Recent advances in spatial econometrics model fitting techniques have made it more desirable to be able to compare results
3 Outline Background Recent advances in spatial econometrics model fitting techniques have made it more desirable to be able to compare results Results should correspond between implementations using different applications
4 Outline Background Recent advances in spatial econometrics model fitting techniques have made it more desirable to be able to compare results Results should correspond between implementations using different applications A broad range of model fitting techniques are provided by the contributed R packages for spatial econometrics
5 Outline Background Recent advances in spatial econometrics model fitting techniques have made it more desirable to be able to compare results Results should correspond between implementations using different applications A broad range of model fitting techniques are provided by the contributed R packages for spatial econometrics These model fitting techniques are associated with methods for estimating impacts and some tests, which will also be presented and compared
6 Background The use of spatial econometrics tools was widened by the ease with which methods and examples presented in Anselin (1988) could be reproduced using SpaceStat TM, written in Gauss TM
7 Background The use of spatial econometrics tools was widened by the ease with which methods and examples presented in Anselin (1988) could be reproduced using SpaceStat TM, written in Gauss TM It was rapidly complemented by the Spatial Econometrics toolbox for Matlab TM, provided as source code together with extensive documentation
8 Background The use of spatial econometrics tools was widened by the ease with which methods and examples presented in Anselin (1988) could be reproduced using SpaceStat TM, written in Gauss TM It was rapidly complemented by the Spatial Econometrics toolbox for Matlab TM, provided as source code together with extensive documentation A suite of commands for spatial data analysis for use with Stata TM was provided by Maurizio Pisati, and macros for Minitab TM and SAS TM were also made available
9 Background The use of spatial econometrics tools was widened by the ease with which methods and examples presented in Anselin (1988) could be reproduced using SpaceStat TM, written in Gauss TM It was rapidly complemented by the Spatial Econometrics toolbox for Matlab TM, provided as source code together with extensive documentation A suite of commands for spatial data analysis for use with Stata TM was provided by Maurizio Pisati, and macros for Minitab TM and SAS TM were also made available The thrust of SpaceStat TM has largely been taken over by GeoDa (Anselin et al. 2006), and more recently by OpenGeoDa
10 Today s software There is now much more software available for spatial econometrics
11 Today s software There is now much more software available for spatial econometrics Stata TM with sppack and Matlab TM with Spatial Econometrics Toolbox are mainstream programmes; the Matlab TM toolbox remains in the public domain, and has a community of contributors
12 Today s software There is now much more software available for spatial econometrics Stata TM with sppack and Matlab TM with Spatial Econometrics Toolbox are mainstream programmes; the Matlab TM toolbox remains in the public domain, and has a community of contributors OpenGeoDa and PySAL are open source, with code hosted on Google, binary versions for common platforms, and a community of users
13 Today s software There is now much more software available for spatial econometrics Stata TM with sppack and Matlab TM with Spatial Econometrics Toolbox are mainstream programmes; the Matlab TM toolbox remains in the public domain, and has a community of contributors OpenGeoDa and PySAL are open source, with code hosted on Google, binary versions for common platforms, and a community of users R with spdep, sphet, McSpatial and other contributed packages is open source, and the packages are cross-platform; the packages also have a community of users and developers
14 Why compare? In the spirit of Rey (2009), this comparison will attempt to examine some features of the implementation of functions for fitting spatial econometrics models
15 Why compare? In the spirit of Rey (2009), this comparison will attempt to examine some features of the implementation of functions for fitting spatial econometrics models Firstly, it may be useful to show which kinds of functions for creating spatial weights, for diagnostics, and for model fitting are available
16 Why compare? In the spirit of Rey (2009), this comparison will attempt to examine some features of the implementation of functions for fitting spatial econometrics models Firstly, it may be useful to show which kinds of functions for creating spatial weights, for diagnostics, and for model fitting are available Next, it is comforting when one can show that fitting the same model on the same data using different implementations gives the same results
17 Why compare? In the spirit of Rey (2009), this comparison will attempt to examine some features of the implementation of functions for fitting spatial econometrics models Firstly, it may be useful to show which kinds of functions for creating spatial weights, for diagnostics, and for model fitting are available Next, it is comforting when one can show that fitting the same model on the same data using different implementations gives the same results Finally, if the results are not the same, it is helpful to be able to show why they vary, possibly because of different design choices in implementation
18 Data set: Ward and Gleditsch 2008 The data set used for comparison here is taken from a Sage volume Spatial Regression Models by Ward and Gleditsch (2008), with political science data for 158 countries
19 Data set: Ward and Gleditsch 2008 The data set used for comparison here is taken from a Sage volume Spatial Regression Models by Ward and Gleditsch (2008), with political science data for 158 countries The model they explore is the relationship between democracy scores (POLITY IV indicators) and the logarithm of country GDP per capita in 2002
20 Data set: Ward and Gleditsch 2008 The data set used for comparison here is taken from a Sage volume Spatial Regression Models by Ward and Gleditsch (2008), with political science data for 158 countries The model they explore is the relationship between democracy scores (POLITY IV indicators) and the logarithm of country GDP per capita in 2002 They treat countries as neighbours with non-zero spatial weights if their borders are closer than 200km from each other; data and weights are available for download from their site
21 Data set: Ward and Gleditsch 2008 The data set used for comparison here is taken from a Sage volume Spatial Regression Models by Ward and Gleditsch (2008), with political science data for 158 countries The model they explore is the relationship between democracy scores (POLITY IV indicators) and the logarithm of country GDP per capita in 2002 They treat countries as neighbours with non-zero spatial weights if their borders are closer than 200km from each other; data and weights are available for download from their site They were among the first to examine the effects of feedback in spatial regression models
22 Spatial weights Creating spatial weights is a necessary step in using areal data, perhaps just to check that there is no remaining spatial patterning in residuals
23 Spatial weights Creating spatial weights is a necessary step in using areal data, perhaps just to check that there is no remaining spatial patterning in residuals The Matlab SE toolbox provides functions for creating nearest neighbour and triangulation contiguity weights, and for reading GAL files
24 Spatial weights Creating spatial weights is a necessary step in using areal data, perhaps just to check that there is no remaining spatial patterning in residuals The Matlab SE toolbox provides functions for creating nearest neighbour and triangulation contiguity weights, and for reading GAL files A number of functions are included in the R spdep package to create neighbours, and from them weights; GAL and GWT files may be read and written
25 Spatial weights Creating spatial weights is a necessary step in using areal data, perhaps just to check that there is no remaining spatial patterning in residuals The Matlab SE toolbox provides functions for creating nearest neighbour and triangulation contiguity weights, and for reading GAL files A number of functions are included in the R spdep package to create neighbours, and from them weights; GAL and GWT files may be read and written Weights may be constructed using functions in GeoDa, and in Pysal using Python programming; GAL and GWT files may be read and written
26 Spatial weights Creating spatial weights is a necessary step in using areal data, perhaps just to check that there is no remaining spatial patterning in residuals The Matlab SE toolbox provides functions for creating nearest neighbour and triangulation contiguity weights, and for reading GAL files A number of functions are included in the R spdep package to create neighbours, and from them weights; GAL and GWT files may be read and written Weights may be constructed using functions in GeoDa, and in Pysal using Python programming; GAL and GWT files may be read and written The Stata spmat command provides for the creation of a number of different kinds of weights, and for file import and export
27 The symptom residual autocorrelation A map of the residuals of a least squares regression may help us to see whether neighbouring observations appear to have residuals od similar value, indicating positive autocorrelation, or possibly dissimilar values, indicating negative autocorrelation: OLS residuals
28 Spatial lag and Durbin models The spatial lag model is not dissimilar to traditional time series models, with the autocorrelation process in the dependent variable controlled by the exogenous matrix of spatial weights W (Ord 1975): y = ρwy + Xβ + ε, where y is an (N 1) vector of observations on a dependent variable taken at each of N locations, X is an (N k) matrix of exogenous variables, β is an (k 1) vector of parameters, ε is an (N 1) vector of disturbances and ρ is a scalar spatial error parameter (but called λ in Stata)
29 Spatial lag and Durbin models The spatial lag model is not dissimilar to traditional time series models, with the autocorrelation process in the dependent variable controlled by the exogenous matrix of spatial weights W (Ord 1975): y = ρwy + Xβ + ε, where y is an (N 1) vector of observations on a dependent variable taken at each of N locations, X is an (N k) matrix of exogenous variables, β is an (k 1) vector of parameters, ε is an (N 1) vector of disturbances and ρ is a scalar spatial error parameter (but called λ in Stata) The spatial Durbin model is the lag model augmented by spatially lagged right hand side variables: y = ρwy + Xβ + WXθ + ε, where θ is an ((k 1) 1) vector of parameters where W is row-standardised, and a (k 1) vector otherwise
30 Spatial lag model log-likelihood function l(β, ρ, σ 2 ) = N 2 ln 2π N 2 ln σ2 + ln I ρw 1 2σ 2 [ y (I ρw) (I X(X X) 1 X )(I ρw)y ] and β = (X X) 1 (I ˆρW)y, where ˆρ is the ML estimate. Unlike the time series case, the logarithm of the determinant of the (N N) asymmetric matrix (I ρw) does not tend to zero with increasing sample size; it constrains the parameter values to their feasible range between the inverses of the smallest and largest eigenvalues of W, since for positive autocorrelation, as ρ 1, ln I ρw
31 ML spatial lag model results The SE toolbox uses a pre-computed grid of log determinant values, choosing the nearest rather than computing exactly at each call to the log likelihood function; this accounts for slightly different coefficient estimates compated to other implementations. There are two implementations of the spatial lag in R in spdep and McSpatial respectively: SE toolbox R spdep R McSpatial Stata spatreg OpenGeoDa (Intercept) log GDP pc ρ SE ρ σ LL
32 ML spatial lag model results The SE toolbox uses a pre-computed grid of log determinant values, choosing the nearest rather than computing exactly at each call to the log likelihood function; this accounts for slightly different coefficient estimates compated to other implementations. There are two implementations of the spatial lag in R in spdep and McSpatial respectively: SE toolbox R spdep R McSpatial Stata spatreg OpenGeoDa (Intercept) log GDP pc ρ SE ρ σ LL
33 ML spatial lag model differences I There are two major discrepancies in the table of results: the first is that the log-likelihood values at the optimimum differ between R McSpatial and the SE toolbox and the rest. SE toolbox R spdep R McSpatial Stata spatreg OpenGeoDa (Intercept) log GDP pc ρ SE ρ σ LL The reason appears to be that π in the log likelihood calculation is not multiplied by 2 in these cases, but is in the remainder. If we convert the R McSpatial value of by subtracting n 2 log(π) (line 65 in file McSpatial/R/sarml.R), and adding n 2 log(2π), we get Similarly, correcting the SE toolbox value of , we get (line 453 file spatial/sar models/sar.m). The same kind of difference appears in other reported SE toolbox log likelihood values
34 ML spatial lag model differences II The other discrepancy is in the coefficient standard errors, given here in full: SE toolbox R spdep R McSpatial Stata spatreg OpenGeoDa SE (Intercept) SE log GDP pc SE ρ We see that R spdep and OpenGeoDa agree, with the SE toolbox function very close the code in spatial/sar models/sar.m following line 202 indicates that for N 500, asymptotic calculations will be used based on Anselin, (1980, 1988). The same source is used in R spdep and OpenGeoDa (probably after line 861, Regression/smile2.cpp), but what about the others?
35 ML spatial lag model differences II The other discrepancy is in the coefficient standard errors, given here in full: SE toolbox R spdep R McSpatial Stata spatreg OpenGeoDa SE (Intercept) SE log GDP pc SE ρ We see that R spdep and OpenGeoDa agree, with the SE toolbox function very close the code in spatial/sar models/sar.m following line 202 indicates that for N 500, asymptotic calculations will be used based on Anselin, (1980, 1988). The same source is used in R spdep and OpenGeoDa (probably after line 861, Regression/smile2.cpp), but what about the others?
36 ML spatial lag model differences II more The standard errors reported by R McSpatial are taken from the Hessian returned by the optimization function nlm. The R spdep function lagsarlm can return a similar numerical Hessian, computed by as a finite difference Hessian by fdhess in nlme, or as the Hessian output by numerical optimization function optim, using by default a quasi-newton method due to Broyden, Fletcher, Goldfarb and Shanno (BFGS). Stata spreg ml can also use a bfgs technique, but the default is a modified Newton-Raphson method nr; spatreg uses optimization method lf for easy fitting of maximum likelihood models. Since the R and Stata BFGS standard errors agree, it is confirmed that Stata is reporting standard errors taken from the Hessian returned by the optimization function, rather than the analytical calculations even for small n. R nlm R fdhess R BFGS Stata BFGS Stata NR Stata lf SE (Intercept) SE log GDP pc SE ρ
37 ML spatial lag model differences II more The standard errors reported by R McSpatial are taken from the Hessian returned by the optimization function nlm. The R spdep function lagsarlm can return a similar numerical Hessian, computed by as a finite difference Hessian by fdhess in nlme, or as the Hessian output by numerical optimization function optim, using by default a quasi-newton method due to Broyden, Fletcher, Goldfarb and Shanno (BFGS). Stata spreg ml can also use a bfgs technique, but the default is a modified Newton-Raphson method nr; spatreg uses optimization method lf for easy fitting of maximum likelihood models. Since the R and Stata BFGS standard errors agree, it is confirmed that Stata is reporting standard errors taken from the Hessian returned by the optimization function, rather than the analytical calculations even for small n. R nlm R fdhess R BFGS Stata BFGS Stata NR Stata lf SE (Intercept) SE log GDP pc SE ρ
38 ML spatial lag model differences II more The standard errors reported by R McSpatial are taken from the Hessian returned by the optimization function nlm. The R spdep function lagsarlm can return a similar numerical Hessian, computed by as a finite difference Hessian by fdhess in nlme, or as the Hessian output by numerical optimization function optim, using by default a quasi-newton method due to Broyden, Fletcher, Goldfarb and Shanno (BFGS). Stata spreg ml can also use a bfgs technique, but the default is a modified Newton-Raphson method nr; spatreg uses optimization method lf for easy fitting of maximum likelihood models. Since the R and Stata BFGS standard errors agree, it is confirmed that Stata is reporting standard errors taken from the Hessian returned by the optimization function, rather than the analytical calculations even for small n. R nlm R fdhess R BFGS Stata BFGS Stata NR Stata lf SE (Intercept) SE log GDP pc SE ρ
39 Maximum likelihood fitting differences The differences identified in the spatial lag case follow through for the other model specifications examined here.
40 Maximum likelihood fitting differences The differences identified in the spatial lag case follow through for the other model specifications examined here. The Matlab SE toolbox uses a grid rather than a line search/optimization to fit the spatial coefficient(s), so they usually agree only for the first few digits.
41 Maximum likelihood fitting differences The differences identified in the spatial lag case follow through for the other model specifications examined here. The Matlab SE toolbox uses a grid rather than a line search/optimization to fit the spatial coefficient(s), so they usually agree only for the first few digits. The Matlab SE toolbox also reports a log likelihood value using π rather than 2π.
42 Maximum likelihood fitting differences The differences identified in the spatial lag case follow through for the other model specifications examined here. The Matlab SE toolbox uses a grid rather than a line search/optimization to fit the spatial coefficient(s), so they usually agree only for the first few digits. The Matlab SE toolbox also reports a log likelihood value using π rather than 2π. Stata (both spreg ml and spatreg) reports coefficient standard errors taken from the coefficient covariance matrix (Hessian) used in optimization, rather than analytical values reported for small n by R spdep functions, Matlab SE toolbox functions, and OpenGeoDa.
43 ML spatial Durbin model results Only R and the Matlab SE toolbox provide the spatial Durbin model directly, Stata and OpenGeoDa can fit it after creating the lagged right-hand side variable(s) by hand. Here, all applications agree except SE toolbox, again because of the gridded log determinant value causing the numerical optimization to exit before finding the exact optimum. This is justified because the SE toolbox is intended to provide for Bayesian methods, and in that context the distribution of values is more important than point estimates: SE toolbox R spdep Stata spatreg OpenGeoDa (Intercept) log GDP pc lag log GDP pc ρ SE ρ σ LL
44 Spatial error model There are a number of alternative forms of spatial regression models; here we will also consider the spatial error model (also known as the simultaneous autoregressive (SAR) model); the model may be written as (Ord 1975): y = Xβ + u, u = λwu + ε, where y is an (N 1) vector of observations on a dependent variable taken at each of N locations, X is an (N k) matrix of exogenous variables, β is an (k 1) vector of parameters, ε is an (N 1) vector of disturbances and λ is a scalar spatial error parameter (except by Kelejian and Prucha, who term it ρ),
45 Spatial error model There are a number of alternative forms of spatial regression models; here we will also consider the spatial error model (also known as the simultaneous autoregressive (SAR) model); the model may be written as (Ord 1975): y = Xβ + u, u = λwu + ε, where y is an (N 1) vector of observations on a dependent variable taken at each of N locations, X is an (N k) matrix of exogenous variables, β is an (k 1) vector of parameters, ε is an (N 1) vector of disturbances and λ is a scalar spatial error parameter (except by Kelejian and Prucha, who term it ρ), and u is a spatially autocorrelated disturbance vector with constant variance and covariance terms specified by a fixed spatial weights matrix and a single coefficient λ: u N(0, σ 2 (I λw) 1 (I λw ) 1 )
46 Spatial error model log-likelihood function The log-likelihood function for the spatial error model: l(β, λ, σ 2 ) = n 2 ln(2π) n 2 ln(σ2 ) + ln( I λw ) 1 2σ 2 [ (y Xβ) (I λw) (I λw)(y Xβ) ]
47 Spatial error model log-likelihood function The log-likelihood function for the spatial error model: l(β, λ, σ 2 ) = n 2 ln(2π) n 2 ln(σ2 ) + ln( I λw ) 1 2σ 2 [ (y Xβ) (I λw) (I λw)(y Xβ) ] As we can see, the problem is one of balancing the log determinant term ln( I λw ) against the sum of squares term. When λ approaches the ends of its feasible range, the log determinant term may swamp the sum of squares term
48 ML spatial error model results Once again, the line search with exact calculation of the log determinant for R, Stata and OpenGeoDa agrees fully. There are minor differences in the standard errors between R and OpenGeoDa on the one hand and Stata on the other, because of the use of analytical standard errors for small n in OpenGeoDa and R, with Stata using a numerical Hessian. The SE toolbox estimates differ somewhat, with λ truncated to its gridded value: SE toolbox R spdep Stata spatreg OpenGeoDa (Intercept) log GDP pc λ SE λ σ LL
49 ML general spatial model results The general model includes two spatial processes and faces identification issues: y = ρwy + Xβ + u, u = λwu + ε, It is however used in some analyses, often in GM estimators, so the ML version is useful for comparison, but does need care in selecting starting values for numerical optimization; here R spdep and Stata spreg ml agree in all but coefficient standard errors (spdep and SE toolbox use analytical calculations, spreg uses the optimization Hessian): SE toolbox R spdep Stata (Intercept) log GDP pc ρ SE ρ λ SE λ σ LL
50 ML general spatial model results The general model includes two spatial processes and faces identification issues: y = ρwy + Xβ + u, u = λwu + ε, It is however used in some analyses, often in GM estimators, so the ML version is useful for comparison, but does need care in selecting starting values for numerical optimization; here R spdep and Stata spreg ml agree in all but coefficient standard errors (spdep and SE toolbox use analytical calculations, spreg uses the optimization Hessian): SE toolbox R spdep Stata (Intercept) log GDP pc ρ SE ρ λ SE λ σ LL
51 ML general model standard error differences The R spdep function sacsarlm can return a similar numerical Hessian, computed by as a finite difference Hessian by fdhess in nlme. Stata spreg ml can use a bfgs technique, but the default is a modified Newton-Raphson method nr. Since the R finite difference values and Stata BFGS standard errors agree closely, it appears that Stata is reporting standard errors taken from the Hessian returned by the optimization function, rather than the analytical calculations even for small n. R asymptotic R fdhess Stata BFGS Stata NR SE (Intercept) SE log GDP pc SE ρ SE λ
52 ML general model standard error differences The R spdep function sacsarlm can return a similar numerical Hessian, computed by as a finite difference Hessian by fdhess in nlme. Stata spreg ml can use a bfgs technique, but the default is a modified Newton-Raphson method nr. Since the R finite difference values and Stata BFGS standard errors agree closely, it appears that Stata is reporting standard errors taken from the Hessian returned by the optimization function, rather than the analytical calculations even for small n. R asymptotic R fdhess Stata BFGS Stata NR SE (Intercept) SE log GDP pc SE ρ SE λ
53 The recent introduction of Stata TM and GeoDaSpace functions makes it helpful to compare them with SE toolbox and R functions
54 The recent introduction of Stata TM and GeoDaSpace functions makes it helpful to compare them with SE toolbox and R functions Within R, some functions have been contributed to spdep by Luc Anselin, and modified by the authors, and others are in sphet, which now uses the function wrapper spreg
55 The recent introduction of Stata TM and GeoDaSpace functions makes it helpful to compare them with SE toolbox and R functions Within R, some functions have been contributed to spdep by Luc Anselin, and modified by the authors, and others are in sphet, which now uses the function wrapper spreg The functions use different parts of the literature as bases for implementation, and the consequences of these choices will be made clear here
56 The recent introduction of Stata TM and GeoDaSpace functions makes it helpful to compare them with SE toolbox and R functions Within R, some functions have been contributed to spdep by Luc Anselin, and modified by the authors, and others are in sphet, which now uses the function wrapper spreg The functions use different parts of the literature as bases for implementation, and the consequences of these choices will be made clear here Once again, we examine spatial lag, spatial error, and general spatial models
57 GMM spatial lag models Using two stage least squares with Wy instrumented by [WX, WWX], all the functions yield the same coefficient estimates. In the two R and SE toolbox functions, the error variance is calculated as σ 2 = e e n k, while in the other two implementations is simply calculated as e e n : SE toolbox R spdep R sphet spreg Stata GeoDaSpace (Intercept) (2.4576) (2.4576) (2.4576) (2.4341) (2.4341) log GDP pc (0.3783) (0.3783) (0.3783) (0.3747) (0.3747) ρ (0.2283) (0.2283) (0.2283) (0.2261) (0.2261)
58 GMM spatial lag models Using two stage least squares with Wy instrumented by [WX, WWX], all the functions yield the same coefficient estimates. In the two R and SE toolbox functions, the error variance is calculated as σ 2 = e e n k, while in the other two implementations is simply calculated as e e n : SE toolbox R spdep R sphet spreg Stata GeoDaSpace (Intercept) (2.4576) (2.4576) (2.4576) (2.4341) (2.4341) log GDP pc (0.3783) (0.3783) (0.3783) (0.3747) (0.3747) ρ (0.2283) (0.2283) (0.2283) (0.2261) (0.2261)
59 The heteroskedastic error case I White standard errors may be calculated in most of the functions directly, in which the asymptotic VC matrix can be estimated consistently by the sandwich form: (Ẑ Ẑ) 1 (Ẑ ˆΣẐ)(Ẑ Ẑ) 1, where ˆΣ is a diagonal matrix whose elements are the ei 2 ; the results are the same in all cases: R spdep R sphet Stata GeoDaSpace (Intercept) log GDP pc ρ
60 The heteroskedastic error case II GeoDaSpace and sphet also implement the Kelejian and Prucha (2007) HAC estimator of the variance covariance matrix. Here we compare standard error estimates using a Triangular kernel with a variable bandwidth of the six nearest neighbours. The available options for the kernel function in R are the Epanechnikov, Triangular, Bisquare, Parzen, Tukey-Hanning and Quadratic Spectral. The options available in GeoDa space are the Uniform, Triangular, Epanechnikov, Quartic and Gaussian. GeoDa space only allows for the implementation of adaptive kernel. R sphet spreg GeoDaSpace (Intercept) log GDP pc ρ
61 GMM spatial error models In the GMM spatial error model, we depend on the first stage residuals, the implementation of the moment coment conditions, and the tuning of the optimiser finding the spatial parameter λ, as well as defintions for finding the standard error of λ. We see that there are three cases, the first for the SE toolbox, spdep, and GeoDaSpace-2, using the Kelejian and Prucha (1999) moment conditions (OLS first stage), and the standard error of λ from Pruch (2004): SE toolbox R spdep R sphet Stata GeoDaSpace-1 GeoDaSpace-2 (Intercept) (3.0403) (3.0403) (3.0319) (3.2603) (3.0788) (3.0403) log GDP pc (0.3767) (0.3767) (0.3758) (0.4072) (0.3796) (0.3767) λ (0.1820) (0.1820) (0.0661) (0.0668) (0.0660)
62 GMM spatial error models The next case is for the sphet and GeoDaSpace-1 implementations, using Drukker et al., which use a TSLS first stage. The results are close, and differences seem to come from numerical optimization. The final case is Stata, which appears to use a restricted instrument set work is continuing to establish why this is chosen: SE toolbox R spdep R sphet Stata GeoDaSpace-1 GeoDaSpace-2 (Intercept) (3.0403) (3.0403) (3.0319) (3.2603) (3.0788) (3.0403) log GDP pc (0.3767) (0.3767) (0.3758) (0.4072) (0.3796) (0.3767) λ (0.1820) (0.1820) (0.0661) (0.0668) (0.0660)
63 GMM spatial error models The next case is for the sphet and GeoDaSpace-1 implementations, using Drukker et al., which use a TSLS first stage. The results are close, and differences seem to come from numerical optimization. The final case is Stata, which appears to use a restricted instrument set work is continuing to establish why this is chosen: SE toolbox R spdep R sphet Stata GeoDaSpace-1 GeoDaSpace-2 (Intercept) (3.0403) (3.0403) (3.0319) (3.2603) (3.0788) (3.0403) log GDP pc (0.3767) (0.3767) (0.3758) (0.4072) (0.3796) (0.3767) λ (0.1820) (0.1820) (0.0661) (0.0668) (0.0660)
64 Heteroskedasticity Results from sphet and GeoDaSpace are quite similar and the very minor differences (in the estimated value of the spatial parameter) seem to be due to differences in the optimizers. This confirms our intuition on the error model with homoskedastic errors. The standard error results can be made closer by making the same implementation choices, which differ slightly with regard to simplifications. Stata differs as before: R sphet-spreg Stata GeoDaSpace (Intercept) (2.9878) (3.1150) (2.9888) log GDP pc (0.3613) (0.3664) (0.3614) ρ (0.0743) (0.0754) (0.0742)
65 GMM implementations of the general (SARAR) model There are various implementations of the GMM general model. Some of them are based on the Kelejian and Prucha (1999) moment conditions (SE toolbox, gstsls in spdep and GeoDaSpace-2), the others are based on the Drukker, Egger and Prucha moments conditions (sphet, Stata and GeoDaSpace-1) with big differences in λ (ML: ρ 0.78, λ 0.45): SE toolbox R spdep R sphet Stata GeoDaSpace-1 GeoDaSpace-2 (Intercept) (2.6631) (2.3185) (2.2101) (2.2101) (2.0800) (2.2963) log GDP pc (0.3556) (0.3600) (0.3314) (0.3314) (0.3255) (0.3566) ρ (0.2033) (0.2072) (0.1814) (0.1814) (0.1763) (0.2053) λ (0.0363) (0.2502) (0.2502) (0.2224)
66 The general model under heteroskedasticity Three implementation are available: one from sphet, one from Stata, and the one from GeoDaSpace. The results are the same for all of the implementations. Again, sphet and GeDaSpace have the option of performing step 1.c from Arraiz et al. Even in this case, the results match, and, therefore, are not reported. R sphet-spreg Stata GeoDaSpace (Intercept) (2.2119) (2.2119) (2.2119) log GDP pc (0.3527) (0.3527) (0.3527) ρ (0.1845) (0.1845) (0.1845) λ (0.2562) (0.2562) (0.2562)
67 Interpreting spatial lag, Durbin and general models It has emerged over time, however, that the spatial dependence in the parameter ρ feeds back
68 Interpreting spatial lag, Durbin and general models It has emerged over time, however, that the spatial dependence in the parameter ρ feeds back This feedback comes from the fact that the reduced form model is y = (I ρw) 1 Xβ + (I ρw) 1 ε
69 Interpreting spatial lag, Durbin and general models It has emerged over time, however, that the spatial dependence in the parameter ρ feeds back This feedback comes from the fact that the reduced form model is y = (I ρw) 1 Xβ + (I ρw) 1 ε In the spatial lag model, y i / x jr = ((I ρw) 1 Iβ r ) ij, where I is the N N identity matrix, and (I ρw) 1 is known to be dense
70 Interpreting spatial lag, Durbin and general models It has emerged over time, however, that the spatial dependence in the parameter ρ feeds back This feedback comes from the fact that the reduced form model is y = (I ρw) 1 Xβ + (I ρw) 1 ε In the spatial lag model, y i / x jr = ((I ρw) 1 Iβ r ) ij, where I is the N N identity matrix, and (I ρw) 1 is known to be dense In the spatial Durbin model, y i / x jr = ((I ρw) 1 Iβ r Wθ r ) ij
71 Implementing impact measures The awkward S r (W) = ((I ρw) 1 Iβ r ) matrix term needed to calculate impact measures for the lag model, and S r (W) = ((I ρw) 1 (Iβ r Wθ r )) for the spatial Durbin model, may be approximated using traces of powers of the spatial weights matrix as well as analytically
72 Implementing impact measures The awkward S r (W) = ((I ρw) 1 Iβ r ) matrix term needed to calculate impact measures for the lag model, and S r (W) = ((I ρw) 1 (Iβ r Wθ r )) for the spatial Durbin model, may be approximated using traces of powers of the spatial weights matrix as well as analytically The average direct impacts are represented by the sum of the diagonal elements of the matrix divided by N for each exogenous variable
73 Implementing impact measures The awkward S r (W) = ((I ρw) 1 Iβ r ) matrix term needed to calculate impact measures for the lag model, and S r (W) = ((I ρw) 1 (Iβ r Wθ r )) for the spatial Durbin model, may be approximated using traces of powers of the spatial weights matrix as well as analytically The average direct impacts are represented by the sum of the diagonal elements of the matrix divided by N for each exogenous variable The average total impacts are the sum of all matrix elements divided by N for each exogenous variable
74 Implementing impact measures The awkward S r (W) = ((I ρw) 1 Iβ r ) matrix term needed to calculate impact measures for the lag model, and S r (W) = ((I ρw) 1 (Iβ r Wθ r )) for the spatial Durbin model, may be approximated using traces of powers of the spatial weights matrix as well as analytically The average direct impacts are represented by the sum of the diagonal elements of the matrix divided by N for each exogenous variable The average total impacts are the sum of all matrix elements divided by N for each exogenous variable The average indirect impacts are the differences between the direct and total impact vectors
75 Total impacts Total impacts are defined as the sum of the elements of S r (W) divided by N. Here we only have a single variable, so the table of total impacts is simple. The SE toolbox implementation uses Monte Carlo simulation to provide a measure of the significance of the impacts, and by default reports the mean of the simulated values, which may differ a little from the computed value. SE toolbox (simulated) SE toolbox (computed) R spdep Stata Spatial lag Spatial Durbin General These are only headline values, and do not do justice to the possibilities for interpretation offered by this analytical advance.
76 Finding total impacts in Stata In Stata, we use the difference in predictions from the reduced form model when incrementing a chosen right-hand side variable, or even a single observation on that variable:. spreg ml democracy x, id(id) dlmat(w). predict y0. generate x_orig = x. quietly replace x = x + 1. predict y1. generate deltay = y1-y0. mean deltay. quietly replace x = x_orig We can do the same in R, using EXP to increment the variable in the scope of the objects: > EXP <- exp(0) > form <- formula(democracy ~ + log((gdp_2002/population) * + EXP)) > sldv.lag <- lagsarlm(form, data = sldv, + listw = lw) > p0 <- predict(sldv.lag, newdata = sldv, + listw = lw) > EXP <- exp(1) > p1 <- predict(sldv.lag, newdata = sldv, + listw = lw) > d <- p1 - p0 > mean(d) [1]
77 Spatial lag model impacts We use Monte Carlo methods to infer from the impacts, by drawing from the fitted model using the variance-covariance matrix of fitted coefficients: > set.seed(120823) > imp.lag <- impacts(sldv.lag, tr = trmat, R = 1999) > summary(imp.lag, short = TRUE, zstats = TRUE) Impact measures (lag, trace): Direct Indirect Total log((gdp_2002/population) * EXP) ======================================================== Simulation results (asymptotic variance matrix): ======================================================== Simulated z-values: Direct Indirect Total log((gdp_2002/population) * EXP) Simulated p-values: Direct Indirect Total log((gdp_2002/population) * EXP) > invirw <- invirw(lw, rho = coef(sldv.lag)[1]) > N <- nrow(sldv) > Sr <- invirw %*% (diag(n) * coef(sldv.lag)[3]) > c(direct = sum(diag(sr))/n, indirect = sum(sr)/n - sum(diag(sr))/n, + total = sum(sr)/n) direct indirect total
78 Spatial Durbin model impacts The spatial Durbin impacts are less significant than those of the spatial lag model: > EXP <- exp(0) > sldv.sd <- lagsarlm(form, data = sldv, listw = lw, type = "mixed") > imp.sd <- impacts(sldv.sd, tr = trmat, R = 1999) > summary(imp.sd, short = TRUE, zstats = TRUE) Impact measures (mixed, trace): Direct Indirect Total log((gdp_2002/population) * EXP) ======================================================== Simulation results (asymptotic variance matrix): ======================================================== Simulated z-values: Direct Indirect Total log((gdp_2002/population) * EXP) Simulated p-values: Direct Indirect Total log((gdp_2002/population) * EXP) > invirw <- invirw(lw, rho = coef(sldv.sd)[1]) > Sr <- invirw %*% ((diag(n) * coef(sldv.sd)[3]) + (W * coef(sldv.sd)[4])) > c(direct = sum(diag(sr))/n, indirect = sum(sr)/n - sum(diag(sr))/n, + total = sum(sr)/n) direct indirect total
79 General model impacts The general model impacts are similar to those of the spatial lag model: > EXP <- exp(0) > sldv.sac <- sacsarlm(form, data = sldv, listw = lw) > imp.sac <- impacts(sldv.sac, tr = trmat, R = 1999) > summary(imp.sac, short = TRUE, zstats = TRUE) Impact measures (sac, trace): Direct Indirect Total log((gdp_2002/population) * EXP) ======================================================== Simulation results (asymptotic variance matrix): ======================================================== Simulated z-values: Direct Indirect Total log((gdp_2002/population) * EXP) Simulated p-values: Direct Indirect Total log((gdp_2002/population) * EXP) > invirw <- invirw(lw, rho = coef(sldv.sac)[1]) > Sr <- invirw %*% ((diag(n) * coef(sldv.sac)[4])) > c(direct = sum(diag(sr))/n, indirect = sum(sr)/n - sum(diag(sr))/n, + total = sum(sr)/n) direct indirect total
80 Impact measures At present, there is no provision for measures of impact in OpenGeoDa or Pysal
81 Impact measures At present, there is no provision for measures of impact in OpenGeoDa or Pysal The total impact (emanating effect, equilibrium effect) can be calculated in Stata, but not broken down into direct and indirect
82 Impact measures At present, there is no provision for measures of impact in OpenGeoDa or Pysal The total impact (emanating effect, equilibrium effect) can be calculated in Stata, but not broken down into direct and indirect Only SE toolbox and R provide full support with Monte Carlo simulation for inference
83 Impact measures At present, there is no provision for measures of impact in OpenGeoDa or Pysal The total impact (emanating effect, equilibrium effect) can be calculated in Stata, but not broken down into direct and indirect Only SE toolbox and R provide full support with Monte Carlo simulation for inference They draw samples from the fitted model using the coefficient values and covariance matrix, and present summaries of the sample values
84 Distributions of general model impact measures Once one has the samples, it is possible to show how the distributions shift. In this case, the direct impacts lie further from zero than the coefficient, followed by indirect impacts even further from zero, with the total impacts shifted substantially beyond the shape of the distribution of the coefficient: Direct Indirect Total Coefficient Log GDP pc
85 Conclusions Background In this case, impact measures were not needed, because the LR and Hausman tests pointed to the spatial error specification; a recent paper by Pace and Zhu (2012) points to a enhanced error Durbin model as being of promise (it didn t help here).
86 Conclusions Background In this case, impact measures were not needed, because the LR and Hausman tests pointed to the spatial error specification; a recent paper by Pace and Zhu (2012) points to a enhanced error Durbin model as being of promise (it didn t help here). We have not considered Bayesian estimation methods, which will be covered in a separate study, where the SE toolbox is the only alternative so far, but an R GSoC project has been carried out in 2012
87 Conclusions Background In this case, impact measures were not needed, because the LR and Hausman tests pointed to the spatial error specification; a recent paper by Pace and Zhu (2012) points to a enhanced error Durbin model as being of promise (it didn t help here). We have not considered Bayesian estimation methods, which will be covered in a separate study, where the SE toolbox is the only alternative so far, but an R GSoC project has been carried out in 2012 The arrival of Stata s sppack opens up the alternatives a lot, but its spatial weights are dense or banded, limiting maximum likelihood estimation to smaller data sets
88 Conclusions Background In this case, impact measures were not needed, because the LR and Hausman tests pointed to the spatial error specification; a recent paper by Pace and Zhu (2012) points to a enhanced error Durbin model as being of promise (it didn t help here). We have not considered Bayesian estimation methods, which will be covered in a separate study, where the SE toolbox is the only alternative so far, but an R GSoC project has been carried out in 2012 The arrival of Stata s sppack opens up the alternatives a lot, but its spatial weights are dense or banded, limiting maximum likelihood estimation to smaller data sets Estimating models with maximum likelihood for large data sets is possible in the SE toolbox, OpenGeoDa and R using sparse matrix methods; GM models are not as limited by the size of data sets given care in avoiding handling n n matrices
Regional Research Institute
Regional Research Institute Working Paper Series Comparing Implementations of Estimation Methods for Spatial Econometrics By: Roger Bivand, Norwegian School of Economics and Gianfranco Piras, West Virginia
More informationJournal of Statistical Software
JSS Journal of Statistical Software January 2015, Volume 63, Issue 18. http://www.jstatsoft.org/ Comparing Implementations of Estimation Methods for Spatial Econometrics Roger Bivand Norwegian School of
More informationsplm: econometric analysis of spatial panel data
splm: econometric analysis of spatial panel data Giovanni Millo 1 Gianfranco Piras 2 1 Research Dept., Generali S.p.A. and DiSES, Univ. of Trieste 2 REAL, UIUC user! Conference Rennes, July 8th 2009 Introduction
More informationAnalyzing spatial autoregressive models using Stata
Analyzing spatial autoregressive models using Stata David M. Drukker StataCorp Summer North American Stata Users Group meeting July 24-25, 2008 Part of joint work with Ingmar Prucha and Harry Kelejian
More informationLecture 6: Hypothesis Testing
Lecture 6: Hypothesis Testing Mauricio Sarrias Universidad Católica del Norte November 6, 2017 1 Moran s I Statistic Mandatory Reading Moran s I based on Cliff and Ord (1972) Kelijan and Prucha (2001)
More informationSpatial Regression. 11. Spatial Two Stage Least Squares. Luc Anselin. Copyright 2017 by Luc Anselin, All Rights Reserved
Spatial Regression 11. Spatial Two Stage Least Squares Luc Anselin http://spatial.uchicago.edu 1 endogeneity and instruments spatial 2SLS best and optimal estimators HAC standard errors 2 Endogeneity and
More informationSpatial Regression. 9. Specification Tests (1) Luc Anselin. Copyright 2017 by Luc Anselin, All Rights Reserved
Spatial Regression 9. Specification Tests (1) Luc Anselin http://spatial.uchicago.edu 1 basic concepts types of tests Moran s I classic ML-based tests LM tests 2 Basic Concepts 3 The Logic of Specification
More informationGMM Estimation of Spatial Error Autocorrelation with and without Heteroskedasticity
GMM Estimation of Spatial Error Autocorrelation with and without Heteroskedasticity Luc Anselin July 14, 2011 1 Background This note documents the steps needed for an efficient GMM estimation of the regression
More informationSpatial Econometrics. Wykªad 6: Multi-source spatial models. Andrzej Torój. Institute of Econometrics Department of Applied Econometrics
Spatial Econometrics Wykªad 6: Multi-source spatial models (6) Spatial Econometrics 1 / 21 Outline 1 Multi-source models 2 SARAR model 3 SDM (Durbin) model 4 SDEM model 5 Exercises (6) Spatial Econometrics
More informationLecture 7: Spatial Econometric Modeling of Origin-Destination flows
Lecture 7: Spatial Econometric Modeling of Origin-Destination flows James P. LeSage Department of Economics University of Toledo Toledo, Ohio 43606 e-mail: jlesage@spatial-econometrics.com June 2005 The
More informationSpatial Tools for Econometric and Exploratory Analysis
Spatial Tools for Econometric and Exploratory Analysis Michael F. Goodchild University of California, Santa Barbara Luc Anselin University of Illinois at Urbana-Champaign http://csiss.org Outline A Quick
More informationsphet: Spatial Models with Heteroskedastic Innovations in R
sphet: Spatial Models with Heteroskedastic Innovations in R Gianfranco Piras Cornell University Abstract This introduction to the R package sphet is a (slightly) modified version of Piras (2010), published
More informationTesting Random Effects in Two-Way Spatial Panel Data Models
Testing Random Effects in Two-Way Spatial Panel Data Models Nicolas Debarsy May 27, 2010 Abstract This paper proposes an alternative testing procedure to the Hausman test statistic to help the applied
More informationW-BASED VS LATENT VARIABLES SPATIAL AUTOREGRESSIVE MODELS: EVIDENCE FROM MONTE CARLO SIMULATIONS
1 W-BASED VS LATENT VARIABLES SPATIAL AUTOREGRESSIVE MODELS: EVIDENCE FROM MONTE CARLO SIMULATIONS An Liu University of Groningen Henk Folmer University of Groningen Wageningen University Han Oud Radboud
More informationSpatial Econometrics
Spatial Econometrics Lecture 5: Single-source model of spatial regression. Combining GIS and regional analysis (5) Spatial Econometrics 1 / 47 Outline 1 Linear model vs SAR/SLM (Spatial Lag) Linear model
More informationLecture 3: Spatial Analysis with Stata
Lecture 3: Spatial Analysis with Stata Paul A. Raschky, University of St Gallen, October 2017 1 / 41 Today s Lecture 1. Importing Spatial Data 2. Spatial Autocorrelation 2.1 Spatial Weight Matrix 3. Spatial
More informationGreene, Econometric Analysis (6th ed, 2008)
EC771: Econometrics, Spring 2010 Greene, Econometric Analysis (6th ed, 2008) Chapter 17: Maximum Likelihood Estimation The preferred estimator in a wide variety of econometric settings is that derived
More informationAfter "Raising the Bar'': applied maximum likelihood estimation of families of models in spatial econometrics ( )
Estadística Española Volumen 54, número 177 / 2012, pp. 71-88 After "Raising the Bar'': applied maximum likelihood estimation of families of models in spatial econometrics ( ) Roger Bivand NHH Norwegian
More informationA SPATIAL CLIFF-ORD-TYPE MODEL WITH HETEROSKEDASTIC INNOVATIONS: SMALL AND LARGE SAMPLE RESULTS
JOURNAL OF REGIONAL SCIENCE, VOL. 50, NO. 2, 2010, pp. 592 614 A SPATIAL CLIFF-ORD-TYPE MODEL WITH HETEROSKEDASTIC INNOVATIONS: SMALL AND LARGE SAMPLE RESULTS Irani Arraiz Inter-American Development Bank,
More informationJournal of Statistical Software
JSS Journal of Statistical Software April 2012, Volume 47, Issue 1. http://www.jstatsoft.org/ splm: Spatial Panel Data Models in R Giovanni Millo Generali SpA Gianfranco Piras West Virginia University
More informationRegional Science and Urban Economics
Regional Science and Urban Economics 46 (2014) 103 115 Contents lists available at ScienceDirect Regional Science and Urban Economics journal homepage: www.elsevier.com/locate/regec On the finite sample
More informationOutline. Overview of Issues. Spatial Regression. Luc Anselin
Spatial Regression Luc Anselin University of Illinois, Urbana-Champaign http://www.spacestat.com Outline Overview of Issues Spatial Regression Specifications Space-Time Models Spatial Latent Variable Models
More informationComputing the Jacobian in Gaussian spatial models: an illustrated comparison of available methods
Computing the Jacobian in Gaussian spatial models: an illustrated comparison of available methods 1 1 1 1 1 1 Abstract When fitting spatial regression models by maximum likelihood using spatial weights
More informationVisualize and interactively design weight matrices
Visualize and interactively design weight matrices Angelos Mimis *1 1 Department of Economic and Regional Development, Panteion University of Athens, Greece Tel.: +30 6936670414 October 29, 2014 Summary
More informationApplied Spatial Econometrics. Professor Bernard Fingleton Director of Research Department of Land Economy University of Cambridge
Applied Spatial Econometrics Professor Bernard Fingleton Director of Research Department of Land Economy University of Cambridge Topics to discuss Regression and spatial dependence Residual Spatial autocorrelation
More informationPOLI 8501 Introduction to Maximum Likelihood Estimation
POLI 8501 Introduction to Maximum Likelihood Estimation Maximum Likelihood Intuition Consider a model that looks like this: Y i N(µ, σ 2 ) So: E(Y ) = µ V ar(y ) = σ 2 Suppose you have some data on Y,
More informationInterpreting dynamic space-time panel data models
Interpreting dynamic space-time panel data models Nicolas Debarsy CERPE, University of Namur, Rempart de la vierge, 8, 5000 Namur, Belgium E-mail: ndebarsy@fundp.ac.be Cem Ertur LEO, Université d Orléans
More informationFinite Sample Properties of Moran s I Test for Spatial Autocorrelation in Probit and Tobit Models - Empirical Evidence
Finite Sample Properties of Moran s I Test for Spatial Autocorrelation in Probit and Tobit Models - Empirical Evidence Pedro V. Amaral and Luc Anselin 2011 Working Paper Number 07 Finite Sample Properties
More informationGeographically weighted regression approach for origin-destination flows
Geographically weighted regression approach for origin-destination flows Kazuki Tamesue 1 and Morito Tsutsumi 2 1 Graduate School of Information and Engineering, University of Tsukuba 1-1-1 Tennodai, Tsukuba,
More informationHETEROSKEDASTICITY, TEMPORAL AND SPATIAL CORRELATION MATTER
ACTA UNIVERSITATIS AGRICULTURAE ET SILVICULTURAE MENDELIANAE BRUNENSIS Volume LXI 239 Number 7, 2013 http://dx.doi.org/10.11118/actaun201361072151 HETEROSKEDASTICITY, TEMPORAL AND SPATIAL CORRELATION MATTER
More informationSPACE Workshop NSF NCGIA CSISS UCGIS SDSU. Aldstadt, Getis, Jankowski, Rey, Weeks SDSU F. Goodchild, M. Goodchild, Janelle, Rebich UCSB
SPACE Workshop NSF NCGIA CSISS UCGIS SDSU Aldstadt, Getis, Jankowski, Rey, Weeks SDSU F. Goodchild, M. Goodchild, Janelle, Rebich UCSB August 2-8, 2004 San Diego State University Some Examples of Spatial
More informationSpatial Relationships in Rural Land Markets with Emphasis on a Flexible. Weights Matrix
Spatial Relationships in Rural Land Markets with Emphasis on a Flexible Weights Matrix Patricia Soto, Lonnie Vandeveer, and Steve Henning Department of Agricultural Economics and Agribusiness Louisiana
More informationSpatial inference. Spatial inference. Accounting for spatial correlation. Multivariate normal distributions
Spatial inference I will start with a simple model, using species diversity data Strong spatial dependence, Î = 0.79 what is the mean diversity? How precise is our estimate? Sampling discussion: The 64
More information1 Overview. 2 Data Files. 3 Estimation Programs
1 Overview The programs made available on this web page are sample programs for the computation of estimators introduced in Kelejian and Prucha (1999). In particular we provide sample programs for the
More informationSIMULATION AND APPLICATION OF THE SPATIAL AUTOREGRESSIVE GEOGRAPHICALLY WEIGHTED REGRESSION MODEL (SAR-GWR)
SIMULATION AND APPLICATION OF THE SPATIAL AUTOREGRESSIVE GEOGRAPHICALLY WEIGHTED REGRESSION MODEL (SAR-GWR) I. Gede Nyoman Mindra Jaya 1, Budi Nurani Ruchjana 2, Bertho Tantular 1, Zulhanif 1 and Yudhie
More informationthe error term could vary over the observations, in ways that are related
Heteroskedasticity We now consider the implications of relaxing the assumption that the conditional variance Var(u i x i ) = σ 2 is common to all observations i = 1,..., n In many applications, we may
More informationRAO s SCORE TEST IN SPATIAL ECONOMETRICS
RAO s SCORE TEST IN SPATIAL ECONOMETRICS Luc Anselin Regional Economics Applications Laboratory (REAL) and Department of Agricultural and Consumer Economics University of Illinois, Urbana-Champaign Urbana,
More informationNo
28 1 Vol. 28 No. 1 2015 1 Research of Finance and Education Jan. 2015 1 1 2 1. 200241 2. 9700AV 1998-2010 30 F830 A 2095-0098 2015 01-0043 - 010 Schumpeter 1911 2014-09 - 19 No. 40671074 2011 3010 1987-1983
More informationSpatial Regression Models: Identification strategy using STATA TATIANE MENEZES PIMES/UFPE
Spatial Regression Models: Identification strategy using STATA TATIANE MENEZES PIMES/UFPE Intruduction Spatial regression models are usually intended to estimate parameters related to the interaction of
More informationA SPATIAL ANALYSIS OF A RURAL LAND MARKET USING ALTERNATIVE SPATIAL WEIGHT MATRICES
A Spatial Analysis of a Rural Land Market Using Alternative Spatial Weight Matrices A SPATIAL ANALYSIS OF A RURAL LAND MARKET USING ALTERNATIVE SPATIAL WEIGHT MATRICES Patricia Soto, Louisiana State University
More informationUsing AMOEBA to Create a Spatial Weights Matrix and Identify Spatial Clusters, and a Comparison to Other Clustering Algorithms
Using AMOEBA to Create a Spatial Weights Matrix and Identify Spatial Clusters, and a Comparison to Other Clustering Algorithms Arthur Getis* and Jared Aldstadt** *San Diego State University **SDSU/UCSB
More informationLecture 4: Maximum Likelihood Estimator of SLM and SEM
Lecture 4: Maximum Likelihood Estimator of SLM and SEM Mauricio Sarrias Universidad Católica del Norte October 16, 2017 1 Spatial Lag Model Mandatory Reading Log-Likelihood function The score function
More informationSpatial Autocorrelation and Interactions between Surface Temperature Trends and Socioeconomic Changes
Spatial Autocorrelation and Interactions between Surface Temperature Trends and Socioeconomic Changes Ross McKitrick Department of Economics University of Guelph December, 00 1 1 1 1 Spatial Autocorrelation
More informationGMM estimation of spatial panels
MRA Munich ersonal ReEc Archive GMM estimation of spatial panels Francesco Moscone and Elisa Tosetti Brunel University 7. April 009 Online at http://mpra.ub.uni-muenchen.de/637/ MRA aper No. 637, posted
More informationarxiv: v1 [stat.co] 3 Mar 2017
Estimating Spatial Econometrics Models with Integrated Nested Laplace Approximation Virgilio Gómez-Rubio a,, Roger S. Bivand b, Håvard Rue c arxiv:1703.01273v1 [stat.co] 3 Mar 2017 a Department of Mathematics,
More informationChapter 2 Linear Spatial Dependence Models for Cross-Section Data
Chapter 2 Linear Spatial Dependence Models for Cross-Section Data Abstract This chapter gives an overview of all linear spatial econometric models with different combinations of interaction effects that
More informationA command for estimating spatial-autoregressive models with spatial-autoregressive disturbances and additional endogenous variables
The Stata Journal (2013) 13, Number 2, pp. 287 301 A command for estimating spatial-autoregressive models with spatial-autoregressive disturbances and additional endogenous variables David M. Drukker StataCorp
More informationAn estimate of the long-run covariance matrix, Ω, is necessary to calculate asymptotic
Chapter 6 ESTIMATION OF THE LONG-RUN COVARIANCE MATRIX An estimate of the long-run covariance matrix, Ω, is necessary to calculate asymptotic standard errors for the OLS and linear IV estimators presented
More informationIntroduction to PySAL and Web Based Spatial Statistics
1 Introduction to PySAL and Web Based Spatial Statistics Myung-Hwa Hwang GeoDa Center for Geospatial Analysis and Computation School of Geographical Sciences and Urban Planning Arizona State University
More informationHeteroskedasticity in Panel Data
Essex Summer School in Social Science Data Analysis Panel Data Analysis for Comparative Research Heteroskedasticity in Panel Data Christopher Adolph Department of Political Science and Center for Statistics
More informationHeteroskedasticity in Panel Data
Essex Summer School in Social Science Data Analysis Panel Data Analysis for Comparative Research Heteroskedasticity in Panel Data Christopher Adolph Department of Political Science and Center for Statistics
More informationHeteroskedasticity. We now consider the implications of relaxing the assumption that the conditional
Heteroskedasticity We now consider the implications of relaxing the assumption that the conditional variance V (u i x i ) = σ 2 is common to all observations i = 1,..., In many applications, we may suspect
More informationMotivation Non-linear Rational Expectations The Permanent Income Hypothesis The Log of Gravity Non-linear IV Estimation Summary.
Econometrics I Department of Economics Universidad Carlos III de Madrid Master in Industrial Economics and Markets Outline Motivation 1 Motivation 2 3 4 5 Motivation Hansen's contributions GMM was developed
More informationSpatial Regression. 6. Specification Spatial Heterogeneity. Luc Anselin.
Spatial Regression 6. Specification Spatial Heterogeneity Luc Anselin http://spatial.uchicago.edu 1 homogeneity and heterogeneity spatial regimes spatially varying coefficients spatial random effects 2
More informationISSN Article
Econometrics 2014, 2, 217-249; doi:10.3390/econometrics2040217 OPEN ACCESS econometrics ISSN 2225-1146 www.mdpi.com/journal/econometrics Article The Biggest Myth in Spatial Econometrics James P. LeSage
More informationProceedings of the 8th WSEAS International Conference on APPLIED MATHEMATICS, Tenerife, Spain, December 16-18, 2005 (pp )
Proceedings of the 8th WSEAS International Conference on APPLIED MATHEMATICS, Tenerife, Spain, December 16-18, 005 (pp159-164) Spatial Econometrics Modeling of Poverty FERDINAND J. PARAGUAS 1 AND ANTON
More informationChristopher Dougherty London School of Economics and Political Science
Introduction to Econometrics FIFTH EDITION Christopher Dougherty London School of Economics and Political Science OXFORD UNIVERSITY PRESS Contents INTRODU CTION 1 Why study econometrics? 1 Aim of this
More informationMaximum Likelihood (ML) Estimation
Econometrics 2 Fall 2004 Maximum Likelihood (ML) Estimation Heino Bohn Nielsen 1of32 Outline of the Lecture (1) Introduction. (2) ML estimation defined. (3) ExampleI:Binomialtrials. (4) Example II: Linear
More informationChoice of Spectral Density Estimator in Ng-Perron Test: Comparative Analysis
MPRA Munich Personal RePEc Archive Choice of Spectral Density Estimator in Ng-Perron Test: Comparative Analysis Muhammad Irfan Malik and Atiq-ur- Rehman International Institute of Islamic Economics, International
More informationInstrumental Variables/Method of
Instrumental Variables/Method of 80 Moments Estimation Ingmar R. Prucha Contents 80.1 Introduction... 1597 80.2 A Primer on GMM Estimation... 1599 80.2.1 Model Specification and Moment Conditions... 1599
More informationLuc Anselin and Nancy Lozano-Gracia
Errors in variables and spatial effects in hedonic house price models of ambient air quality Luc Anselin and Nancy Lozano-Gracia Presented by Julia Beckhusen and Kosuke Tamura February 29, 2008 AGEC 691T:
More informationCreating and Managing a W Matrix
Creating and Managing a W Matrix Carlos Hurtado Department of Economics University of Illinois at Urbana-Champaign hrtdmrt2@illinois.edu Junel 22th, 2016 C. Hurtado (UIUC - Economics) Spatial Econometrics
More informationCOLUMN. Spatial Analysis in R: Part 2 Performing spatial regression modeling in R with ACS data
Spatial Demography 2013 1(2): 219-226 http://spatialdemography.org OPEN ACCESS via Creative Commons 3.0 ISSN 2164-7070 (online) COLUMN Spatial Analysis in R: Part 2 Performing spatial regression modeling
More informationA Spatial Cliff-Ord-type Model with Heteroskedastic Innovations: Small and Large Sample Results 1
A Spatial Cliff-Ord-type Model with Heteroskedastic Innovations: Small and Large Sample Results 1 Irani Arraiz 2,DavidM.Drukker 3, Harry H. Kelejian 4 and Ingmar R. Prucha 5 August 21, 2007 1 Our thanks
More informationThe exact bias of S 2 in linear panel regressions with spatial autocorrelation SFB 823. Discussion Paper. Christoph Hanck, Walter Krämer
SFB 83 The exact bias of S in linear panel regressions with spatial autocorrelation Discussion Paper Christoph Hanck, Walter Krämer Nr. 8/00 The exact bias of S in linear panel regressions with spatial
More informationLecture 2: Spatial Models
Lecture 2: Spatial Models Mauricio Sarrias Universidad Católica del Norte August 28, 2017 1 Taxonomy of Models Mandatory Reading Spatial Lag Model Reduced Form and Parameter Space Expected Value and Variance
More informationSpatial Autocorrelation (2) Spatial Weights
Spatial Autocorrelation (2) Spatial Weights Luc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign http://sal.agecon.uiuc.edu Outline
More informationSpatial Effects and Externalities
Spatial Effects and Externalities Philip A. Viton November 5, Philip A. Viton CRP 66 Spatial () Externalities November 5, / 5 Introduction If you do certain things to your property for example, paint your
More informationRecent Advances in the Field of Trade Theory and Policy Analysis Using Micro-Level Data
Recent Advances in the Field of Trade Theory and Policy Analysis Using Micro-Level Data July 2012 Bangkok, Thailand Cosimo Beverelli (World Trade Organization) 1 Content a) Classical regression model b)
More informationSpatial Econometrics The basics
Spatial Econometrics The basics Interaction effects spatial econometric models Direct and indirect (spatial spillover) effects Bootlegging effect Basic numerical findings Identification issues Estimation
More informationNetwork data in regression framework
13-14 July 2009 University of Salerno (Italy) Network data in regression framework Maria ProsperinaVitale Department of Economics and Statistics University of Salerno (Italy) mvitale@unisa.it - Theoretical
More informationPitfalls in higher order model extensions of basic spatial regression methodology
Pitfalls in higher order model extensions of basic spatial regression methodology James P. LeSage Department of Finance and Economics Texas State University San Marcos San Marcos, TX 78666 jlesage@spatial-econometrics.com
More informationIntroduction to Spatial Statistics and Modeling for Regional Analysis
Introduction to Spatial Statistics and Modeling for Regional Analysis Dr. Xinyue Ye, Assistant Professor Center for Regional Development (Department of Commerce EDA University Center) & School of Earth,
More informationAgricultural and Applied Economics 637 Applied Econometrics II
Agricultural and Applied Economics 637 Applied Econometrics II Assignment 1 Review of GLS Heteroskedasity and Autocorrelation (Due: Feb. 4, 2011) In this assignment you are asked to develop relatively
More informationLinear models. Linear models are computationally convenient and remain widely used in. applied econometric research
Linear models Linear models are computationally convenient and remain widely used in applied econometric research Our main focus in these lectures will be on single equation linear models of the form y
More informationA Course in Applied Econometrics Lecture 7: Cluster Sampling. Jeff Wooldridge IRP Lectures, UW Madison, August 2008
A Course in Applied Econometrics Lecture 7: Cluster Sampling Jeff Wooldridge IRP Lectures, UW Madison, August 2008 1. The Linear Model with Cluster Effects 2. Estimation with a Small Number of roups and
More informationOmitted Variable Biases of OLS and Spatial Lag Models
Omitted Variable Biases of OLS and Spatial Lag Models R. Kelley Pace and James P. LeSage 1 Introduction Numerous authors have suggested that omitted variables affect spatial regression methods less than
More information1 Motivation for Instrumental Variable (IV) Regression
ECON 370: IV & 2SLS 1 Instrumental Variables Estimation and Two Stage Least Squares Econometric Methods, ECON 370 Let s get back to the thiking in terms of cross sectional (or pooled cross sectional) data
More informationSpatial Filtering with EViews and MATLAB
AUSTRIAN JOURNAL OF STATISTICS Volume 36 (2007), Number 1, 17 26 Spatial Filtering with EViews and MATLAB Robert Ferstl Vienna University of Economics and Business Administration Abstract: This article
More informationGMM Estimation of the Spatial Autoregressive Model in a System of Interrelated Networks
GMM Estimation of the Spatial Autoregressive Model in a System of Interrelated etworks Yan Bao May, 2010 1 Introduction In this paper, I extend the generalized method of moments framework based on linear
More informationA correlated randon effects spatial Durbin model
A correlated randon effects spatial Durbin model K Miranda, O Martínez-Ibañez, and M Manjón-Antolín QURE-CREIP Department of Economics, Rovira i Virgili University Abstract We consider a correlated random
More informationSpatial Interdependence and Instrumental Variable Models
Spatial Interdependence and Instrumental Variable Models Timm Betz Scott J. Cook Florian M. Hollenbach Conditionally accepted, Political Science Research and Methods Abstract Instrumental variable (IV)
More informationSingle Equation Linear GMM with Serially Correlated Moment Conditions
Single Equation Linear GMM with Serially Correlated Moment Conditions Eric Zivot October 28, 2009 Univariate Time Series Let {y t } be an ergodic-stationary time series with E[y t ]=μ and var(y t )
More informationP1: JYD /... CB495-08Drv CB495/Train KEY BOARDED March 24, :7 Char Count= 0 Part II Estimation 183
Part II Estimation 8 Numerical Maximization 8.1 Motivation Most estimation involves maximization of some function, such as the likelihood function, the simulated likelihood function, or squared moment
More informationShort T Panels - Review
Short T Panels - Review We have looked at methods for estimating parameters on time-varying explanatory variables consistently in panels with many cross-section observation units but a small number of
More informationInstrumental variables estimation using heteroskedasticity-based instruments
Instrumental variables estimation using heteroskedasticity-based instruments Christopher F Baum, Arthur Lewbel, Mark E Schaffer, Oleksandr Talavera Boston College/DIW Berlin, Boston College, Heriot Watt
More informationSpatial Panel Data Analysis
Spatial Panel Data Analysis J.Paul Elhorst Faculty of Economics and Business, University of Groningen, the Netherlands Elhorst, J.P. (2017) Spatial Panel Data Analysis. In: Shekhar S., Xiong H., Zhou X.
More informationLecture 7 Autoregressive Processes in Space
Lecture 7 Autoregressive Processes in Space Dennis Sun Stanford University Stats 253 July 8, 2015 1 Last Time 2 Autoregressive Processes in Space 3 Estimating Parameters 4 Testing for Spatial Autocorrelation
More informationG. S. Maddala Kajal Lahiri. WILEY A John Wiley and Sons, Ltd., Publication
G. S. Maddala Kajal Lahiri WILEY A John Wiley and Sons, Ltd., Publication TEMT Foreword Preface to the Fourth Edition xvii xix Part I Introduction and the Linear Regression Model 1 CHAPTER 1 What is Econometrics?
More informationTitle. Description. var intro Introduction to vector autoregressive models
Title var intro Introduction to vector autoregressive models Description Stata has a suite of commands for fitting, forecasting, interpreting, and performing inference on vector autoregressive (VAR) models
More informationGravity Models, PPML Estimation and the Bias of the Robust Standard Errors
Gravity Models, PPML Estimation and the Bias of the Robust Standard Errors Michael Pfaffermayr August 23, 2018 Abstract In gravity models with exporter and importer dummies the robust standard errors of
More informationEconometrics. 9) Heteroscedasticity and autocorrelation
30C00200 Econometrics 9) Heteroscedasticity and autocorrelation Timo Kuosmanen Professor, Ph.D. http://nomepre.net/index.php/timokuosmanen Today s topics Heteroscedasticity Possible causes Testing for
More informationInstrumental Variables and GMM: Estimation and Testing. Steven Stillman, New Zealand Department of Labour
Instrumental Variables and GMM: Estimation and Testing Christopher F Baum, Boston College Mark E. Schaffer, Heriot Watt University Steven Stillman, New Zealand Department of Labour March 2003 Stata Journal,
More informationxtdpdqml: Quasi-maximum likelihood estimation of linear dynamic short-t panel data models
xtdpdqml: Quasi-maximum likelihood estimation of linear dynamic short-t panel data models Sebastian Kripfganz University of Exeter Business School, Department of Economics, Exeter, UK UK Stata Users Group
More informationSpatial Regression. 13. Spatial Panels (1) Luc Anselin. Copyright 2017 by Luc Anselin, All Rights Reserved
Spatial Regression 13. Spatial Panels (1) Luc Anselin http://spatial.uchicago.edu 1 basic concepts dynamic panels pooled spatial panels 2 Basic Concepts 3 Data Structures 4 Two-Dimensional Data cross-section/space
More informationReading Assignment. Distributed Lag and Autoregressive Models. Chapter 17. Kennedy: Chapters 10 and 13. AREC-ECON 535 Lec G 1
Reading Assignment Distributed Lag and Autoregressive Models Chapter 17. Kennedy: Chapters 10 and 13. AREC-ECON 535 Lec G 1 Distributed Lag and Autoregressive Models Distributed lag model: y t = α + β
More informationAfter Raising the Bar : applied maximum likelihood estimation of families of models in spatial econometrics
INSTITUTT FOR SAMFUNNSØKONOMI DEPARTMENT OF ECONOMICS SAM 22 2011 ISSN: 0804-6824 November 2011 Discussion paper After Raising the Bar : applied maximum likelihood estimation of families of models in spatial
More informationMax. Likelihood Estimation. Outline. Econometrics II. Ricardo Mora. Notes. Notes
Maximum Likelihood Estimation Econometrics II Department of Economics Universidad Carlos III de Madrid Máster Universitario en Desarrollo y Crecimiento Económico Outline 1 3 4 General Approaches to Parameter
More informationIntroduction to Eco n o m et rics
2008 AGI-Information Management Consultants May be used for personal purporses only or by libraries associated to dandelon.com network. Introduction to Eco n o m et rics Third Edition G.S. Maddala Formerly
More informationON THE NEGATION OF THE UNIFORMITY OF SPACE RESEARCH ANNOUNCEMENT
MIDWEST STUDENT SUMMIT ON SPACE, HEALTH AND POPULATION ECONOMICS APRIL 18-19, 2007 PURDUE UNIVERSITY ON THE NEGATION OF THE UNIFORMITY OF SPACE RESEARCH ANNOUNCEMENT Benoit Delbecq Agricultural Economics
More information