require(geoRglm) SP3107 = read.table("SP3107semponto.TXT", head=F) SP3107Geo<-as.geodata(SP3107, head=F, coords=1:2, data.col=4, units.m.col=3) bor <- SP3107Geo$coords[chull(SP3107Geo$coords),] SP3107Geo$borders <- bor plot(SP3107Geo) mcmc.5 <- mcmc.control(S.scale = 0.6, thin=20,n.iter=50000, burn.in=1000, S.start="random") model.5 <- list(cov.pars=c(0.6, 0.01), beta=1, family="poisson") outmcmc.5 <- glsm.mcmc(SP3107Geo, units.m=SP3107Geo$units.m, model= model.5, mcmc.input= mcmc.5) mcmcobj.5 <- prepare.likfit.glsm(outmcmc.5) lik.5 <- likfit.glsm(mcmcobj.5, ini.phi = 1, fix.nugget.rel = FALSE) print(lik.5) summary(lik.5) lik.5.mat.nugget <- likfit.glsm(mcmcobj.5, ini.phi =1,cov.model = "matern", nugget.rel = 0.385) print(lik.5.mat.nugget) summary(lik.5.mat.nugget) lik.5.mat.nugget