View source: R/fit_TylerCauchy.R
fit_Cauchy  R Documentation 
Estimate parameters of a multivariate elliptical distribution, namely, the mean vector and the covariance matrix, to fit data. Any data sample with NAs will be simply dropped. The estimation is based on the maximum likelihood estimation (MLE) under a Cauchy distribution and the algorithm is obtained from the majorizationminimization (MM) optimization framework. The Cauchy distribution does not have secondorder moments and the algorithm actually estimates the scatter matrix. Nevertheless, assuming that the observed data has secondorder moments, the covariance matrix is returned by computing the missing scaling factor with a very effective method.
fit_Cauchy(
X,
initial = NULL,
max_iter = 200,
ptol = 0.001,
ftol = Inf,
return_iterates = FALSE,
verbose = FALSE
)
X 
Data matrix containing the multivariate time series (each column is one time series). 
initial 
List of initial values of the parameters for the iterative estimation method. Possible elements include:

max_iter 
Integer indicating the maximum number of iterations for the iterative estimation
method (default is 
ptol 
Positive number indicating the relative tolerance for the change of the variables
to determine convergence of the iterative method (default is 
ftol 
Positive number indicating the relative tolerance for the change of the loglikelihood
value to determine convergence of the iterative method (default is 
return_iterates 
Logical value indicating whether to record the values of the parameters (and possibly the
loglikelihood if 
verbose 
Logical value indicating whether to allow the function to print messages (default is 
A list containing possibly the following elements:
mu 
Mean vector estimate. 
cov 
Covariance matrix estimate. 
scatter 
Scatter matrix estimate. 
converged 
Boolean denoting whether the algorithm has converged ( 
num_iterations 
Number of iterations executed. 
cpu_time 
Elapsed CPU time. 
log_likelihood 
Value of loglikelihood after converge of the estimation algorithm (if 
iterates_record 
Iterates of the parameters ( 
Daniel P. Palomar
Ying Sun, Prabhu Babu, and Daniel P. Palomar, "Regularized Robust Estimation of Mean and Covariance Matrix Under HeavyTailed Distributions," IEEE Trans. on Signal Processing, vol. 63, no. 12, pp. 30963109, June 2015.
fit_Tyler
and fit_mvt
library(mvtnorm) # to generate heavytailed data
library(fitHeavyTail)
X < rmvt(n = 1000, df = 6) # generate Student's t data
fit_Cauchy(X)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.