[R-sig-Geo] DATA for GLSM/package geoRglm

[carolina.lang at ifop.cl]

Dear List, 

I am working with binary spatial data collected on hydroacoustic survey, for this i am trying to run a script using GLSM. But this does not work launch me an ERROR, maybe somebody has had the same ERROR, that i do not understand how to fix (last line ). 

SCRIPT 
bin=as.geodata(data.bin,coords.col=1:2,data.col=3,units.m.col=4) # 


model2 <- list(cov.pars = c(1, 0.11), beta = 1, family = "binomial")# 
mcmc.tune <- mcmc.control(S.scale = 0.43, thin = 5,n.iter=20)# 

test.tune <- glsm.mcmc(bin, model= model2, mcmc.input = mcmc.tune) 
#MODELO ESPACIAL 

bin.tune=test.tune 
bin.pre.lf<-prepare.likfit.glsm(bin.tune); str(bin.pre.lf) 
bin.lf<-likfit.glsm(bin.pre.lf,ini.phi=0.1,hessian=TRUE,cov.model="exp")# here ERROR 

results 
-------------------------------------------------------------------- 
likfit.glsm: likelihood maximisation using the function optim. 
phi = 0.1 tausq.rel = 0 
Min. 1st Qu. Median Mean 3rd Qu. Max. 

0 Inf Inf Inf Inf Inf 
Error in .func.val(SivS, SivD, DivD, beta.hat[ll], sigmasq.hat[ll], log.f.c) : Some function values are not finite 







I really hope somebody can help me 




Cheers, CArolina Lang A. 






	
www.ifop.cl 




------------------------------------------------------------------------------------------------
Certificación ISO 9001/2008: Sistema de Datos Bio-Pesqueros ( Arica, Iquique, Coquimbo, Valparaíso, San Antonio, Talcahuano y Calbuco, pesquerías industriales y artesanales)


     


	[[alternative HTML version deleted]]