Chapter Contents |
Previous |
Next |
The MODEL Procedure |
Equality restrictions can be written as a vector function
Inequality restrictions are either active or inactive. When an inequality restriction is active, it is treated as an equality restriction. All inactive inequality restrictions can be written as a vector function
Strict inequalities, such as , are transformed into inequalities as ,where the tolerance is controlled by the EPSILON= option on the FIT statement and defaults to 10^{-8}. The ith inequality restriction becomes active if F_{i} < 0 and remains active until its lagrange multiplier becomes negative. Lagrange multipliers are computed for all the nonredundant equality restrictions and all the active inequality restrictions.
For the following, assume the vector contains all the current active restrictions. The constraint matrix A is
The covariance matrix for the restricted parameter estimates is computed as
where H is Hessian or approximation to the Hessian of the objective function ( for OLS), and Z is the last (np - nc) columns of Q. Q is from an LQ factorization of the constraint matrix, nc is the number of active constraints, and np is the number of parameters. Refer to Gill, Murray, and Wright (1981) for more details on LQ factorization. The covariance column in Table 14.1 summarizes the Hessian approximation used for each estimation method.
The covariance matrix for the Lagrange multipliers is computed as
The p-value reported for a restriction is computed from a beta distribution rather than a t-distribution because the numerator and the denominator of the t-ratio for an estimated Lagrange multiplier are not independent.
The Lagrange multipliers for the active restrictions are printed with the parameter estimates. The Lagrange multiplier estimates are computed using the relationship
where the dimension of the constraint matrix A is the number of constraints by the number of parameters, is the vector of Lagrange multipliers, and g is the gradient of the objective function at the final estimates.
The final gradient includes the effects of the estimated S matrix. For example, for OLS the final gradient would be:
where r is the residual vector. Note that when nonlinear restrictions are imposed, the convergence measure R may have values greater than one for some iterations.
Chapter Contents |
Previous |
Next |
Top |
Copyright © 1999 by SAS Institute Inc., Cary, NC, USA. All rights reserved.