MMboot_loccov: Fast and Robust Bootstrap for MM-estimates of Location and...

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

View source: R/MMboot_loccov.R

Description

Calculates bootstrapped MM-estimates of multivariate location and scatter using the Fast and Robust Bootstrap method.

Usage

1
MMboot_loccov(Y, R = 999, ests = MMest_loccov(Y))

Arguments

Y

matrix or data frame.

R

number of bootstrap samples. Default is R=999.

ests

original MM-estimates as returned by MMest_loccov().

Details

This function is called by FRBpcaMM and FRBhotellingMM, it is typically not to be used on its own. It requires the MM-estimates of multivariate location and scatter/shape (the result of MMest_loccov applied on Y), supplied through the argument ests. If ests is not provided, MMest_loccov calls the implementation of the multivariate MM-estimates in package rrcov of Todorov and Filzmoser (2009) with default arguments.

For multivariate data the fast and robust bootstrap was developed by Salibian-Barrera, Van Aelst and Willems (2006).

The value centered gives a matrix with R columns and 2*(p+p*p) rows (p is the number of variables in Y), containing the recalculated estimates of the MM-location, MM-shape, S-covariance and S-location. Each column represents a different bootstrap sample. The first p rows are the MM-location estimates, the next p*p rows are the MM-shape estimates (vectorized). Then the next p*p rows are the S-covariance estimates (vectorized) and the final p rows are the S-location estimates. The estimates are centered by the original estimates, which are also returned through MMest in vectorized form.

Value

A list containing:

centered

recalculated MM- and S-estimates of location and scatter (centered by original estimates), see Details

MMest

original MM- and S-estimates of location and scatter, see Details

Author(s)

Gert Willems, Ella Roelant and Stefan Van Aelst

References

See Also

FRBpcaMM, FRBhotellingMM, Sboot_loccov

Examples

1
2
3
Y <- matrix(rnorm(50*5), ncol=5)
MMests <- MMest_loccov(Y) 
bootresult <- MMboot_loccov(Y, R = 1000, ests = MMests)

FRB documentation built on May 29, 2017, 5:45 p.m.

Related to MMboot_loccov in FRB...