The call that produced the object x is printed, followed by a five-column
matrix with columns NonZero, Lin, NonLin,
columns say how many nonzero, linear and nonlinear terms there are.
the percent deviance explained (relative to the null deviance).
# NOT RUN {set.seed(1)
n <- 100; p <- 12
x <- matrix(rnorm(n * p), n, p)
beta <- matrix(c(rep(2, 3), rep(0, 9)), ncol = 1)
y <- x %*% beta + x[, 4]^2 + rnorm(n)
fit <- rgam(x, y)
print(fit)
# }