getCharSTM2u: Computing Some Overall Characteristics in 'compContourM2u'

View source: R/getCharSTM2u.R

getCharSTM2uR Documentation

Computing Some Overall Characteristics in compContourM2u

Description

The function computes some overall characteristics of directional regression quantiles in the output of compContourM2u, namely the list COutST$CharST. It makes possible to obtain some useful information without saving any file on the disk, and it can be easily modified by the users according to their wishes.

Usage

getCharSTM2u(Tau, N, M, P, BriefDQMat, CharST, IsFirst)

Arguments

Tau

the quantile level in (0, 0.5).

N

the number of observations.

M

the dimension of responses.

P

the dimension of regressors including the intercept.

BriefDQMat

the method-specific matrix containing the rows of a potential individual output file corresponding to CTechST$BriefOutputI = 1. See the details below.

CharST

the output list, updated with each run of the function.

IsFirst

the indicator equal to one in the first run of getCharSTM2u (when CharST is initialized) and equal to zero otherwise.

Details

This function is called inside compContourM2u. First, it is called with BriefDQMat = NULL, CharST = NULL and IsFirst = 1 to initialize the output list CharST, and then it is called with IsFirst = 0 successively for the content of each potential output file corresponding to CTechST$BriefOutputI = 1, i.e., even if the output file(s) are not stored on the disk owing to CTechST$OutSaveI = 0.

It still remains to describe in detail the content of possible output files, describing the optimal conic segmentation of the directional space that lies behind the optimization problem involved.

If CTechST$BriefOutputI = 1, then the rows of such files are vectors of length 1+1+M+P*M+M of the form c(ConeID, Nu, UVec, vec(ACOMat), MuBRow) where

ConeID

is the number/order of the cone related to the line. If M > 2, then a cone can appear in the output repeatedly (under different numbers).

Nu

is the number of negative residuals corresponding to the interior directions of the cone.

UVec

is a normalized vector of the cone. It is usually its vertex direction but it may also be its interior vector pointing to a vertex of the artificial intersection of the cone with the bounding box [-1,1]^M. The max normalization is used if the breadth-first search algorithm is employed and the L2 normalization is used in the other case (when M = 2 and CTechST$D2SpecI = 1).

ACOMat

is the matrix describing AVec, AVec = ACOMat%*%UVec.

MuBRow

is the constant vector of the Lagrange multipliers corresponding to BVec. Its inner product with UVec is equal to the optimal value Psi of the objective function for that direction.

Recall that c(BVec, AVec) stands for the coefficients of the regression quantile hyperplane associated with UVec and always BVec = UVec.

If CTechST$BriefOutputI = 0, then the rows of the potential output file(s) are longer (of length 1+1+P*M+M+P+P) because they contain two more vectors appended at the end. The rows are of the form c(ConeID, Nu, UVec, vec(ACOMat), MuBRow, MuR0Row, IZ) where

MuRORow

is the constant vector of the Lagrange multipliers corresponding to zero residuals associated with the interior of the cone. That is to say that all directions from the interior of the cone result in the regression Tau-quantile hyperplanes containing the same P observations.

IZ

is the vector containing original indices of the P observations with zero residuals for all directions from the interior of the cone.

Value

getCharSTM2u returns a list with the following components:

NUESkip

the number of (skipped) directions (and corresponding hyperplanes) artificially induced by intersecting the cones with the [-1,1]^M bounding box

NAZSkip

the number of (skipped) hyperplanes (and corresponding directions) not counted in NUESkip and with at least one coordinate of AVec zero.

NBZSkip

the number of (skipped) hyperplanes (and corresponding directions) not counted in NUESkip and with at least one coordinate of BVec zero.

HypMat

(for M > 4) the component is missing
(for M <= 4) the matrix with M + P columns containing (in rows) all the distinct regression Tau-quantile hyperplane coefficients c(BVec, AVec) rounded to the eighth decimal digit and sorted lexicographically. This matrix can be used for the computation of the regression Tau-quantile contour.

CharMaxMat

the matrix with the (slightly rounded) maxima of certain directional regression Tau-quantile characteristics over all remaining vertex directions.
If P = 1, then CharMaxMat has only two rows:
c(UVec, max(Psi)), and
c(UVec, max(|MuBRow|)),
respectively.
If P > 1, then the rows of CharMaxMat are as follows:
c(UVec, max(|Psi|)),
c(UVec, max(MuBRow)),
c(UVec, max(|c(a_2,...,a_P)|)),
c(UVec, max(|a_2|)),
...,
c(UVec, max(|a_P|)),
respectively. If P = 2, then the last row is missing for not being included twice.

CharMinMat

the matrix with the (slightly rounded) minima of certain directional regression Tau-quantile characteristics over all remaining vertex directions.
If P = 1, then CharMinMat has only two rows:
c(UVec, min(Psi)), and
c(UVec, min(|MuBRow|)),
respectively.
If P > 1, then CharMinMat has three rows:
c(UVec, min(Psi)),
c(UVec, min(|MuBRow|)), and
c(UVec, min(|c(a_2,...,a_P)|)),
respectively.

Note that || symbolizes the Euclidean norm, and that the vertices (UVec) in the rows of CharMaxMat and CharMinMat are generally different and denote (one of) the direction(s) where the row maximum or minimum is attained.

Examples

##Run print(getCharSTM2u) to examine the default setting.

modQR documentation built on May 11, 2022, 5:18 p.m.

Related to getCharSTM2u in modQR...