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: ker...@sccwrp.org
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
