Learn R Programming

robustHD (version 0.5.0)

tslars: (Robust) least angle regression for time series data

Description

(Robustly) sequence groups of candidate predictors and their respective lagged values according to their predictive content and find the optimal model along the sequence. Note that lagged values of the response are included as a predictor group as well.

Usage

tslars(x, ...)

## S3 method for class 'formula': tslars(formula, data, ...)

## S3 method for class 'default': tslars(x, y, h = 1, pMax = 3, sMax = NA, fit = TRUE, s = c(0, sMax), crit = "BIC", ncores = 1, cl = NULL, model = TRUE, ...)

rtslars(x, ...)

## S3 method for class 'formula': rtslars(formula, data, ...)

## S3 method for class 'default': rtslars(x, y, h = 1, pMax = 3, sMax = NA, centerFun = median, scaleFun = mad, regFun = lmrob, regArgs = list(), combine = c("min", "euclidean", "mahalanobis"), winsorize = FALSE, const = 2, prob = 0.95, fit = TRUE, s = c(0, sMax), crit = "BIC", ncores = 1, cl = NULL, seed = NULL, model = TRUE, ...)

Arguments

formula
a formula describing the full model.
data
an optional data frame, list or environment (or object coercible to a data frame by as.data.frame) containing the variables in the model. If not found in data, the variables are taken from
x
a numeric matrix or data frame containing the candidate predictor series.
y
a numeric vector containing the response series.
h
an integer giving the forecast horizon (defaults to 1).
pMax
an integer giving the maximum number of lags in the model (defaults to 3).
sMax
an integer giving the number of predictor series to be sequenced. If it is NA (the default), predictor groups are sequenced as long as there are twice as many observations as predictor variables.
centerFun
a function to compute a robust estimate for the center (defaults to median).
scaleFun
a function to compute a robust estimate for the scale (defaults to mad).
regFun
a function to compute robust linear regressions that can be interpreted as weighted least squares (defaults to lmrob).
regArgs
a list of arguments to be passed to regFun.
combine
a character string specifying how to combine the data cleaning weights from the robust regressions with each predictor group. Possible values are "min" for taking the minimum weight for each observation, "euclidean" for w
winsorize
a logical indicating whether to clean the data by multivariate winsorization.
const
numeric; tuning constant for multivariate winsorization to be used in the initial corralation estimates based on adjusted univariate winsorization (defaults to 2).
prob
numeric; probability for the quantile of the $\chi^{2}$ distribution to be used in multivariate winsorization (defaults to 0.95).
fit
a logical indicating whether to fit submodels along the sequence (TRUE, the default) or to simply return the sequence (FALSE).
s
an integer vector of length two giving the first and last step along the sequence for which to compute submodels. The default is to start with a model containing only an intercept (step 0) and iteratively add all series along the sequence (step <
crit
a character string specifying the optimality criterion to be used for selecting the final model. Currently, only "BIC" for the Bayes information criterion is implemented.
ncores
a positive integer giving the number of processor cores to be used for parallel computing (the default is 1 for no parallelization). If this is set to NA, all available processor cores are used. For each lag length, parallel computi
cl
a parallel cluster for parallel computing as generated by makeCluster. This is preferred over ncores for tasks that are parallelized on the Rlevel, in which case
seed
optional initial seed for the random number generator (see .Random.seed), which is useful because many robust regression functions (including lmrob

Value

References

Gelper, S. and Croux, C. (2010) Time series least angle regression for selecting predictive economic sentiment series. Working paper.

See Also

coef, fitted, plot, predict, residuals, tslarsP, lmrob