The SYSLIN Procedure
Example 35.1 Klein’s Model I Estimated with LIML and 3SLS
(View the complete code for this example.)
This example uses PROC SYSLIN to estimate the classic Klein Model I. For a discussion of this model, see Theil (1971). The following statements read the data:
*---------------------------Klein's Model I----------------------------*
| By L.R. Klein, Economic Fluctuations in the United States, 1921-1941 |
| (1950), NY: John Wiley. A macro-economic model of the U.S. with |
| three behavioral equations, and several identities. See Theil, p.456.|
*----------------------------------------------------------------------*;
data klein;
input year c p w i x wp g t k wsum;
date=mdy(1,1,year);
format date monyy.;
y =c+i+g-t;
yr =year-1931;
klag=lag(k);
plag=lag(p);
xlag=lag(x);
label year='Year'
date='Date'
c ='Consumption'
p ='Profits'
w ='Private Wage Bill'
i ='Investment'
k ='Capital Stock'
y ='National Income'
x ='Private Production'
wsum='Total Wage Bill'
wp ='Govt Wage Bill'
g ='Govt Demand'
i ='Taxes'
klag='Capital Stock Lagged'
plag='Profits Lagged'
xlag='Private Product Lagged'
yr ='YEAR-1931';
datalines;
1920 . 12.7 . . 44.9 . . . 182.8 .
1921 41.9 12.4 25.5 -0.2 45.6 2.7 3.9 7.7 182.6 28.2
1922 45.0 16.9 29.3 1.9 50.1 2.9 3.2 3.9 184.5 32.2
1923 49.2 18.4 34.1 5.2 57.2 2.9 2.8 4.7 189.7 37.0
1924 50.6 19.4 33.9 3.0 57.1 3.1 3.5 3.8 192.7 37.0
1925 52.6 20.1 35.4 5.1 61.0 3.2 3.3 5.5 197.8 38.6
1926 55.1 19.6 37.4 5.6 64.0 3.3 3.3 7.0 203.4 40.7
1927 56.2 19.8 37.9 4.2 64.4 3.6 4.0 6.7 207.6 41.5
1928 57.3 21.1 39.2 3.0 64.5 3.7 4.2 4.2 210.6 42.9
1929 57.8 21.7 41.3 5.1 67.0 4.0 4.1 4.0 215.7 45.3
1930 55.0 15.6 37.9 1.0 61.2 4.2 5.2 7.7 216.7 42.1
... more lines ...
The following statements estimate the Klein model using the limited information maximum likelihood method. In addition, the parameter estimates are written to a SAS data set with the OUTEST= option.
proc syslin data=klein outest=b liml; endogenous c p w i x wsum k y; instruments klag plag xlag wp g t yr; consume: model c = p plag wsum; invest: model i = p plag klag; labor: model w = x xlag yr; run;
proc print data=b; run;
The PROC SYSLIN estimates are shown in Output 35.1.1 through Output 35.1.3.
Output 35.1.1: LIML Estimates for Consumption
Output 35.1.2: LIML Estimates for Investments
Output 35.1.3: LIML Estimates for Labor
The OUTEST= data set is shown in part in Output 35.1.4. Note that the data set contains the parameter estimates and root mean squared errors, _SIGMA_, for the first-stage instrumental regressions as well as the parameter estimates and for the LIML estimates for the three structural equations.
Output 35.1.4: The OUTEST= Data Set
The following statements estimate the model using the 3SLS method. The reduced form estimates are produced by the REDUCED option; IDENTITY statements are used to make the model complete.
proc syslin data=klein 3sls reduced; endogenous c p w i x wsum k y; instruments klag plag xlag wp g t yr; consume: model c = p plag wsum; invest: model i = p plag klag; labor: model w = x xlag yr; product: identity x = c + i + g; income: identity y = c + i + g - t; profit: identity p = y - w; stock: identity k = klag + i; wage: identity wsum = w + wp; run;
The preliminary 2SLS results and estimated cross-model covariance matrix are not shown. The 3SLS estimates are shown in Output 35.1.5 through Output 35.1.7. The reduced form estimates are shown in Output 35.1.8 through Output 35.1.11.
Output 35.1.5: 3SLS Estimates for Consumption
Output 35.1.6: 3SLS Estimates for Investments
Output 35.1.7: 3SLS Estimates for Labor
Output 35.1.8: Reduced Form Estimates
Output 35.1.9: Reduced Form Estimates
Output 35.1.10: Reduced Form Estimates
| Inverse Endogenous Variables | ||||||||
|---|---|---|---|---|---|---|---|---|
| CONSUME | INVEST | LABOR | PRODUCT | INCOME | PROFIT | STOCK | WAGE | |
| c | 1.634654 | 0.634654 | 1.095657 | 0.438802 | 0.195852 | 0.195852 | 0 | 1.291509 |
| p | 0.972364 | 0.972364 | -0.34048 | -0.13636 | 1.108721 | 1.108721 | 0 | 0.768246 |
| w | 0.649572 | 0.649572 | 1.440585 | 0.576943 | 0.072629 | 0.072629 | 0 | 0.513215 |
| i | -0.01272 | 0.987282 | 0.004453 | 0.001783 | -0.0145 | -0.0145 | 0 | -0.01005 |
| x | 1.621936 | 1.621936 | 1.10011 | 1.440585 | 0.181351 | 0.181351 | 0 | 1.281461 |
| wsum | 0.649572 | 0.649572 | 1.440585 | 0.576943 | 0.072629 | 0.072629 | 0 | 1.513215 |
| k | -0.01272 | 0.987282 | 0.004453 | 0.001783 | -0.0145 | -0.0145 | 1 | -0.01005 |
| y | 1.621936 | 1.621936 | 1.10011 | 0.440585 | 1.181351 | 0.181351 | 0 | 1.281461 |
Output 35.1.11: Reduced Form Estimates
| Reduced Form | ||||||||
|---|---|---|---|---|---|---|---|---|
| Intercept | plag | klag | xlag | yr | g | t | wp | |
| c | 46.7273 | 0.746307 | -0.12366 | 0.198633 | 0.163991 | 0.634654 | -0.19585 | 1.291509 |
| p | 42.77363 | 0.893474 | -0.18946 | -0.06173 | -0.05096 | 0.972364 | -1.10872 | 0.768246 |
| w | 31.57207 | 0.596871 | -0.12657 | 0.261165 | 0.215618 | 0.649572 | -0.07263 | 0.513215 |
| i | 27.6184 | 0.744038 | -0.19237 | 0.000807 | 0.000667 | -0.01272 | 0.014501 | -0.01005 |
| x | 74.3457 | 1.490345 | -0.31603 | 0.19944 | 0.164658 | 1.621936 | -0.18135 | 1.281461 |
| wsum | 31.57207 | 0.596871 | -0.12657 | 0.261165 | 0.215618 | 0.649572 | -0.07263 | 1.513215 |
| k | 27.6184 | 0.744038 | 0.80763 | 0.000807 | 0.000667 | -0.01272 | 0.014501 | -0.01005 |
| y | 74.3457 | 1.490345 | -0.31603 | 0.19944 | 0.164658 | 1.621936 | -1.18135 | 1.281461 |