[R] likfit geodata

Duarte freddyduarte at gmail.com
Fri Nov 7 05:19:21 CET 2008


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)



More information about the R-help mailing list