https://github.com/cran/tgp
Tip revision: f21e634797d9739743606285b5f0d584fac4f24a authored by Robert B. Gramacy on 07 October 2008, 00:00:00 UTC
version 2.1-3
version 2.1-3
Tip revision: f21e634
fried.R
###################################################
### chunk number 1:
###################################################
library(tgp)
##options(width=65)
seed <- 0; set.seed(seed)
###################################################
### 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), pred.n=FALSE, verb=0)
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, pred.n=FALSE, verb=0)
fr.bgpllm.mse <- sqrt(mean((fr.bgpllm$ZZ.mean - ff$Ytrue)^2))
fr.bgpllm.mse
###################################################
### chunk number 5:
###################################################
XX1 <- matrix(rep(0,10), nrow=1)
fr.bgpllm.tr <- bgpllm(X=X, Z=Z, XX=XX1, pred.n=FALSE, trace=TRUE, m0r1=FALSE, verb=0)
###################################################
### chunk number 6:
###################################################
trace <- fr.bgpllm.tr$trace$XX[[1]]
apply(trace[,27:36], 2, mean)
###################################################
### chunk number 7:
###################################################
mean(fr.bgpllm.tr$trace$linarea$ba)
###################################################
### chunk number 8:
###################################################
summary(trace[,9:10])
###################################################
### chunk number 9:
###################################################
apply(trace[,11:15], 2, mean)