[R-sig-Geo] logistic GWR

[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"))
-- 
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/logistic-GWR-tp4115587p5236339.html
Sent from the R-sig-geo mailing list archive at Nabble.com.