Semiparametric Geographically Weighted Regression *

* Release 1.0.80 (GWR 4.0.80) *

* 12 March 2014 *
* (Originally coded by T. Nakaya: 1 Nov 2009) *
* *
* Tomoki Nakaya(1), Martin Charlton(2), Paul Lewis(2), *
* Jing Yao (3), A. Stewart Fotheringham (3), Chris Brunsdon (2) *
* (c) GWR4 development team *
* (1) Ritsumeikan University, (2) National University of Ireland, Maynooth, *
* (3) University of St. Andrews *
*****************************************************************************

Program began at 4/2/2024 10:09:16 PM

*****************************************************************************
Session:
Session control file: C:\Users\User\Documents\SEM 4\SPASIAL\kanippeler.ctl
*****************************************************************************
Data filename: C:\Users\User\Documents\SEM 4\SPASIAL\data jatim.csv
Number of areas/points: 38

Model settings---------------------------------
Model type: Poisson
Geographic kernel: adaptive Gaussian
Method for optimal bandwidth search: interval search
Criterion for optimal bandwidth: AIC
Number of varying coefficients: 6
Number of fixed coefficients: 0

Modelling options---------------------------------
Standardisation of independent variables: OFF
Testing geographical variability of local coefficients: On
Local to Global Variable selection: OFF
Global to Local Variable selection: OFF
Prediction at non-regression points: OFF

Variable settings---------------------------------
Area key: field1: Kab/
Easting (x-coord): field12 : longitude
Northing (y-coord): field11: latitude
Cartesian coordinates: Euclidean distance
Dependent variable: field2: Y
Offset variable is not specified
Intercept: varying (Local) intercept
Independent variable with varying (Local) coefficient: field3: X1
Independent variable with varying (Local) coefficient: field5: X3
Independent variable with varying (Local) coefficient: field6: X4
Independent variable with varying (Local) coefficient: field7: X5
Independent variable with varying (Local) coefficient: field8: X6
*****************************************************************************

*****************************************************************************
Global regression result
*****************************************************************************
< Diagnostic information >
Number of parameters: 6
Deviance: 680.595663
Classic AIC: 692.595663
AICc: 695.305341
BIC/MDL: 702.421180
Percent deviance explained 0.547074

Variable Estimate Standard Error z(Est/SE) Exp(Est)

-------------------- --------------- --------------- ---------------

---------------
Intercept 3.489366 0.044527 78.365705
32.765156
X1 0.004184 0.000312 13.397191
1.004193
X3 -0.000016 0.000003 -4.595240
0.999984
X4 0.000027 0.000004 6.604269
1.000027
X5 0.000124 0.000042 2.933706
1.000124
X6 0.000925 0.000544 1.700937
1.000925

*****************************************************************************
GWR (Geographically weighted regression) bandwidth selection
*****************************************************************************

Bandwidth search <interval search> min, max, step

0, 30, 0.5
Bandwdith: 30.000 Dev: 641.295 trace(Hat): 7.875 Criterion: 657.046
Valid_fit
Bandwdith: 29.000 Dev: 638.930 trace(Hat): 7.938 Criterion: 654.807
Valid_fit
Bandwdith: 29.000 Dev: 638.930 trace(Hat): 7.938 Criterion: 654.807
Valid_fit
Bandwdith: 28.000 Dev: 633.390 trace(Hat): 8.075 Criterion: 649.539
Valid_fit
Bandwdith: 28.000 Dev: 633.390 trace(Hat): 8.075 Criterion: 649.539
Valid_fit
Bandwdith: 27.000 Dev: 621.044 trace(Hat): 8.373 Criterion: 637.790
Valid_fit
Bandwdith: 27.000 Dev: 621.044 trace(Hat): 8.373 Criterion: 637.790
Valid_fit
Bandwdith: 26.000 Dev: 616.497 trace(Hat): 8.543 Criterion: 633.582
Valid_fit
Bandwdith: 26.000 Dev: 616.497 trace(Hat): 8.543 Criterion: 633.582
Valid_fit
Bandwdith: 25.000 Dev: 615.468 trace(Hat): 8.627 Criterion: 632.722
Valid_fit
Bandwdith: 25.000 Dev: 615.468 trace(Hat): 8.627 Criterion: 632.722
Valid_fit
Bandwdith: 24.000 Dev: 602.415 trace(Hat): 8.900 Criterion: 620.214
Valid_fit
Bandwdith: 24.000 Dev: 602.415 trace(Hat): 8.900 Criterion: 620.214
Valid_fit
Bandwdith: 23.000 Dev: 595.837 trace(Hat): 9.047 Criterion: 613.931
Valid_fit
Bandwdith: 23.000 Dev: 595.837 trace(Hat): 9.047 Criterion: 613.931
Valid_fit
Bandwdith: 22.000 Dev: 587.096 trace(Hat): 9.199 Criterion: 605.495
Valid_fit
Bandwdith: 22.000 Dev: 587.096 trace(Hat): 9.199 Criterion: 605.495
Valid_fit
Bandwdith: 21.000 Dev: 579.398 trace(Hat): 9.401 Criterion: 598.199
Valid_fit
Bandwdith: 21.000 Dev: 579.398 trace(Hat): 9.401 Criterion: 598.199
Valid_fit
Bandwdith: 20.000 Dev: 573.433 trace(Hat): 9.589 Criterion: 592.612
Valid_fit
Bandwdith: 20.000 Dev: 573.433 trace(Hat): 9.589 Criterion: 592.612
Valid_fit
Bandwdith: 19.000 Dev: 555.328 trace(Hat): 9.885 Criterion: 575.098
Valid_fit
Bandwdith: 19.000 Dev: 555.328 trace(Hat): 9.885 Criterion: 575.098
Valid_fit
Bandwdith: 18.000 Dev: 544.480 trace(Hat): 10.092 Criterion: 564.665
Valid_fit
Bandwdith: 18.000 Dev: 544.480 trace(Hat): 10.092 Criterion: 564.665
Valid_fit
Bandwdith: 17.000 Dev: 525.234 trace(Hat): 10.390 Criterion: 546.013
Valid_fit
Bandwdith: 17.000 Dev: 525.234 trace(Hat): 10.390 Criterion: 546.013
Valid_fit
Bandwdith: 16.000 Dev: 515.443 trace(Hat): 10.744 Criterion: 536.931
Valid_fit
Bandwdith: 16.000 Dev: 515.443 trace(Hat): 10.744 Criterion: 536.931
Valid_fit
Bandwdith: 15.000 Dev: 504.180 trace(Hat): 11.018 Criterion: 526.216
Valid_fit
Bandwdith: 15.000 Dev: 504.180 trace(Hat): 11.018 Criterion: 526.216
Valid_fit
Bandwdith: 14.000 Dev: 483.014 trace(Hat): 11.417 Criterion: 505.849
Valid_fit
Bandwdith: 14.000 Dev: 483.014 trace(Hat): 11.417 Criterion: 505.849
Valid_fit
Bandwdith: 13.000 Dev: 470.596 trace(Hat): 11.802 Criterion: 494.201
Valid_fit
Bandwdith: 13.000 Dev: 470.596 trace(Hat): 11.802 Criterion: 494.201
Valid_fit
Bandwdith: 12.000 Dev: 449.530 trace(Hat): 12.161 Criterion: 473.852
Valid_fit
Bandwdith: 12.000 Dev: 449.530 trace(Hat): 12.161 Criterion: 473.852
Valid_fit
Bandwdith: 11.000 Dev: 443.203 trace(Hat): 12.779 Criterion: 468.762
Valid_fit
Bandwdith: 11.000 Dev: 443.203 trace(Hat): 12.779 Criterion: 468.762
Valid_fit
Bandwdith: 10.000 Dev: 432.130 trace(Hat): 13.217 Criterion: 458.564
Valid_fit
Bandwdith: 10.000 Dev: 432.130 trace(Hat): 13.217 Criterion: 458.564
Valid_fit
Bandwdith: 9.000 Dev: 417.483 trace(Hat): 13.763 Criterion: 445.009
Valid_fit
Bandwdith: 9.000 Dev: 417.483 trace(Hat): 13.763 Criterion: 445.009
Valid_fit
Bandwdith: 8.000 Dev: 396.940 trace(Hat): 14.640 Criterion: 426.221
Valid_fit
Bandwdith: 8.000 Dev: 396.940 trace(Hat): 14.640 Criterion: 426.221
Valid_fit
Bandwdith: 7.000 Dev: 363.160 trace(Hat): 15.546 Criterion: 394.252
Valid_fit
 Bandwdith: 7.000 Dev: 363.160 trace(Hat): 15.546 Criterion: 394.252
Valid_fit
Bandwdith: 6.000 Dev: 335.584 trace(Hat): 16.969 Criterion: 369.522
Valid_fit
Bandwdith: 6.000 Dev: 335.584 trace(Hat): 16.969 Criterion: 369.522
Valid_fit
Bandwdith: 5.000 Dev: 315.547 trace(Hat): 19.218 Criterion: 353.983
Valid_fit
Bandwdith: 5.000 Dev: 315.547 trace(Hat): 19.218 Criterion: 353.983
Valid_fit
Bandwdith: 4.000 Dev: 230.857 trace(Hat): 21.531 Criterion: 273.920
Valid_fit
Bandwdith: 4.000 Dev: 230.857 trace(Hat): 21.531 Criterion: 273.920
Valid_fit
Bandwdith: 3.000 Dev: 99.210 trace(Hat): 28.152 Criterion: 155.513
Valid_fit
Bandwdith: 3.000 Dev: 99.210 trace(Hat): 28.152 Criterion: 155.513
Valid_fit
Bandwdith: 2.000 Dev: NaN trace(Hat): NaN Criterion: NaN
Invalid_fit
Bandwdith: 2.000 Dev: NaN trace(Hat): NaN Criterion: NaN
Invalid_fit
Bandwdith: 1.000 Dev: NaN trace(Hat): NaN Criterion: NaN
Invalid_fit
Bandwdith: 1.000 Dev: NaN trace(Hat): NaN Criterion: NaN
Invalid_fit

Error in the initial weight calculation loop

Index was outside the bounds of the array. Bandwdith: 0.000 Dev: NaN
trace(Hat): NaN Criterion: NaN Invalid_fit

Error in the initial weight calculation loop

Index was outside the bounds of the array. Bandwdith: 0.000 Dev: NaN
trace(Hat): NaN Criterion: NaN Invalid_fit
Best bandwidth size 3.500
Minimum AIC 155.513

*****************************************************************************
GWR (Geographically weighted regression) result
*****************************************************************************
Bandwidth and geographic ranges
Bandwidth size: 3.000000
Coordinate Min Max Range
--------------- --------------- --------------- ---------------
X-coord 108.800194 124.466957 15.666763
Y-coord -8.749315 0.882681 9.631996

Diagnostic information
Effective number of parameters (model: trace(S)): 28.151665
Effective number of parameters (variance: trace(S'WSW^-1)): 24.900117
Degree of freedom (model: n - trace(S)): 9.848335
Degree of freedom (residual: n - 2trace(S) + trace(S'WSW^-1)): 6.596786
Deviance: 99.210081
Classic AIC: 155.513411
AICc: 341.009998
BIC/MDL: 201.614189
Percent deviance explained 0.933977

***********************************************************
 << Geographically varying (Local) coefficients >>
***********************************************************
Estimates of varying coefficients have been saved in the following file.
Listwise output file: C:\Users\User\Documents\SEM 4\SPASIAL\
kanippeler_listwise.csv

Summary statistics for varying (Local) coefficients

Variable Mean STD
-------------------- --------------- ---------------
Intercept 2.756183 0.810792
X1 -0.000648 0.011742
X3 -0.000045 0.000097
X4 0.000022 0.000072
X5 0.001613 0.002088
X6 -0.000673 0.012040

Variable Min Max Range

-------------------- --------------- --------------- ---------------
Intercept 0.588017 4.116024 3.528007
X1 -0.048152 0.011198 0.059350
X3 -0.000493 0.000037 0.000530
X4 -0.000093 0.000319 0.000412
X5 -0.000897 0.009131 0.010028
X6 -0.026305 0.043889 0.070194

Variable Lwr Quartile Median Upr Quartile

-------------------- --------------- --------------- ---------------
Intercept 2.144282 3.062276 3.463754
X1 -0.005330 0.001256 0.006244
X3 -0.000057 -0.000016 -0.000001
X4 -0.000007 0.000009 0.000029
X5 -0.000072 0.001410 0.002788
X6 -0.004277 0.001003 0.004889

Variable Interquartile R Robust STD

-------------------- --------------- ---------------
Intercept 1.319471 0.978111
X1 0.011574 0.008580
X3 0.000055 0.000041
X4 0.000036 0.000027
X5 0.002860 0.002120
X6 0.009165 0.006794
(Note: Robust STD is given by (interquartile range / 1.349) )

*****************************************************************************
GWR Analysis of Deviance Table
*****************************************************************************
Source Deviance DOF Deviance/DOF
------------ ------------------- ---------- ----------------
Global model 680.596 32.000 21.269
GWR model 99.210 6.597 15.039
Difference 581.386 25.403 22.886

*************************************************************************
Geographical variability tests of local coefficients
*************************************************************************
Variable Diff of deviance Diff of DOF DIFF of Criterion
-------------------- ------------------ ---------------- -----------------
Intercept 44.822383 2.729285 -39.363812
X1 31.393805 1.714230 -27.965346
X3 32.957049 1.522827 -29.911395

< Warning Backfitting iteration reached MAXITER >

X4 29.866262 2.212602 -25.441058
X5 73.943221 2.121505 -69.700212
X6 NaN NaN NaN
-------------------- ------------------ ---------------- -----------------
Note: positive value of diff-Criterion (AICc, AIC, BIC/MDL or CV) suggests no
spatial variability in terms of model selection criteria.
Chi-square test: in case of no spatial variability, [Diff of deviance] follows
the Chi-square distribution (DOF is the diff of DOF).

**********************************************