The HPQLIM Procedure

HETERO Statement

HETERO dependent variables $\sim $ exogenous variables </ options > ;

The HETERO statement specifies variables that are related to the heteroscedasticity of the residuals and the way that these variables are used to model the error variance. PROC HPQLIM supports the following heteroscedastic regression model:

\[  y_{i} = \mb {x} _{i}’{\bbeta } + {\epsilon }_{i}  \]
\[  {\epsilon }_{i} \sim \mr {N}(0, {\sigma }_{i}^{2})  \]

For more information about the specification of functional forms, see the section Heteroscedasticity. The following options specify the functional forms of heteroscedasticity:

LINK=EXP | LINEAR

specifies the functional form.

EXP

specifies the exponential link function:

$\displaystyle  \sigma _{i}^{2}  $
$\displaystyle  =  $
$\displaystyle  \sigma ^{2}(1+{\exp }( \Strong{z} ^{}_{i}{\bgamma }))  $

LINEAR

specifies the linear link function:

$\displaystyle  \sigma _{i}^{2}  $
$\displaystyle  =  $
$\displaystyle  \sigma ^{2}(1+\Strong{z} ^{}_{i}{\bgamma })  $

The default is LINK=EXP.

NOCONST

specifies that there be no constant in the linear or exponential heteroscedasticity model:

$\displaystyle  \sigma _{i}^{2}  $
$\displaystyle  =  $
$\displaystyle  \sigma ^{2}( \Strong{z} ^{}_{i}{\bgamma })  $
$\displaystyle \sigma _{i}^{2}  $
$\displaystyle  =  $
$\displaystyle  \sigma ^{2} {\exp }(\Strong{z} ^{}_{i}{\bgamma })  $

This option is ignored if you do not specify the LINK= option.

SQUARE

estimates the model by using the square of the linear heteroscedasticity function. For example, you can specify the following heteroscedasticity function:

\[  \sigma _{i}^{2} = \sigma ^2(1+(\mb {z} ^{}_{i}{\bgamma })^{2})  \]
   model y = x1 x2 / censored(lb=0);
   hetero y ~ z1 / link=linear square;

The SQUARE option does not apply to the exponential heteroscedasticity function because the square of an exponential function of $\mb {z} ^{}_{i}{\bgamma }$ is the same as the exponential of $2\mb {z} ^{}_{i}{\bgamma }$. Hence, the only difference is that all $\bgamma $ estimates are divided by two.

This option is ignored if you do not specify the LINK= option. You cannot use the HETERO statement within a Bayesian framework.