Joint Hypothesis Tests with Complex Alternatives, the Chi-Bar-Square Statistic |
Silvapulle and Sen (2004) propose a test statistic for testing hypotheses where the null or the alternative hypothesis or both involve inequalities. You can test special cases of these hypotheses with the JOINT option in the ESTIMATE and the LSMESTIMATE statement. Consider the estimable functions
and the hypotheses
and
. The alternative hypothesis defines a convex cone
at the origin. Suppose that under the null hypothesis
follows a multivariate normal distribution with mean
and variance
. The restricted alternative prevents you from using the usual
or chi-square test machinery, since the distribution of the test statistic under the alternative might not follow the usual rules. Silvapulle and Sen (2004) coined a statistic that takes into account the projection of the observed estimate onto the convex cone formed by the alternative parameter space. This test statistic is called the chi-bar-square statistic, and
-values are obtained by simulation; see, in particular, Chapter 3.4 in Silvapulle and Sen (2004).
Briefly, let be a multivariate normal random variable with mean
and variance matrix
. The chi-bar-square statistic is the random variable
![]() |
![]() |
![]() |
|||
![]() |
![]() |
![]() |
and it can be motivated by a geometric argument. The quadratic form in is the
-projection of
onto the cone
. Suppose that this projected point is
. If
, then
and
. If
is completely outside of the cone
, then
is a point on the surface of the cone. Similarly,
is the length of the segment from the origin to
in the
-space with norm
. If you apply the Pythagorean theorem, you can see that the chi-bar-square statistic measures the length of the segment from the origin to the projected point
in
.
To calculate -values for chi-bar-square statistics, a simulation-based approach is taken. Consider again the set of
estimable functions
with estimate
and variance
.
First, the observed value of the statistic is computed as
![]() |
Then, independent random samples
are drawn from an
distribution and the following chi-bar-statistics are computed for the sample:
![]() |
![]() |
![]() |
|||
![]() |
![]() |
![]() |
|||
![]() |
![]() |
![]() |
The -value is estimated by the fraction of simulated statistics that are greater than or equal to the observed value
.
Notice that unless is interior to the cone
, finding the value of
requires the solution to a quadratic optimization problem. When
is large, or when many simulations are requested, the computation of
-values for chi-bar-square statistics might require considerable computing time.