The GAM Procedure |
MODEL Statement |
The MODEL statement specifies the dependent variable and the independent effects you want to use to model its values. Specify the independent parametric variables inside the parentheses of PARAM( ). The parametric variables can be either classification variables or continuous variables. Classification variables must be declared with a CLASS statement. Interactions between variables can also be included as parametric effects. Multiple PARAM() statements are allowed in the MODEL statement. The syntax for the specification of effects is the same as for the GLM procedure ( Chapter 39, The GLM Procedure ).
Only continuous variables can be specified in smoothing effects. Any number of smoothing effects can be specified, as follows:
Smoothing Effect |
Meaning |
---|---|
SPLINE(variable <, DF=number>) |
fit a smoothing spline with the variable and with DF=number |
LOESS(variable <, DF=number>) |
fit a local regression with the variable and with DF=number |
SPLINE2(variable1, variable2 <,DF=number>) |
fit a bivariate thin-plate smoothing spline with variable1 and variable2 and with DF=number |
The number specified in the DF= option must be positive. If you specify neither the DF= option nor the METHOD=GCV in the MODEL statement, then the default used is DF=4. Note that for univariate spline components, a degree of freedom is removed by default to account for the linear portion of the model, so the value displayed in the “Fit Summary” and “Analysis of Deviance tables” will be one less than the value you specify.
Both parametric effects and smoothing effects are optional, but at least one of them must be present.
If only parametric variables are present, PROC GAM fits a parametric linear model by using the terms inside the parentheses of PARAM( ). If only smoothing effects are present, PROC GAM fits a nonparametric additive model. If both types of effect are present, PROC GAM fits a semiparametric model by using the parametric effects as the linear part of the model.
Table 36.1 shows how to specify various models for a dependent variable y and independent variables x, x1, and x2.
Type of Model |
Syntax for model |
Mathematical Form |
Parametric |
y=param(x1 x2) |
|
Nonparametric |
y=spline(x) |
|
Nonparametric |
y=loess(x) |
|
Semiparametric |
y=spline(x1) param(x2) |
|
Additive |
y=spline(x1) spline(x2) |
|
Thin-plate spline |
y=spline2(x1,x2) |
|
Response variable options determine how the GAM procedure models probabilities for binary data.
You can specify the following options by enclosing them in parentheses after the response variable. See the section CLASS Statement for more detail.
reverses the order of the response categories. If both the DESCENDING and ORDER= options are specified, PROC GAM orders the response categories according to the ORDER= option and then reverses that order.
designates the first ordered category as the event.
designates the last ordered category as the event.
model Y (event='1') = X;
specifies the sort order for the levels of the response variable. When ORDER=FORMATTED (the default) for numeric variables for which you have supplied no explicit format (that is, for which there is no corresponding FORMAT statement in the current PROC GAM run or in the DATA step that created the data set), the levels are ordered by their internal (numeric) value. If you specify the ORDER= option in the MODEL statement and the ORDER= option in the CLASS statement, the former takes precedence. The following table shows the interpretation of the ORDER= values:
Value of ORDER= |
Levels Sorted By |
---|---|
DATA |
order of appearance in the input data set |
FORMATTED |
external formatted value, except for numeric variables with no explicit format, which are sorted by their unformatted (internal) value |
FREQ |
descending frequency count; levels with the most observations come first in the order |
INTERNAL |
unformatted value |
By default, ORDER=FORMATTED. For the FORMATTED and INTERNAL values, the sort order is machine-dependent.
For more information about sort order, see the chapter on the SORT procedure in the Base SAS Procedures Guide and the discussion of BY-group processing in SAS Language Reference: Concepts.
specifies the significance level of the confidence limits on the final nonparametric component estimates when you request confidence limits to be included in the output data set. Specify number as a value between 0 and 1. The default value is 0.05. See the section OUTPUT Statement for more information about the OUTPUT statement.
specifies the type of method to be used to produce the "Analysis of Deviance" table for smoothing effects. The available choices are as follows:
specifies that PROC GAM perform tests by fitting nested GAM models. This is the default choice if you do not specify the ANODEV= option. This choice requires fitting separate GAM models where one smoothing term is omitted from each model.
specifies that PROC GAM perform approximate tests of smoothing effects. To test each smoothing effect, a weighted least squares model is fitted to the remaining parametric part of the model while keeping other nonlinear smoothers fixed. For details, see Hastie (1991). This choice requires only a single GAM fitting to be performed, which reduces the time of the procedure.
requests that the procedure not produce the "Analysis of Deviance" table for smoothing effects.
specifies the distribution family used in the model. The choices for distribution-id are displayed in Table 36.2. See Distribution Family and Canonical Link for more information.
DIST= |
Distribution |
Response Data Type |
GAUSSIAN | GAUS | NORM |
Normal (Gaussian) |
continuous variables |
BINOMIAL | LOGI | BIN |
Binomial |
binary variables |
POISSON | POIS | LOGL |
Poisson |
nonnegative discrete variables |
GAMMA | GAMM |
Gamma |
positive continuous variables |
IGAUSSIAN | IGAU | INVG |
Inverse gaussian |
positive continuous variables |
Canonical link functions are used with those distributions. Although theoretically, alternative links are possible, with nonparametric regression models the final fit is relatively insensitive to the precise choice of link functions. Therefore, only the canonical link for each distribution family is implemented in PROC GAM. The loess smoother is not available for DIST=BINOMIAL when the number of trials is greater than 1.
specifies the convergence criterion for the backfitting algorithm. The default value is 1E.
specifies the convergence criterion for the local scoring algorithm. The default value is 1E.
produces an iteration summary table for the smoothing effects when doing backfitting and local scoring.
specifies the maximum number of iterations for the backfitting algorithm. The default value is .
specifies the maximum number of iterations for the local scoring algorithm. The default value is .
specifies that the value of the smoothing parameter should be selected by generalized cross validation. If you specify both METHOD=GCV and the DF= option for the smoothing effects, the user-specified DF is used, and the METHOD=GCV option is ignored. See the section Selection of Smoothing Parameters for more details on the GCV method.
Copyright © 2009 by SAS Institute Inc., Cary, NC, USA. All rights reserved.