subselect: xi2 – R documentation

Become an expert in R — Interactive courses, Cheat Sheets, certificates and more!

xi2

Computes the Xi squared coefficient for a multivariate linear hypothesis

Description

Computes the Xi squared index of "effect magnitude". The maximization of this criterion is equivalent to the maximization of the traditional test statistic, the Bartllet-Pillai trace.

Usage

xi2.coef(mat, H, r, indices,
tolval=10*.Machine$double.eps, tolsym=1000*.Machine$double.eps)

Arguments

`mat`	the Variance or Total sums of squares and products matrix for the full data set.
`H`	the Effect description sums of squares and products matrix (defined in the same way as the `mat` matrix).
`r`	the Expected rank of the H matrix. See the `Details` below.
`indices`	a numerical vector, matrix or 3-d array of integers giving the indices of the variables in the subset. If a matrix is specified, each row is taken to represent a different k-variable subset. If a 3-d array is given, it is assumed that the third dimension corresponds to different cardinalities.
`tolval`	the tolerance level to be used in checks for ill-conditioning and positive-definiteness of the 'total' and 'effects' (H) matrices. Values smaller than `tolval` are considered equivalent to zero.
`tolsym`	the tolerance level for symmetry of the covariance/correlation/total matrix and for the effects (`H`) matrix. If corresponding matrix entries differ by more than this value, the input matrices will be considered asymmetric and execution will be aborted. If corresponding entries are different, but by less than this value, the input matrix will be replaced by its symmetric part, i.e., input matrix A becomes (A+t(A))/2.

Details

Different kinds of statistical methodologies are considered within the framework, of a multivariate linear model:

X = A B + U

where X is the (nxp) data matrix of original variables, A is a known (nxp) design matrix, B an (qxp) matrix of unknown parameters and U an (nxp) matrix of residual vectors. The Xi squared index is related to the traditional test statistic (Bartllet-Pillai trace) and measures the contribution of each subset to an Effect characterized by the violation of a linear hypothesis of the form C B = 0, where C is a known cofficient matrix of rank r. The Bartllet-Pillai trace (P) is given by: P=tr(HT^{-1}) where H is the Effect matrix and T is the Total matrix. The Xi squared index is related to Bartllet-Pillai trace (P) by:

ξ2 = P/r

where r is the rank of H matrix.

The fact that indices can be a matrix or 3-d array allows for the computation of the Xi squared values of subsets produced by the search functions anneal, genetic, improve and eleaps (whose output option $subsets are matrices or 3-d arrays), using a different criterion (see the example below).

Value

The value of the xi2 coefficient.

Examples

## ---------------------------------------------------------------

## 1) A Linear Discriminant Analysis example with a very small data set. 
## We considered the Iris data and three groups, 
## defined by species (setosa, versicolor and virginica).

data(iris)
irisHmat <- ldaHmat(iris[1:4],iris$Species)
xi2.coef(irisHmat$mat,H=irisHmat$H,r=2,c(1,3))
## [1] 0.4942503

## ---------------------------------------------------------------

## 2) An example computing the value of the xi_2 criterion for two subsets 
## produced when the anneal function attempted to optimize the tau_2 
## criterion (using an absurdly small number of iterations).

tauresults<-anneal(irisHmat$mat,2,nsol=2,niter=2,criterion="tau2",
H=irisHmat$H,r=2)
xi2.coef(irisHmat$mat,H=irisHmat$H,r=2,tauresults$subsets)

##              Card.2
##Solution 1 0.5718811
##Solution 2 0.5232262

## ---------------------------------------------------------------

subselect

Selecting Variable Subsets

v0.15.2

GPL (>= 2)

Authors

Jorge Orestes Cerdeira [aut], Pedro Duarte Silva [aut], Jorge Cadima [aut, cre], Manuel Minhoto [aut]

Initial release

2020-03-04

xi2