Files
crankshaft/doc/21_gwr.md
T

5.9 KiB
Raw Blame History

Regression

Geographically weighted regression

Can currently estimate Gaussian, Poisson, and logistic models (built on a GLM framework). GWR object prepares model input. Fit method performs estimation and returns a GWR Results object.

Arguments

Name Type Description
subquery text SQL query that expose the data to be analyzed (e.g., SELECT * FROM regression_inputs). This query must have the geometry column name (see the optional geom_col for default), the id column name (see id_col), dependent and independent column names.
dep_var text name of the dependent variable in the regression model
ind_vars text[] Text array of independent used in the model to describe the dependent variable
bw (optional) numeric bandwidth value consisting of either a distance or N nearest neighbors. Defaults to calculate an optimal bandwidth.
fixed (optional) boolean True for distance based kernel function and False for adaptive (nearest neighbor) kernel function (default). Defaults to false.
kernel text Type of kernel function used to weight observations. One of gaussian, bisquare (default), or exponential.

Returns

Column Name Type Description
coeffs JSON JSON object with parameter estimates for each of the dependent variables. The keys of the JSON object are the dependent variables, with values corresponding to the parameter estimate.
stand_errs JSON Standard errors for each of the dependent variables. The keys of the JSON object are the dependent variables, with values corresponding to the respective standard errors.
t_vals JSON T-values for each of the dependent variables. The keys of the JSON object are the dependent variable names, with values corresponding to the respective t-value.
predicted numeric predicted value of y
residuals numeric residuals of the response
r_squared numeric R-squared for the parameter fit
bandwidth numeric bandwidth value consisting of either a distance or N nearest neighbors
rowid int row id of the original row

Example Usage

SELECT
  g.cartodb_id,
  g.the_geom,
  g.the_geom_webmercator,
  (gwr.coeffs->>'pctblack')::numeric as coeff_pctblack,
  (gwr.coeffs->>'pctrural')::numeric as coeff_pctrural,
  (gwr.coeffs->>'pcteld')::numeric as coeff_pcteld,
  (gwr.coeffs->>'pctpov')::numeric as coeff_pctpov,
  gwr.residuals
FROM cdb_crankshaft.CDB_GWR('select * from g_utm'::text, 'pctbach'::text, Array['pctblack', 'pctrural', 'pcteld', 'pctpov']) As gwr
JOIN g_utm as g
on g.cartodb_id = gwr.rowid

Note: See PostgreSQL syntax for parsing JSON objects.

Advanced reading

I think it would be good to have some of the early papers and "the GWR book" as the most base references.

GWR for prediction

GWR in application

Ali, Kamar, Mark D. Partridge, and M. Rose Olfert. 2007. "Can Geographically Weighted Regressions Improve Regional Analysis and Policy Making?" International Regional Science Review 30 (3): 300–329. doi:10.1177/0160017607301609. https://www.researchgate.net/publication/249682503_Can_Geographically_Weighted_Regressions_Improve_Regional_Analysis_and_Policy_Making

Lu, Binbin, Martin Charlton, and A. Stewart Fotheringhama. 2011. "Geographically Weighted Regression Using a Non-Euclidean Distance Metric with a Study on London House Price Data." Procedia Environmental Sciences, Spatial Statistics 2011: Mapping Global Change, 7: 92–97. doi:10.1016/j.proenv.2011.07.017. https://www.researchgate.net/publication/261960122_Geographically_weighted_regression_with_a_non-Euclidean_distance_metric_A_case_study_using_hedonic_house_price_data