Skip to content

logistic GWR

1 message · Seth-3

#
Yes, it is very possible.  

memory.limit(size=2800)
library(sp)
library(spgwr)
spatdata<-read.table(file="C:\\ctlogreg_rfordata\\data\\650GRD_AOIM1KNONDEV85DRYLNDCNSH2OMSK_85to06chcommameancent.txt",header=TRUE,sep=",")


ggwr5k<-ggwr(devchng~
  agricden          +
grassden          +
agric             +
barren            +
grass             +
wtlnd             +
allcitydis.w      +
devden            +
dis.citge100k.wab +
dis.prim.wab      +
dis.secloc.wab    +
dist.rampm.wab    +
elev.wab          +
fema100yr.w01     +
medinc.wab        +
munifeddep.w01    +
popinvdis.wab     +
sewer01.w         +
slope.deg.wab     
,
data=spatdata,coords=cbind(spatdata$eastdis.mab,spatdata$northdist.mab),
bandwidth=5000,gweight=gwr.Gauss,family=binomial(link="logit"),longlat=FALSE,type=c("response"))