The QLIM Procedure

Tests on Parameters

In general, the hypothesis tested can be written as

$\text{[math]}$

where $\text{[math]}$ is an r by 1 vector valued function of the parameters $\text{[math]}$ given by the r expressions specified in the TEST statement.

Let $\text{[math]}$ be the estimate of the covariance matrix of $\text{[math]}$ . Let $\text{[math]}$ be the unconstrained estimate of $\text{[math]}$ and $\text{[math]}$ be the constrained estimate of $\text{[math]}$ such that $\text{[math]}$ . Let

$\text{[math]}$

Using this notation, the test statistics for the three kinds of tests are computed as follows.

The Wald test statistic is defined as

$\text{[math]}$

The Wald test is not invariant to reparameterization of the model (Gregory 1985; Gallant 1987, p. 219). For more information about the theoretical properties of the Wald test, see Phillips and Park (1988).

The Lagrange multiplier test statistic is

$\text{[math]}$

where $\text{[math]}$ is the vector of Lagrange multipliers from the computation of the restricted estimate $\text{[math]}$ .

The likelihood ratio test statistic is

$\text{[math]}$

where $\text{[math]}$ represents the constrained estimate of $\text{[math]}$ and $\text{[math]}$ is the concentrated log-likelihood value.

For each kind of test, under the null hypothesis the test statistic is asymptotically distributed as a $\text{[math]}$ random variable with r degrees of freedom, where r is the number of expressions in the TEST statement. The p-values reported for the tests are computed from the $\text{[math]}$ distribution and are only asymptotically valid.

Monte Carlo simulations suggest that the asymptotic distribution of the Wald test is a poorer approximation to its small sample distribution than that of the other two tests. However, the Wald test has the lowest computational cost, since it does not require computation of the constrained estimate $\text{[math]}$ .

The following is an example of using the TEST statement to perform a likelihood ratio test:

   proc qlim;
      model y = x1 x2 x3;
      test x1 = 0, x2 * .5 + 2 * x3 = 0 /lr;
   run;