https://github.com/cran/tgp
Tip revision: 5456136e4cca0a4d47d7305f7a4812f3fde04a36 authored by Robert B. Gramacy on 10 January 2006, 18:06:09 UTC
version 1.0-2
version 1.0-2
Tip revision: 5456136
.gdb_history
run
plot(out) # plot the surface
tgp.trees(out) # plot the MAP trees
# try a treed GP
out <- btgp(X=mcycle[,1], Z=mcycle[,2], m0r1=TRUE)
plot(out) # plot the surface
tgp.trees(out) # plot the MAP trees
# try a treed GP LLM
# best to use the "b0" beta linear prior to capture common
# common linear process throughout all regions
out <- btgpllm(X=mcycle[,1], Z=mcycle[,2], bprior="b0",
m0r1=TRUE)
plot(out) # plot the surface
backtrace
:443
up
up
p nn
up
p n
p nn
up
p r
up
p r
down
down
p n
p this
p *this
up
q
q
run
p <- tgp.default.params(2)
p$gamma <- c(0,0.2,0.7) # force no llm
backtrace
q
run
p i
up
run
q