env: Fit the response envelope model

In Renvlp: Computing Envelope Estimators

Fit the response envelope model

Description

Fit the response envelope model in multivariate linear regression with dimension u.

Usage

env(X, Y, u, asy = TRUE, init = NULL)

Arguments

X	Predictors. An n by p matrix, p is the number of predictors. The predictors can be univariate or multivariate, discrete or continuous.
Y	Multivariate responses. An n by r matrix, r is the number of responses and n is number of observations. The responses must be continuous variables.
u	Dimension of the envelope. An integer between 0 and r.
asy	Flag for computing the asymptotic variance of the envelope estimator. The default is TRUE. When p and r are large, computing the asymptotic variance can take much time and memory. If only the envelope estimators are needed, the flag can be set to asy = FALSE.
init	The user-specified value of Gamma for the envelope subspace. An r by u matrix. The default is the one generated by function envMU.

Details

This function fits the envelope model to the responses and predictors,

Y = \mu + \Gamma\eta X+\varepsilon, \Sigma=\Gamma\Omega\Gamma'+\Gamma_{0}\Omega_{0}\Gamma'_{0}

using the maximum likelihood estimation. When the dimension of the envelope is between 1 and r-1, the starting value and blockwise coordinate descent algorithm in Cook et al. (2016) is implemented. When the dimension is r, then the envelope model degenerates to the standard multivariate linear regression. When the dimension is 0, it means that X and Y are uncorrelated, and the fitting is different.

Value

The output is a list that contains the following components:

beta	The envelope estimator of the regression coefficients.
Sigma	The envelope estimator of the error covariance matrix.
Gamma	An orthonormal basis of the envelope subspace.
Gamma0	An orthonormal basis of the complement of the envelope subspace.
eta	The coordinates of beta with respect to Gamma.
Omega	The coordinates of Sigma with respect to Gamma.
Omega0	The coordinates of Sigma with respect to Gamma0.
mu	The estimated intercept.
loglik	The maximized log likelihood function.
covMatrix	The asymptotic covariance of vec(beta). The covariance matrix returned are asymptotic. For the actual standard errors, multiply by 1 / n.
asySE	The asymptotic standard error for elements in beta under the envelope model. The standard errors returned are asymptotic, for actual standard errors, multiply by 1 / sqrt(n).
ratio	The asymptotic standard error ratio of the standard multivariate linear regression estimator over the envelope estimator, for each element in beta.
n	The number of observations in the data.

References

Cook, R. D., Li, B. and Chiaromente, F. (2010). Envelope Models for Parsimonious and Efficient Multivariate Linear Regression (with discussion). Statist. Sinica 20, 927- 1010.

Cook, R. D., Forzani, L. and Su, Z. (2016) A Note on Fast Envelope Estimation. Journal of Multivariate Analysis. 150, 42-54.

Examples

data(wheatprotein)
X <- wheatprotein[, 8]
Y <- wheatprotein[, 1:6]
u <- u.env(X, Y)
u

m <- env(X, Y, 1)
m
m$beta

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