https://github.com/cran/tgp
Revision 80e1597735d81e97c796d89aefee8213d9c1b44d authored by Robert B. Gramacy on 17 April 2008, 18:06:09 UTC, committed by cran-robot on 17 April 2008, 18:06:09 UTC
1 parent c3e1134
Tip revision: 80e1597735d81e97c796d89aefee8213d9c1b44d authored by Robert B. Gramacy on 17 April 2008, 18:06:09 UTC
version 2.1-2
version 2.1-2
Tip revision: 80e1597
pred.R
###################################################
### chunk number 1:
###################################################
library(tgp)
library(maptree)
##options(width=65)
seed <- 0; set.seed(seed)
###################################################
### chunk number 2:
###################################################
library(MASS)
out <- btgpllm(X=mcycle[,1], Z=mcycle[,2], bprior="b0",
m0r1=TRUE, pred.n=FALSE, verb=0)
save(out, file="out.Rsave")
out <- NULL
###################################################
### chunk number 3:
###################################################
load("out.Rsave")
XX <- seq(2.4, 56.7, length=200)
out.kp <- predict(out, XX=XX, pred.n=FALSE)
###################################################
### chunk number 4:
###################################################
out.p <- predict(out, XX=XX, pred.n=FALSE, BTE=c(0,1000,1))
###################################################
### chunk number 5:
###################################################
out2 <- predict(out, XX, pred.n=FALSE, BTE=c(0,2000,2),
krige=FALSE)
###################################################
### chunk number 6: kp
###################################################
plot(out.kp, center="km", as="ks2")
###################################################
### chunk number 7:
###################################################
rl <- readline("press RETURN to continue: ")
graphics.off()
###################################################
### chunk number 8: p
###################################################
plot(out.p)
###################################################
### chunk number 9:
###################################################
rl <- readline("press RETURN to continue: ")
graphics.off()
###################################################
### chunk number 10: 2
###################################################
plot(out2)
###################################################
### chunk number 11:
###################################################
rl <- readline("press RETURN to continue: ")
graphics.off()
###################################################
### chunk number 12:
###################################################
unlink("out.Rsave")
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...