Example 11.2 Analyzing Demand for Liquid Assets: Random Effects
Feige (1964) provides data on the demand for liquid assets. The data are for six states and the District of Columbia (CA, DC, FL, IL, NY, TX, and WA) and were collected each year from 1949 to 1959. All variables are log-transformed.
The following statements create the Assets data table. These statements assume that your libref is named mylib, but you can substitute any appropriately defined libref.
data mylib.Assets;
length state $ 2;
input state $ year d t s y rd rt rs;
label d = 'Per Capita Demand Deposits'
t = 'Per Capita Time Deposits'
s = 'Per Capita S & L Association Shares'
y = 'Permanent Per Capita Personal Income'
rd = 'Service Charge on Demand Deposits'
rt = 'Interest on Time Deposits'
rs = 'Interest on S & L Association Shares';
datalines;
CA 1949 6.2785 6.1924 4.4998 7.2056 -1.0700 0.1080 1.0664
CA 1950 6.4019 6.2106 4.6821 7.2889 -1.0106 0.1501 1.0767
CA 1951 6.5058 6.2729 4.8598 7.3827 -1.0024 0.4008 1.1291
CA 1952 6.4785 6.2729 5.0039 7.4000 -0.9970 0.4492 1.1227
CA 1953 6.4118 6.2538 5.1761 7.4200 -0.8916 0.4662 1.2110
CA 1954 6.4520 6.2971 5.3613 7.4478 -0.6951 0.4756 1.1924
... more lines ...
The data contain per capita consumptions for three liquid assets: demand deposits such as checking, time deposits, and savings and loan (S&L) shares. You posit a linear model for per capita demand deposits, with random effects for states.
The following statements fit a one-way random-effects model:
proc cpanel data = mylib.Assets;
id state year;
model d = y rd rt rs / ranone;
run;
The regression results are provided in Output 11.2.1.
The "Variance Component Estimates" table provides the estimated variances and standard deviations for the cross-sectional (state) effects in addition to the overall errors. A majority of the overall error variance can be attributed to differences between states, not differences within states.
The "Hausman Test for Random Effects" table shows the result of a Hausman specification test. The null hypothesis is that state effects can be treated as random (random-effects model) and that they do not need to be estimated directly (fixed-effects model). The test results favor the random-effects specification that is used to generate this output.
Output 11.2.1: Demand for Demand Deposits, One-Way Random-Effects Model
The CPANEL Procedure
One-Way Random Effects (RanOne)
Swamy and Arora Variance Components
Dependent Variable: d (Per Capita Demand Deposits)
| RanOne |
|---|
| Model Based |
|---|
| ASSETS |
|---|
| 77 |
|---|
| 7 |
|---|
| 11 |
|---|
| 0.0939 | 72 |
| 0.0013 | 0.0361 |
| 0.7651 | |
| 0.046394 | 0.215393 |
| 0.00134 | 0.036608 |
| 1 | -1.64325 | 0.6841 | -2.40 | 0.0189 | Intercept |
| 1 | 1.136776 | 0.0998 | 11.39 | <.0001 | Permanent Per Capita Personal Income |
| 1 | -0.26428 | 0.0517 | -5.11 | <.0001 | Service Charge on Demand Deposits |
| 1 | 0.032076 | 0.0291 | 1.10 | 0.2746 | Interest on Time Deposits |
| 1 | -0.41213 | 0.1192 | -3.46 | 0.0009 | Interest on S & L Association Shares |
The parameter estimate for the variable Y is greater than 1, indicating that demand is elastic to income—income has a more than proportional positive association with the demand for demand deposits. The coefficient on the variable RD indicates that demand deposits increase significantly as the service charge is reduced.
The variables RT and RS represent positive aspects of competing products, and you would expect these variables to affect demand negatively. The coefficient for RS meets that expectation, but the coefficient for RT is not significant.
The previous analysis used the default Swamy-Arora method to estimate the variance components. The CPANEL procedure supports four other methods, and you might be interested in how the different methods affect the analysis.
The following statements fit the model by using all five methods and include a COMPARE statement to compare the results:
proc cpanel data = mylib.Assets;
id state year;
sa: model d = y rd rt rs / ranone vcomp = sa;
wh: model d = y rd rt rs / ranone vcomp = wh;
wk: model d = y rd rt rs / ranone vcomp = wk;
fb: model d = y rd rt rs / ranone vcomp = fb;
nl: model d = y rd rt rs / ranone vcomp = nl;
compare / mstat(varcs varerr);
run;
The tables produced by the COMPARE statement are provided in Output 11.2.2.
Output 11.2.2: Comparison of Variance-Component Methods, Assets Data
The CPANEL Procedure
Model Comparison
Model: NL
Dependent Variable: d (Per Capita Demand Deposits)
| 0.0464 | 0.0315 | 0.0315 | 0.0291 | 0.0327 |
| 0.001340 | 0.000107 | 0.001340 | 0.001340 | 0.001149 |
You conclude that how you estimate variance components has little bearing on the regression results.