[R-sig-Geo] geoR- "likfit" not estimating range parameters

[Kerry Ritter]

Hi. I am using the library "geoR" and am having trouble getting REML 
estimates for the range and  different from those specified in the 
intial values.

My statement is....

reml.fit=likfit(SDmatchC.geo, ini=c(100,2020), fix.nugget = FALSE, 
nugget = 5,  fix.psiA = FALSE, psiA = 90, fix.psiR = FALSE, psiR = 1.4,
       cov.model = "gaussian", lik.method = "REML")

The results is
 > reml.fit
likfit: estimated model parameters:
      beta      tausq    sigmasq        phi       psiA       psiR
"  41.513" "  12.083" " 267.910" "2019.996" "   5.806" "   5.998"
Practical Range with cor=0.05 for asymptotic range: 3496.247

Note that if I replace 2020 with 3020, I get "3019.99" for phi.  So it 
is not really estimating the
range parameters.

What am I missing?

Thanks,
Kerry


-- 


**********************
Kerry Ritter, Ph.D.
statistician
Southern California Coastal Water Research Project
3535 Harbor Blvd., Suite 110

work: 714-755-3210
cell: 714-420-3346
fax:  714-755-3299

email: kerryr at sccwrp.org