Description Usage Arguments Details Value Author(s) References See Also Examples

Based on a partial Newton-Raphson approach offers this function two ways to compute the symmetrized M-estimator of scatter. The user can choose if all pairwise differences are choosen and stored in the memory or if the computation and storage of this large matrix is to be avoided.

1 | ```
MVTMLEsymm(X, nu = 1, nmax = 500, eps = 1e-06, maxiter = 100, perm = FALSE)
``` |

`X` |
numeric data matrix or dataframe with more rows than columns. Missing values are not allowed. |

`nu` |
assumed degrees of freedom of the t-distribution, must be larger than 0. Default is '1'. |

`nmax` |
integer, if the sample size n (number of rows of |

`eps` |
convergence tolerance, which means that the algorithm stops when the Frobenius norm of the gradient is smaller than eps. |

`maxiter` |
maximum number of iterations. |

`perm` |
logical. If TRUE the rows of |

To get a good initial value for the algorithm, the estimator is first computed based on the pairwise differences of
successive observations. Therefore the order of the rows of `X`

is supposed to be random. If this is not the case, the data
should be first permuted using the argument `perm`

.

In case `maxiter`

is reached before convergence, the estimate at that iteration is returned and a warning is given.

A list containing:

`Sigma` |
Estimated scatter matrix. |

`iter` |
Number of iterations of the algorithm. |

Lutz Duembgen and Klaus Nordhausen

Duembgen, L., Nordhausen, K. and Schuhmacher, H. (2016), New algorithms for M-estimation of multivariate location and scatter, *Journal of Multivariate Analysis*, **144**, 200–217. doi: 10.1016/j.jmva.2015.11.009

1 2 | ```
MVTMLEsymm(longley)
MVTMLEsymm(longley, nmax=10)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.