Apply this function over all m, j not equal (order matters) to calculate each portion of the integral we need for the lower bound for calc_b2.
performIntegralUpper2(
x,
muVec,
sigMat,
lBounds1,
uBounds1,
lBounds2,
uBounds2,
lBounds3,
uBounds3,
lBounds4,
uBounds4
)
Apply over this 2*1 vector, the elements that will be the largest and second largest in magnitude.
Mean vector of test statistics under the alternative (assuming it's MVN).
Covariance matrix of test statistics under the alternative (assuming it's MVN).
A 3J-2 vector of lower bounds for the first integral (see paper).
A 3J-2 vector of upper bounds for the second integral (see paper).
A 3J-2 vector of lower bounds for the first integral (see paper).
A J3J-2 vector of upper bounds for the second integral (see paper).
A 3J-2 vector of lower bounds for the third integral (see paper).
A 3J-2 vector of upper bounds for the third integral (see paper).
A 3J-2 vector of lower bounds for the fourth integral (see paper).
A 3J-2 vector of upper bounds for the fourth integral (see paper).
The value of the integration.