https://github.com/cran/fields
Tip revision: 8ad4e6ec8732718f8759ef4472695e9bfb5ead73 authored by Douglas Nychka on 23 April 2019, 11:50:03 UTC
version 9.7
version 9.7
Tip revision: 8ad4e6e
mKrigMLETest.Rout.save
R Under development (unstable) (2019-04-06 r76327) -- "Unsuffered Consequences"
Copyright (C) 2019 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin17.7.0 (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> # fields is a package for analysis of spatial data written for
> # the R software environment .
> # Copyright (C) 2018
> # University Corporation for Atmospheric Research (UCAR)
> # Contact: Douglas Nychka, nychka@ucar.edu,
> # National Center for Atmospheric Research,
> # PO Box 3000, Boulder, CO 80307-3000
> #
> # This program is free software; you can redistribute it and/or modify
> # it under the terms of the GNU General Public License as published by
> # the Free Software Foundation; either version 2 of the License, or
> # (at your option) any later version.
> # This program is distributed in the hope that it will be useful,
> # but WITHOUT ANY WARRANTY; without even the implied warranty of
> # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
> # GNU General Public License for more details.
>
> suppressMessages(library( fields ))
> options( echo=FALSE)
PASSED test at tolerance 1e-08
PASSED test at tolerance 1e-08
PASSED test at tolerance 1e-08
PASSED test at tolerance 1e-08
PASSED test at tolerance 1e-08
PASSED test at tolerance 0.002
PASSED test at tolerance 1e-08
PASSED test at tolerance 1e-08
PASSED test at tolerance 0.02
PASSED test at tolerance 0.02
PASSED test at tolerance 0.002
PASSED test at tolerance 0.02
PASSED test at tolerance 0.02
PASSED test at tolerance 0.002
PASSED test at tolerance 0.02
PASSED test at tolerance 0.02
PASSED test at tolerance 0.002
lnProfileLike.FULL lambda theta sigmaMLE
-61.8145417 0.1015556 0.1002546 0.3185669
rhoMLE funEval gradEval
0.9993032 8.0000000 6.0000000
lnProfileREML.FULL lambda theta sigmaMLE
-61.8145417 0.1015556 0.1002546 0.3185669
rhoMLE funEval gradEval
0.9993032 8.0000000 6.0000000
lnProfileLike.FULL lambda theta sigmaMLE
-61.8143732 0.1015557 0.1002514 0.3185628
rhoMLE funEval gradEval
0.9992768 8.0000000 6.0000000
1 2 3 4 5 6 7 8 9 10 all done with mKrigMLEGrid tests
>
> test.for.zero( trueTheta, mean(hold[,3]), tol=2e-3,tag="Monte Carlo theta")
Testing: Monte Carlo theta
PASSED test at tolerance 0.002
> test.for.zero( trueLambda, mean(hold[,2]), tol=2e-2,tag="Monte Carlo theta")
Testing: Monte Carlo theta
PASSED test at tolerance 0.02
>
>
> proc.time()
user system elapsed
10.217 1.767 10.274