tests/exact_mle_given_locs_test.R

In ecrc/exageostatR: R Package Demonstrates the R / C Language Interface for Exageostat

#
# Copyright (c) 2017-2023, King Abdullah University of Science and Technology
# All rights reserved.
#
# ExaGeoStat-R is a software package provided by KAUST
#
#
#
# @file exact_mle_test2.R
# ExaGeoStat R wrapper test example
#
# @version 1.2.0
#
# @author Sameh Abdulah
# @date 2019-01-17
library("exageostatr")                                                  #Load ExaGeoStat-R lib.
sigma_sq        = 1                                                     #Initial variance.
beta            = 0.1                                                   #Initial smoothness.
nu   	        = 0.5                                                   #Initial range.
dmetric         = "euclidean"                                           #0 --> Euclidean distance, 1--> great circle distance.
n               = 1600                                                  #n*n locations grid.
x               = rnorm(n = 1600, mean = 39.74, sd = 25.09)     	#x measurements of n locations.
y               = rnorm(n = 1600, mean = 80.45, sd = 100.19)    	#y measurements of n locations.
theta = c(sigma_sq, beta, nu)
seed = 0
# Initiate exageostat instance
exageostat_init(hardware = list (ncores=2, ngpus=0, ts=320, lts=0,  pgrid=1, qgrid=1))
# Generate Z observation vector based on given locations
data      = simulate_data_exact("ugsm-s", theta, dmetric, n, seed)

# Estimate MLE parameters (Exact)
result    = exact_mle(data, "ugsm-s", dmetric, optimization = 
                      list(clb = c(0.001, 0.001, 0.001), cub = c(5, 5,5 ), tol = 1e-4, max_iters = 3))
# Finalize exageostat instance
exageostat_finalize()