testcoef.rrenv.apweights: Hypothesis test of the coefficients of the reduced rank...

In Renvlp: Computing Envelope Estimators

Hypothesis test of the coefficients of the reduced rank envelope model with nonconstant error variance

Description

This function tests the null hypothesis L * beta * R = A versus the alternative hypothesis L * beta * R ~= A, where beta is estimated under the reduced rank envelope model with nonconstant error variance.

Usage

testcoef.rrenv.apweights(m, L, R, A)

Arguments

m	A list containing estimators and other statistics inherited from rrenv.apweights.
L	The matrix multiplied to beta on the left. It is a d1 by r matrix, while d1 is less than or equal to r.
R	The matrix multiplied to beta on the right. It is a p by d2 matrix, while d2 is less than or equal to p.
A	The matrix on the right hand side of the equation. It is a d1 by d2 matrix.

Note that inputs L, R and A must be matrices, if not, use as.matrix to convert them.

Details

This function tests for hypothesis H0: L beta R = A, versus Ha: L beta R != A. The beta is estimated by the reduced rank envelope model that accommodates nonconstant error variance. If L = Ir, R = Ip and A = 0, then the test is equivalent to the standard F test on if beta = 0. The test statistic used is vec(L beta R - A) hatSigma^-1 vec(L beta R - A)^T, where beta is the envelope estimator and hatSigma is the estimated asymptotic covariance of vec(L beta R - A). The reference distribution is chi-squared distribution with degrees of freedom d1 * d2.

Value

The output is a list that contains following components.

chisqStatistic	The test statistic.
dof	The degrees of freedom of the reference chi-squared distribution.
pValue	p-value of the test.
covMatrix	The covariance matrix of vec(L beta R).

Examples

data(vehicles)
X <- vehicles[, 1:11]  
Y <- vehicles[, 12:15]
X <- scale(X)
Y <- scale(Y)  # The scales of Y are vastly different, so scaling is reasonable here
## Not run: m <- rrenv.apweights(X, Y, 3, 2)
## Not run: m

L <- diag(4)
R <- matrix(1, 11, 1)
A <- matrix(0, 4, 1)

## Not run: test.res <- testcoef.rrenv.apweights(m, L, R, A)
## Not run: test.res

Renvlp documentation built on Oct. 11, 2023, 1:06 a.m.