https://github.com/cran/tgp
Tip revision: 885fb776f364ad29c761ea5b4d1a8b33d47f1d77 authored by Robert B. Gramacy on 24 December 2005, 01:09:16 UTC
version 1.0
version 1.0
Tip revision: 885fb77
fried.R
###################################################
### chunk number 1:
###################################################
library(tgp)
library(akima)
library(maptree)
###################################################
### chunk number 2:
###################################################
f <- friedman.1.data(200)
ff <- friedman.1.data(1000)
X <- f[,1:10]; Z <- f$Y
XX <- ff[,1:10]
###################################################
### chunk number 3:
###################################################
fr.btlm <- btlm(X=X, Z=Z, XX=XX, tree=c(0.95,2,10), m0r1=TRUE)
fr.btlm.mse <- sqrt(mean((fr.btlm$ZZ.mean - ff$Ytrue)^2))
fr.btlm.mse
###################################################
### chunk number 4:
###################################################
fr.bgpllm <- bgpllm(X=X, Z=Z, XX=XX, m0r1=TRUE)
fr.bgpllm.mse <- sqrt(mean((fr.bgpllm$ZZ.mean - ff$Ytrue)^2))
fr.bgpllm.mse