The HPPANEL Procedure

One-Way Fixed-Effects Model

The specification for the one-way fixed-effects model is

\[ u_{it}={\gamma }_{i} + {\epsilon }_{it} \]

where the ${{\gamma }_{i}}$ are nonrandom parameters to be estimated.

Let ${\mb{Q} _{0}=\mr{diag}(\mb{E} _{T_{i}})}$, with ${{\bar{\mb{J}}}_{T_{i}}=\mb{J} _{T_{i}}/\mi{T} _{i}}$ and $\mb{E}_{T_{i}}=\mb{I} _{T_{i}}-\bar{\mb{J}}_{T_{i}}$, where $\mb{J} _{T_{i}}$ is a matrix of $T_\mi {i} $ ones.

The matrix $\mb{Q} _{0}$ represents the within transformation. In the one-way model, the within transformation is the conversion of the raw data to deviations from a cross section’s mean. The vector $\mi{\tilde{\mb{x}}_\mi {it} }$ is a row of the general matrix $\mb{X} _{s}$, where the subscripted s implies that the constant (column of ones) is missing.

Let ${\tilde{\mb{X} }_{s}=\mb{Q} _{0}\mb{X} _{s}}$ and ${\tilde{\mb{y} }=\mb{Q} _{0}\mb{y} }$. The estimator of the slope coefficients is given by

\[ \tilde{{\beta }}_{s}=( \tilde{\mb{X} }^{'}_{s} \tilde{\mb{X} }_{s})^{-1} \tilde{\mb{X} }^{'}_{s}\tilde{\mb{y} } \]

After the slope estimates have been calculated, the estimation of an intercept or the cross-sectional fixed effects is handled as follows. First, you obtain the cross-sectional effects:

\[ \gamma _\mi {i} = \bar{\mi{y}}_\mi {i \cdot } - \tilde{{\beta }}_{s}\bar{\mb{\mi{x}}}_\mi {i \cdot } \quad \mbox{for}\quad \mi{i} = 1 \ldots \mi{N} \]

If the NOINT option is specified, then the dummy variables’ coefficients are set equal to the fixed effects. If you want an intercept, then the ith dummy variable is obtained from the following expression:

\[ D_\mi {i} = \gamma _\mi {i}- \gamma _\mi {N} \quad \mbox{for}\quad \mi{i} = 1 \ldots \mi{N} -1 \]

The intercept is the Nth fixed effect $\gamma _\mi {N} $.

The within-model sum of squared errors is

\[ \mr{SSE} = \sum _\mi {i = 1} ^\mi {N} \sum _\mi {t = 1} ^\mi {T_\mi {i}} (y_\mi {it} - \gamma _\mi {i} - \mb{X} _{s}\tilde{{\beta }}_{s})^{2} \]

The estimated error variance can be written as

\[ \hat{{\sigma }}_{{\epsilon }}^{2}= \mr{SSE} / (\mi{M}-\mi{N}-(\mi{K} -1)) \]

Alternatively, an equivalent way to express the error variance is

\[ \hat{{\sigma }}_{{\epsilon }}^{2}=\tilde{\mb{u} }^{'}\mb{Q} _{0} \tilde{\mb{u} } / (\mi{M}-\mi{N}-(\mi{K} -1)) \]

where the residuals ${\tilde{\mb{u} }}$ are given by $\tilde{\mb{u}}= (\mb{I} _{M}-\mb{j} _{M} \mb{j} {’}_{M}/\mi{M} )(\mb{y} -\mb{X} _{s}{\tilde{\beta }}_{s})$ if there is an intercept and by ${\tilde{\mb{u} }= (\mb{y}-\mb{X} _{s}\tilde{\beta }_{s}) }$ if there is not. The drawback is that the formula changes (but the results do not) with the inclusion of a constant.

The variance covariance matrix of ${\tilde{\beta }}_{s}$ is given by

\[ \mr{Var}\left[{\tilde{\beta }}_{s}\right] = \hat{{\sigma }}_{{\epsilon }}^{2}(\tilde{\mb{X} }^{'}_{s} \tilde{\mb{X} }_{s})^{-1} \]

The covariance of the dummy variables and the dummy variables with the ${\tilde{\beta }}_{s}$ depends on whether the intercept is included in the model. For more information, see the section One-Way Fixed-Effects Model in SAS/ETS 14.1 User's Guide.

Alternatively, the FIXONETIME model option estimates a one-way model in which the heterogeneity comes from time effects. This option is analogous to re-sorting the data by time and then by cross section, and then running a FIXONE model. The advantage of using the FIXONETIME option is that sorting is avoided and the model remains labeled correctly.