Extract binary weights that indicate outliers from sparse least trimmed squares regression models.

1 2 3 4 5 |

`object` |
the model fit from which to extract outlier weights. |

`s` |
an integer vector giving the indices of the models for which to
extract outlier weights. If |

`fit` |
a character string specifying for which estimator to extract
outlier weights. Possible values are |

`drop` |
a logical indicating whether to reduce the dimension to a vector in case of only one model. |

`...` |
currently ignored. |

A numeric vector or matrix containing the requested outlier weights.

The weights are *1* for observations with reasonably small
residuals and *0* for observations with large residuals.

Andreas Alfons

`sparseLTS`

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | ```
## generate data
# example is not high-dimensional to keep computation time low
library("mvtnorm")
set.seed(1234) # for reproducibility
n <- 100 # number of observations
p <- 25 # number of variables
beta <- rep.int(c(1, 0), c(5, p-5)) # coefficients
sigma <- 0.5 # controls signal-to-noise ratio
epsilon <- 0.1 # contamination level
Sigma <- 0.5^t(sapply(1:p, function(i, j) abs(i-j), 1:p))
x <- rmvnorm(n, sigma=Sigma) # predictor matrix
e <- rnorm(n) # error terms
i <- 1:ceiling(epsilon*n) # observations to be contaminated
e[i] <- e[i] + 5 # vertical outliers
y <- c(x %*% beta + sigma * e) # response
x[i,] <- x[i,] + 5 # bad leverage points
## sparse LTS over a grid of values for lambda
# fit model
frac <- seq(0.2, 0.05, by = -0.05)
fitGrid <- sparseLTS(x, y, lambda = frac, mode = "fraction")
# extract outlier weights
wt(fitGrid)
head(wt(fitGrid, fit = "both"))
head(wt(fitGrid, s = NULL))
head(wt(fitGrid, fit = "both", s = NULL))
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.