[R-sig-Geo] How do I have to work with Poisson distribution with geoRglm (count data)?

Jimmy Neutron jimmyjmv at hotmail.com
Mon Mar 30 23:13:20 CEST 2015


Dear comRades:
Because of my data, I have just realized that I have to work with Poisson, because I have 'count data' with geographic reference. Then, my my 'geodata' is such as:
       xUTMkm   yUTMkm 'Count'1385 450.1202 1011.425   01386 450.4273 1007.219   01387 450.2584 1011.884   01388 450.1696 1010.261   01389 450.1718 1009.887   01390 450.6981 1004.379   0...
I read the geoRglm structure. What does it mean that I have to make a empiric variogram?.
When I worked with Binomial model (True, False), my script was as following:a2008.posCEROS.spmod<-list(cov.pars=c(1,20),beta=1.0,cov.model="matern",nugget=0,kappa=0.35,family="binomial",link="logit")
Then, I suppose that my Poisson model would be as following:a2008.posCEROS.spmod<-list(cov.pars=c(1,20),beta=1.0,cov.model="matern",nugget=0,kappa=0.35,family="poisson",link="logit")
Do I have to calculate nugget and kappa by some method (like likelihood) from my 'count data' or it is merely a theoretical model?. Is my Poisson model right?.
After that, I will generate MCMC simulations, such as:
a2008.posCEROS.mcmc<-mcmc.control(S.scale=0.582, thin=10) #mcmc marc of change monte carlo, S.scalea2008.posCEROS.tune<-glsm.mcmc(a2008.posCEROSbin, model=a2008.posCEROS.spmod, mcmc.input=a2008.posCEROS.mcmc)
My goal is to predict how many 'Count' do I have in the survey.
Thanks in advance for youR help. 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150330/8b134a4d/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 2008TrialsSuccesses.gif
Type: image/gif
Size: 15953 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150330/8b134a4d/attachment.gif>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 2008ploteo.gif
Type: image/gif
Size: 11859 bytes
Desc: not available
URL: <https://stat.ethz.ch/pipermail/r-sig-geo/attachments/20150330/8b134a4d/attachment-0001.gif>


More information about the R-sig-Geo mailing list