[R-sig-Geo] GWmodel package results

bhartendupandey bhartendu1988pandey at gmail.com
Tue Jul 15 21:07:46 CEST 2014


Dear list, 

I am experimenting with the GWmodel package and I just ran an example code
(with LondonHP data) given in the reference manual for gwr.generalised()
function. Incidently, I ran it in two different computers installed with
different versions of the package. I noted differences between the results
(for both poisson and binomial models) from the two versions:

1. GWmodel_1.2-1
2. GWmodel_1.2-3 

Essentially, I noted that the bandwidth estimated using bw.ggwr() function
from the two versions were different, which is ultimately causing the GWR
model results to vary. Now I wonder what is causing the differences in
bandwidth estimates and more importantly, which version of the package is
reliable? 

I ran the following code:

data(LondonHP)
DM<-gw.dist(dp.locat=coordinates(londonhp))
bw.f1 <- bw.ggwr(BATH2~FLOORSZ,data=londonhp, dMat=DM)
res.poisson<-gwr.generalised(BATH2~FLOORSZ, bw=bw.f1,data=londonhp, dMat=DM)
bw.f2 <- bw.ggwr(BATH2~FLOORSZ,data=londonhp, dMat=DM,family ="binomial")
res.binomial<-gwr.generalised(BATH2~FLOORSZ, bw=bw.f2,data=londonhp,
dMat=DM,
family ="binomial")
 
Here is the sessionInfo() from the two computers: 

Computer 1: 
> sessionInfo()
R version 3.0.1 (2013-05-16)
Platform: x86_64-w64-mingw32/x64 (64-bit)

locale:
[1] LC_COLLATE=English_United States.1252 
[2] LC_CTYPE=English_United States.1252   
[3] LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C                          
[5] LC_TIME=English_United States.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] GWmodel_1.2-3     robustbase_0.91-1 maptools_0.8-30   sp_1.0-15        

loaded via a namespace (and not attached):
[1] DEoptimR_1.0-1  foreign_0.8-53  grid_3.0.1      lattice_0.20-15
[5] tools_3.0.1    

 Computer 2: 
> sessionInfo()
R version 3.0.0 (2013-04-03)
Platform: x86_64-w64-mingw32/x64 (64-bit)

locale:
[1] LC_COLLATE=English_United States.1252 
[2] LC_CTYPE=English_United States.1252   
[3] LC_MONETARY=English_United States.1252
[4] LC_NUMERIC=C                          
[5] LC_TIME=English_United States.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] raster_2.1-16     lattice_0.20-15   GWmodel_1.2-1     robustbase_0.9-10
[5] maptools_0.8-26   sp_1.0-11        

loaded via a namespace (and not attached):
[1] foreign_0.8-53 grid_3.0.0     tools_3.0.0         

loaded via a namespace (and not attached):
[1] DEoptimR_1.0-1  foreign_0.8-53  grid_3.0.1      lattice_0.20-15
tools_3.0.1 

Many thanks for any help and your time.

Best regards,
Bhartendu





--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/GWmodel-package-results-tp7586742.html
Sent from the R-sig-geo mailing list archive at Nabble.com.



More information about the R-sig-Geo mailing list