## ----include=FALSE------------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
library(gamutil)
mdl <- gamutil:::mdl
RhpcBLASctl::blas_set_num_threads(1)
RhpcBLASctl::omp_set_num_threads(1)

## ----setup--------------------------------------------------------------------
library(gamutil)

## -----------------------------------------------------------------------------
print(mdl)

## ----fig.dim=c(5,5)-----------------------------------------------------------
library(mgcv)
vis.gam(mdl, view=c('x1','x2'), plot.type='contour')

## ----fig.dim=c(6,5)-----------------------------------------------------------
plot_contour(mdl, view=c('x1','x2'))

## ----fig.dim=c(6,5)-----------------------------------------------------------
lbr <- seq(-10, 10, by=1.0)
cbr <- seq(-10, 10, by=0.1)
plot_contour(mdl, view=c('x1','x2'),
	     contour.line.breaks=lbr, contour.color.breaks=cbr)

## ----fig.dim=c(6,5)-----------------------------------------------------------
plot.gam(mdl, select=6)

## ----fig.dim=c(6,5)-----------------------------------------------------------
lbr <- seq(-5, 5, by=0.50)
cbr <- seq(-5, 5, by=0.05)
plot_contour(mdl, view=c('x1','x2'), summed=FALSE, terms.size='min',
	     verbose=TRUE, contour.line.breaks=lbr, contour.color.breaks=cbr)

## ----fig.dim=c(6,5)-----------------------------------------------------------
plot_contour(mdl, view=c('x1','x2'), summed=FALSE, terms.size='medium',
	     verbose=TRUE)

## ----fig.dim=c(6,5)-----------------------------------------------------------
plot_contour(mdl, view=c('x1','x2'), summed=FALSE, terms.size='max',
	     verbose=TRUE)

## ----fig.dim=c(6,5)-----------------------------------------------------------
library(ggplot2)
plt <- plot_contour(mdl, view=c('x1','x2'), summed=FALSE, terms.size='medium',
		    verbose=TRUE)
plt <- plt + theme(text=element_text(size=25))
plt