
R version 3.4.0 (2017-04-21) -- "You Stupid Darkness"
Copyright (C) 2017 The R Foundation for Statistical Computing
Platform: x86_64-apple-darwin15.6.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) 2017
>   # 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.
> # Test adapted from fields package, under GPL license
> 
> library( fields )
Loading required package: spam
Loading required package: grid
Spam version 1.4-0 (2016-08-29) is loaded.
Type 'help( Spam)' or 'demo( spam)' for a short introduction 
and overview of this package.
Help for individual functions is also obtained by adding the
suffix '.spam' to the function name, e.g. 'help( chol.spam)'.

Attaching package: 'spam'

The following objects are masked from 'package:base':

    backsolve, forwardsolve

Loading required package: maps
> 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 
 17.036   1.064  18.124 
