gamsel ObjectProduces plots of the estimated functions for specified variables at a given
value of lambda.
# S3 method for gamsel
plot(x, newx, index, which = 1:p, rugplot = TRUE, ylims, ...)Fitted gamsel object.
nobs_new x p matrix giving values of each
predictor at which to plot.
Index of lambda value (i.e., model) for which plotting is desired.
Which values to plot. Default is all variables, i.e.
{1,2,...,nvars}. Besides indices, which can take two special values:
"nonzero" will plot only the nonzero functions, and
"nonlinear" only the nonlinear functions.
If TRUE, a rugplot showing values of x is shown
at the bottom of each fitted function plot.
ylim argument for plotting each curve, which overides
the default which is the range of all the functions.
Optional graphical parameters to plot.
Alexandra Chouldechova and Trevor Hastie
Maintainer: Trevor
Hastie hastie@stanford.edu
A plot of the specified fitted functions is produced. Nothing is returned.
Chouldechova, A. and Hastie, T. (2015) Generalized Additive Model Selection
gamsel, and print.gamsel, summary.gamsel
##set.seed(1211)
##data=gamsel:::gendata(n=500,p=12,k.lin=3,k.nonlin=3,deg=8,sigma=0.5)
data = readRDS(system.file("extdata/gamsel_example.RDS", package = "gamsel"))
attach(data)
bases=pseudo.bases(X,degree=10,df=6)
# Gaussian gam
gamsel.out=gamsel(X,y,bases=bases)
par(mfrow=c(3,4))
plot(gamsel.out,newx=X,index=20)
Run the code above in your browser using DataLab