# (C) Petteri Pajunen 2007 library(tgp) X11(width=5,height=5) nsamples <- 40 standarddev <- 0.5 X <- seq(-2,2, length=nsamples) XX <- seq(-2,2,length=49) Zlin <- 1+2*X+rnorm(length(X),sd=standarddev) Znlin <- 1+2*X+3*X^2+rnorm(length(X),sd=standarddev) # mixture model with linear data lin.gpllm.tr <- bgpllm(X=X,Z=Zlin,XX=XX,bprior="b0",trace=TRUE, verb=0) plot(lin.gpllm.tr, main="FIT MIXTURE MODEL TO 1+2x", layout="surf") pb <- mean(lin.gpllm.tr$trace$linarea$la) cat("POSTERIOR PROBABILITY OF LINEAR MODEL ",100*pb," PERCENT\n") scan(quiet=TRUE) # mixture model with nonlinear data lin.gpllm.tr <- bgpllm(X=X,Z=Znlin,XX=XX,bprior="b0",trace=TRUE, verb=0) plot(lin.gpllm.tr, main="FIT MIXTURE MODEL TO 1+2x+3x^2", layout="surf") pb <- mean(lin.gpllm.tr$trace$linarea$la) cat("POSTERIOR PROBABILITY OF LINEAR MODEL ",100*pb,"PERCENT\n") scan(quiet=TRUE) dev.off()