Kriging::save & Kriging::load

Description

Save/Load a Kriging Model

Usage

  • Python

    # k = Kriging(...)
    k.save("k.json")
    k2 = load("k.json")
    
  • R

    # k = Kriging(...)
    k$save("k.json")
    k2 = load("k.json")
    
  • Matlab/Octave

    % k = Kriging(...)
    k.save("k.json")
    k2 = load("k.json")
    
  • Julia

    # k = Kriging(...)
    save(k, "k.json")
    k2 = load("k.json")
    

Arguments

Argument

Description

filename

Path to the JSON file used for persistence.

Details

save() serialises the fitted model to JSON, and load() reconstructs a compatible Kriging object from that file. This is useful for checkpointing fitted models or exchanging them across wrapper sessions.

Value

save() writes the file and returns invisibly; load() returns the restored object.

Examples

f <- function(x) 1 - 1 / 2 * (sin(12 * x) / (1 + x) + 2 * cos(7 * x) * x^5 + 0.7)
set.seed(123)
X <- as.matrix(runif(10))
y <- f(X)

k <- Kriging(y, X, kernel = "matern3_2")
k

k$save("k.json")
print(load("k.json"))

Results

* data: 10x[0.0455565,0.940467] -> 10x[0.194057,1.00912]
* trend constant (est.): 0.433954
* variance (est.): 0.0873685
* covariance:
  * kernel: matern3_2
  * range (est.): 0.240585
  * fit:
    * objective: LL
    * optim: BFGS
* data: 10x[0.0455565,0.940467] -> 10x[0.194057,1.00912]
* trend constant (est.): 0.433954
* variance (est.): 0.0873685
* covariance:
  * kernel: matern3_2
  * range (est.): 0.240585
  * fit:
    * objective: LL
    * optim: BFGS
{
    "F": {
        "base64_data": "AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8=",
        "n_cols": 1,
        "n_rows": 10,
        "type": "mat"
    },
    "M": {
        "base64_data": "AAAAAAAA8D+DnnzBJDfsPyKj0qxEzMI/VExeoaYW1z8NZ44RK9jWP4YTOPTDYOU/rGAwjOUdyT+tPdoJrWuwvwylmewT6a8/otYuWrdSjT8=",
        "n_cols": 1,
        "n_rows": 10,
        "type": "mat"
    },
    "R": {
        "base64_data": "AAAAAAAA8D+qznH4ygbAP6YuzrEABuk/dGW40r+bsj8MkJdIFYmqP4s0WkGIu94//hGZxWLx3j9q3YLs+pyxP4bBlw5pxNs/aK0oahwC5T+qznH4ygbAPwAAAAAAAPA/XsyYHRYezz8W9aReFDfrP6hlLXsMbOY/7InDlNrxnj+E1mEuzz/cP/pdq0Qmdeo/VZhNMVt23z81BvRDoejTP6YuzrEABuk/XsyYHRYezz8AAAAAAADwP0WyeeRvnMI/uvgy47Lruj8av2b7P+nQPz/CXSXJNuk/t6pQ++uowT8iRDf+Lj7nP8bpQnJqf+4/dGW40r+bsj8W9aReFDfrP0WyeeRvnMI/AAAAAAAA8D8bpKA1l+ntP1MjsT4zVJE/HcM8y8Ss0T8bxOt7D+7vPywcZYhs69M/pBo3WtkvyD8MkJdIFYmqP6hlLXsMbOY/uvgy47Lruj8bpKA1l+ntPwAAAAAAAPA/FaTpVzlEiD9MGWRpVBjKP/nNd/+Bee4/JMGayaaPzT9I7GJIjp7BP4s0WkGIu94/7InDlNrxnj8av2b7P+nQP1MjsT4zVJE/FaTpVzlEiD8AAAAAAADwPxPv6xKOv8E/1muwm8pZkD/eR7A9iyK/P8GDQldfR8o//hGZxWLx3j+E1mEuzz/cPz/CXSXJNuk/HcM8y8Ss0T9MGWRpVBjKPxPv6xKOv8E/AAAAAAAA8D8/KB2AJtXQPy7BK52cmO8/3zFjq0z47D9q3YLs+pyxP/pdq0Qmdeo/t6pQ++uowT8bxOt7D+7vP/nNd/+Bee4/1muwm8pZkD8/KB2AJtXQPwAAAAAAAPA/u0B5tuv80j8a659b2/rGP4bBlw5pxNs/VZhNMVt23z8iRDf+Lj7nPywcZYhs69M/JMGayaaPzT/eR7A9iyK/Py7BK52cmO8/u0B5tuv80j8AAAAAAADwPx0+ZfvhNOs/aK0oahwC5T81BvRDoejTP8bpQnJqf+4/pBo3WtkvyD9I7GJIjp7BP8GDQldfR8o/3zFjq0z47D8a659b2/rGPx0+ZfvhNOs/AAAAAAAA8D8=",
        "n_cols": 10,
        "n_rows": 10,
        "type": "mat"
    },
    "Rinv": {
        "base64_data": "iWZx8Ch3EECacYl1FOp9v6IBZw272xrAOiPUifltjD9eKRkhHf4nPyYnreNRg/C/7YQsVRxh878E9LTdSTWGv5hXT43WQeI/LRgox1oWEkAAaIl1FOp9vzyBb5LZvx5Az9byPx/DvD84UwpovyxEwOQ9Lii77uK/rsgSBZ5ERT+k65T/HjoVQDQ0g0dAhkFAnD1vmyzyGcBe7eIdQKfXv6QBZw272xrA59byPx/DvD/y/5Qgo/08QMl9WjGoVcu/cNf6c4sRZ7+IgKvIEAXpPyTsRm4HojJAIejmhFBaxT9dxpEt0o0hwIZObf4CaEDAAC7UifltjD9BUwpovyxEwBl9WjGoVcu/38jJKuGoiEDMxPnsNPpXQDCFcFNgNlS/p5xqrGYsJMC74xZyuHSKwCK9IlUJhylAJV8nc7N65j8AaBkhHf4nP1k9Lii77uK/BOT6c4sRZ7/MxPnsNPpXQKks74j6LTpATDt2/tEO8b5jiLTiZgbBv/FRkBAgCl7Alu6W5B+LxT+1howgk/iCPyYnreNRg/C/J8kSBZ5ERT+KgKvIEAXpPwGGcFNgNlS/Rj12/tEO8b4Bqq/nts31P6/8mfpWjrs/iZj3dBGUTz9b88fl2/Wpv0uDYlIHuNm/FIUsVRxh87+l65T/HjoVQC3sRm4HojJAq5xqrGYsJMB2iLTiZgbBv638mfpWjrs/ltVaiUMWYEDAWwjRe4QfQOfUb22Kt1jAfkmapPRJScAADrXdSTWGv1Q0g0dAhkFAN+XmhFBaxT+74xZyuHSKwPFRkBAgCl7A8pf3dBGUTz/AWwjRe4QfQLd4H9ozIY1AkSogpu7wI8Bv7g/TXI/hv5xXT43WQeI/qz1vmyzyGcBmxpEt0o0hwCe9IlUJhylAn+6W5B+LxT9y88fl2/Wpv+bUb22Kt1jAkSogpu7wI8Db3KqZ5JVUQHBW3zc33zxAJRgox1oWEkCM7eIdQKfXv4ZObf4CaEDAPF8nc7N65j+qh4wgk/iCP0eDYlIHuNm/fkmapPRJScBp7g/TXI/hv3BW3zc33zxAU5vUQH9aSUA=",
        "n_cols": 10,
        "n_rows": 10,
        "type": "mat"
    },
    "T": {
        "base64_data": "AAAAAAAA8D+qznH4ygbAP6YuzrEABuk/dGW40r+bsj8MkJdIFYmqP4s0WkGIu94//hGZxWLx3j9q3YLs+pyxP4bBlw5pxNs/aK0oahwC5T8AAAAAAAAAAPYyL62Iv+8/2CHFy4K7wj+7kEpZMSPrP/dD6ET/Y+Y/VXoAXnzdnr9fAsIKfZHYP+5KolO+Y+o/KisFUig13D+tn94+O4fNPwAAAAAAAAAAAAAAAAAAAABc19WHVWPjP9Oplr4XC66/AiUZ9sjqr7+AAZ+NPpPGv0Hi1kYdruI/M+p6qzzdrr/9awKuYgnhPwh6HUJZcOU/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAALuXn5yOx4D+ujGPTtnzkP+FyplqTv3W/r2KaPJmruL/BMbKy9sjhP/pywjoq1b6/i1fTYUacnb8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAANlXCLPzndM/ZAC3eEttWD8GcMZPcrWbP2QYCa7eTKQ/WomOxKdQoT9KKxUt8aCAPwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAC2l9fyeXrrP0/B9qb0Z5g/RLvllZ1wDD9ktlNFtJSXP6WllJPs4pE/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAr/cyryCN4D87/X63Ir1aP5Dhy9/6muE/k+E/4gXxyT8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAALyJRjJWyKA/fURpxw+tcz/YoLBnKYBjvwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAys2mO34i/PxKMUdKX9LG/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA5b4Sy6z5wT8=",
        "n_cols": 10,
        "n_rows": 10,
        "type": "mat"
    },
    "X": {
        "base64_data": "AAAAi6tn0j8AAACxyznpPwAAAIqtLNo/AAAAt61B7D8AAKDYThjuPwAAAHYuU6c/AAAAez3m4D8AAEBiso7sPwAAIAtbpeE/AADAAi053T8=",
        "n_cols": 1,
        "n_rows": 10,
        "type": "mat"
    },
    "beta": {
        "base64_data": "6LhstujF2z8=",
        "size": 1,
        "type": "colvec"
    },
    "centerX": {
        "base64_data": "AAAAAAAAAAA=",
        "size": 1,
        "type": "rowvec"
    },
    "centerY": 0.0,
    "circ": {
        "base64_data": "Cu4VN0KM+T8=",
        "n_cols": 1,
        "n_rows": 1,
        "type": "mat"
    },
    "content": "Kriging",
    "covType": "matern3_2",
    "dX": {
        "base64_data": "AAAAAAAAAAAAAIDr9QXgvwAAAPwHFL+/AACA8dcN478AACATeeTkvwAAgHiL+s4/AAAA1p7Jzr8AAMCc3FrjvwAAQIsK49C/AACA7wKjxb8AAIDr9QXgvwAAAAAAAAAAAAAA2OlG2D8AAAAwED+4vwAAgJ4MesO/AACgyZjE5z8AAABsHKfQPwAAAIo1p7q/AACAl8JRzj8AAEBfajrVPwAAAPwHFL+/AAAA2OlG2D8AAAAAAAAAAAAAAOStVt6/AACgE/gB4b8AAEC7R0LXPwAAALA1f76/AACAOrfw3r8AAIAYETzCvwAAAMb7Y6i/AACA8dcN478AAAAwED+4vwAAAOStVt6/AAAAAAAAAAAAAAAaEmqtvwAAoM96zOo/AAAAeOC21j8AAADQKkGDvwAAwFelONU/AABAay5K2z8AACATeeTkvwAAgJ4MesO/AACgE/gB4b8AAAAaEmqtvwAAAAAAAAAAAABA8Ruj7D8AAEC7ImTaPwAAAGbHmag/AAAAm+fl2D8AAICucPfePwAAgHiL+s4/AACgyZjE5z8AAEC7R0LXPwAAoM96zOo/AABA8Ruj7D8AAAAAAAAAAAAAQCcV4t6/AADgen8Z678AAMAjKDDgvwAAADTHTtq/AAAA1p7Jzr8AAABsHKfQPwAAALA1f76/AAAAeOC21j8AAEC7ImTaPwAAQCcV4t6/AAAAAAAAAAAAAIDO6VDXvwAAAASy45e/AAAAzTdNsj8AAMCc3FrjvwAAAIo1p7q/AACAOrfw3r8AAADQKkGDvwAAAGbHmag/AADgen8Z678AAIDO6VDXvwAAAAAAAAAAAABArq7S1T8AAMDBN+TbPwAAQIsK49C/AACAl8JRzj8AAIAYETzCvwAAwFelONU/AAAAm+fl2D8AAMAjKDDgvwAAAASy45e/AABArq7S1T8AAAAAAAAAAAAAAE4kRrg/AACA7wKjxb8AAEBfajrVPwAAAMb7Y6i/AABAay5K2z8AAICucPfePwAAADTHTtq/AAAAzTdNsj8AAMDBN+TbPwAAAE4kRrg/AAAAAAAAAAA=",
        "n_cols": 100,
        "n_rows": 1,
        "type": "mat"
    },
    "est_beta": true,
    "est_sigma2": true,
    "est_theta": true,
    "maxdX": {
        "base64_data": "AABA8Ruj7D8=",
        "size": 1,
        "type": "colvec"
    },
    "noise_model": "None",
    "normalize": false,
    "objective": "LL",
    "optim": "BFGS",
    "regmodel": "constant",
    "scaleX": {
        "base64_data": "AAAAAAAA8D8=",
        "size": 1,
        "type": "rowvec"
    },
    "scaleY": 1.0,
    "sigma2": 0.08736845168977267,
    "star": {
        "base64_data": "",
        "n_cols": 0,
        "n_rows": 0,
        "type": "mat"
    },
    "theta": {
        "base64_data": "yoo1OXvLzj8=",
        "size": 1,
        "type": "colvec"
    },
    "version": 2,
    "y": {
        "base64_data": "Ett2oSGc6D9CynZ5kSvcP1bGVfNZJfA/VPt7kN8T1z9gesnC3NbIPyjBELSFsNk/hJLAxOFZ5T/2Fal4Gp7VP1604hmqpeI/QFngxtBX7T8=",
        "size": 10,
        "type": "colvec"
    },
    "z": {
        "base64_data": "PP2AjFpy1T9qRm7pWnOiv/PExZX10OA/vgBKp0DUsb+ZgKvemybgv/SROyTYQb6/vqAWle442r/fews4XEuoP69aDk0p3Mi/RcxzjXrBqb8=",
        "size": 10,
        "type": "colvec"
    }
}