The MODEL procedure analyzes models in which the relationships among the variables form a system of one or more nonlinear equations. Primary uses of the MODEL procedure are estimation, simulation, and forecasting of nonlinear simultaneous equation models.
PROC MODEL features include the following:
SAS programming statements to define simultaneous systems of nonlinear equations
tools to analyze the structure of the simultaneous equation system
ARIMA, PDL, and other dynamic modeling capabilities
tools to specify and estimate the error covariance structure
tools to estimate and solve ordinary differential equations
the following methods of parameter estimation:
– ordinary least squares (OLS)
– twostage least squares (2SLS)
– seemingly unrelated regression (SUR) and iterative SUR (ITSUR)
– threestage least squares (3SLS) and iterative 3SLS (IT3SLS)
– generalized method of moments (GMM)
– simulated method of moments (SMM)
– full information maximum likelihood (FIML)
– general loglikelihood maximization
simulation and forecasting capabilities
Monte Carlo simulation
goalseeking solutions
A system of equations can be nonlinear in the parameters, nonlinear in the observed variables, or nonlinear in both the parameters and the variables. Nonlinear in the parameters means that the mathematical relationship between the variables and parameters is not required to have a linear form. (A linear model is a special case of a nonlinear model.) A general nonlinear system of equations can be written as












where is an endogenous variable, is an exogenous variable, is a parameter, and is the unknown error. The subscript represents time or some index to the data.
In econometrics literature, the observed variables are either endogenous (dependent) variables or exogenous (independent) variables. This system can be written more succinctly in vector form as

This system of equations is in general form because the error term is by itself on one side of the equality. Systems can also be written in normalized form by placing the endogenous variable on one side of the equality, with each equation defining a predicted value for a unique endogenous variable. A normalized form equation system can be written in vector notation as

PROC MODEL handles equations written in both forms.
Econometric models often explain the current values of the endogenous variables as functions of past values of exogenous and endogenous variables. These past values are referred to as lagged values, and the variable is called lag i of the variable x. Using lagged variables, you can create a dynamic, or timedependent, model. In the preceding model systems, the lagged exogenous and endogenous variables are included as part of the exogenous variables.
If the data are time series, so that t indexes time (see Chapter 3: Working with Time Series Data, for more information on time series), it is possible that depends on or, more generally, the ’s are not identically and independently distributed. If the errors of a model system are autocorrelated, the standard error of the estimates of the parameters of the system will be inflated.
Sometimes the ’s are not identically distributed because the variance of is not constant. This is known as heteroscedasticity. Heteroscedasticity in an estimated model can also inflate the standard error of the estimates of the parameters. Using a weighted estimation can sometimes eliminate this problem. Alternately, a variance model such as GARCH or EGARCH can be estimated to correct for heteroscedasticity. If the proper weighting scheme and the form of the error model is difficult to determine, generalized methods of moments (GMM) estimation can be used to determine parameter estimates with very few assumptions about the form of the error process.
Other problems can also arise when estimating systems of equations. Consider the following system of equations, which is nonlinear in its parameters and cannot be estimated with linear regression:






This system of equations represents a rudimentary predatorprey process with as the prey and as the predator (the second term in both equations is a logistics curve). The two equations must be estimated simultaneously because of the crossdependency of ’s. This crossdependency makes and violate the assumption of independence. Nonlinear ordinary least squares estimation of these equations produce biased and inconsistent parameter estimates. This is called simultaneous equation bias.
One method to remove simultaneous equation bias, in the linear case, is to replace the endogenous variables on the righthand side of the equations with predicted values that are uncorrelated with the error terms. These predicted values can be obtained through a preliminary, or “firststage,” instrumental variable regression. Instrumental variables, which are uncorrelated with the error term, are used as regressors to model the predicted values. The parameter estimates are obtained by a second regression by using the predicted values of the regressors. This process is called twostage least squares.
In the nonlinear case, nonlinear ordinary least squares estimation is performed iteratively by using a linearization of the model with respect to the parameters. The instrumental solution to simultaneous equation bias in the nonlinear case is the same as the linear case, except the linearization of the model with respect to the parameters is predicted by the instrumental regression. Nonlinear twostage least squares is one of several instrumental variables methods available in the MODEL procedure to handle simultaneous equation bias.
When you have a system of several regression equations, the random errors of the equations can be correlated. In this case, the largesample efficiency of the estimation can be improved by using a joint generalized least squares method that takes the crossequation correlations into account. If the equations are not simultaneous (no dependent regressors), then seemingly unrelated regression (SUR) can be used. The SUR method requires an estimate of the crossequation error covariance matrix, . The usual approach is to first fit the equations by using OLS, compute an estimate from the OLS residuals, and then perform the SUR estimation based on . The MODEL procedure estimates by default, or you can supply your own estimate of .
If the equation system is simultaneous, you can combine the 2SLS and SUR methods to take into account both simultaneous equation bias and crossequation correlation of the errors. This is called threestage least squares or 3SLS.
A different approach to the simultaneous equation bias problem is the full information maximum likelihood (FIML) estimation method. FIML does not require instrumental variables, but it assumes that the equation errors have a multivariate normal distribution. 2SLS and 3SLS estimation do not assume a particular distribution for the errors.
Other nonnormal error distribution models can be estimated as well. The centered t distribution with estimated degrees of freedom and nonconstant variance is an additional builtin likelihood function. If the distribution of the equation errors is not normal or t but known, then the log likelihood can be specified by using the ERRORMODEL statement.
Once a nonlinear model has been estimated, it can be used to obtain forecasts. If the model is linear in the variables you want to forecast, a simple linear solve can generate the forecasts. If the system is nonlinear, an iterative procedure must be used. The preceding example system is linear in its endogenous variables. The MODEL procedure’s SOLVE statement is used to forecast nonlinear models.
One of the main purposes of creating models is to obtain an understanding of the relationship among the variables. There are usually only a few variables in a model you can control (for example, the amount of money spent on advertising). Often you want to determine how to change the variables under your control to obtain some target goal. This process is called goal seeking. PROC MODEL allows you to solve for any subset of the variables in a system of equations given values for the remaining variables.
The nonlinearity of a model creates two problems with the forecasts: the forecast errors are not normally distributed with zero mean, and no formula exists to calculate the forecast confidence intervals. PROC MODEL provides Monte Carlo techniques, which, when used with the covariance of the parameters and error covariance matrix, can produce approximate error bounds on the forecasts. The following distributions on the errors are supported for multivariate Monte Carlo simulation:
Cauchy
chisquared
empirical
F
Poisson
t
uniform
A transformation technique is used to create a covariance matrix for generating the correct innovations in a Monte Carlo simulation.