https://github.com/cran/tgp
Tip revision: 799f0308754b1689d24d3e6a886fa5d74144c627 authored by Robert B. Gramacy on 25 November 2008, 00:00:00 UTC
version 2.1-6
version 2.1-6
Tip revision: 799f030
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)