evSSCM2evShape transforms the eigenvalues of the SSCM of an elliptical distribution into that of the shape matrix.
(required) p-dimensional numeric representing the eigenvalues of the SSCM.
(optional) numeric, defines the stopping rule of the approximation procedure, see details.
(optional) numeric, defines the maximal number of iterations, see details.
The eigenvalues of the SSCM given that of the shape matrix can be calculated by evaluations of numerical integrals, see the help of
evShape2evSSCM or Dürre, Tyler, Vogel (2016). There is no expression for the inverse relationshop known. Though one can apply a fixed point iteration to get an approximation of the eigenvalues of the shape matrix. The iteration stops if either the maximal number of iterations is reached, which produces a warning, or if the error between the eigenvalues of the SSCM and the ones calculated from the actual fixed point iteration in L1 norm is smaller than the given tolerance. Since the mapping between the sets of eigenvalues is injective, see Dürre, Tyler, Vogel (2016), this gives a reasonable approximation of the eigenvalues of the shape matrix.
p-dimensional numerical, representing the eigenvalues of the shape matrix. They are standardized to sum to 1.
Dürre, A., Vogel, D., Fried, R. (2015): Spatial sign correlation, Journal of Multivariate Analyis, vol. 135, 89–105. arvix 1403.7635
Dürre, A., Tyler, D. E., Vogel, D. (2016): On the eigenvalues of the spatial sign covariance matrix in more than two dimensions, to appear in: Statistics and Probability Letters. arvix 1512.02863
Calculating the theoretical shape from the theoretical SSCM
Calculating the eigenvalues of the SSCM from the eigenvalues of the shape matrix
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16
# defining eigenvalues of the shape matrix evShape <- seq(from=0,to=1,by=0.1) # standardized to have sum 1 evShape <- evShape/sum(evShape) # calculating the related eigenvalues of the SSCM evSSCM <- evShape2evSSCM(evShape) plot(evShape,evSSCM) # recalculate the eigenvalues of the shape matrix evShape2 <- evSSCM2evShape(evSSCM) # error is negligible sum(abs(evShape-evShape2))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.