Hi, before aplied, the function likfit, for one geodata class: or.geo.lf<-likfit(or.geo,cov.model="matern",ini.cov.pars=c(2.5,1),kappa=3,fix.kappa=FALSE,nugget=1.2,lambda=0.019,fix.lambda=FALSE,hessian=TRUE) got this with a summary> summary(or.geo.lf)Summary of the parameter estimation ----------------------------------- Estimation method: maximum likelihood Parameters of the mean component (trend): beta 5.0808 Parameters of the spatial component: correlation function: matern (estimated) variance parameter sigmasq (partial sill) = 0 (estimated) cor. fct. parameter phi (range parameter) = 0 (estimated) extra parameter kappa = 2.422 anisotropy parameters: (fixed) anisotropy angle = 0 ( 0 degrees ) (fixed) anisotropy ratio = 1 Parameter of the error component: (estimated) nugget = 1.859 Transformation parameter: (estimated) Box-Cox parameter = -0.0473 Practical Range with cor=0.05 for asymptotic range: 0.9999339 Maximised Likelihood: log.L n.params AIC BIC "-900" "6" "1812" "1829" non spatial model: log.L n.params AIC BIC "-900" "3" "1806" "1814" Call: likfit(geodata = or.geo, ini.cov.pars = c(2.5, 1), nugget = 1.2, fix.kappa = FALSE, kappa = 3, fix.lambda = FALSE, lambda = 0.019, cov.model = "matern", hessian = TRUE) ############################################################ this two 0s is it ok or not? (estimated) variance parameter sigmasq (partial sill) = 0 (estimated) cor. fct. parameter phi (range parameter) = 0 Thnks -- Freddy Duarte Mu?oz Licenciado en Ciencias Bi?logo Marino Concepci?n Chile (56 - 82454357) (56 - 41- 2738544)