tyler.shape {ICSNP} | R Documentation |
Iterative algorithm to estimate Tyler's shape matrix.
tyler.shape(X, location = NULL, init = NULL, steps = Inf, eps = 1e-06, maxiter = 100, in.R = FALSE, print.it = FALSE, na.action = na.fail)
X |
numeric data matrix or dataframe. |
location |
if NULL the sample mean is used, otherwise a vector with the location can be specified. |
init |
an optional matrix giving the starting value for the iteration |
steps |
a fixed number of iteration steps to take. See details. |
eps |
convergence tolerance. |
maxiter |
maximum number of iterations. |
in.R |
logical. If TRUE R-code (and not C) is used in the iteration |
print.it |
logical. If TRUE prints the number of iterations, otherwise not. |
na.action |
a function which indicates what should happen when the data contain 'NA's. Default is to fail. |
The most robust M-estimator of shape. It is proportional to the regular covariance matrix for elliptical contoured distributions. The estimate is in such a way standardized, that its determinate is 1.
The algorithm requires an estimate of location, if none is provided, the sample mean is used. Observations which are equal to the location estimate are removed form the data.
The algorithm also allows for a k-step version where the iteration is run for a fixed number of steps instead of until convergence. If steps
is finite that number of steps is taken and maxiter
is ignored.
A matrix.
Klaus Nordhausen, klaus.nordhausen@uta.fi, and Seija Sirkiä, ssirkia@maths.jyu.fi
Tyler, D.E. (1987), A distribution-free M-estimator of scatter, Annals of Statistics, 15, 234–251.
set.seed(654321) cov.matrix <- matrix(c(3,2,1,2,4,-0.5,1,-0.5,2), ncol=3) X <- rmvnorm(100, c(0,0,0), cov.matrix) tyler.shape(X) tyler.shape(X, location=0) cov.matrix/det(cov.matrix)^(1/3) rm(.Random.seed)