bootDominanceAnalysis {dominanceanalysis}

Bootstrap Analysis for Dominance Analysis

Description

Implements a bootstrap procedure as presented by Azen and Budescu (2003). Provides the expected level of dominance of predictor X_i over X_j, as the degree to which the pattern found in the sample is reproduced in the bootstrap samples.

Usage

bootDominanceAnalysis(
  x,
  R,
  constants = c(),
  terms = NULL,
  fit.functions = "default",
  null.model = NULL,
  ...
)

Arguments

x

An object of class lm, glm, or lmer.

R

The number of bootstrap resamples.

constants

A vector of predictors to remain unchanged between models, i.e., variables not subjected to bootstrap analysis.

terms

A vector of terms to be analyzed. By default, terms are obtained from the model.

fit.functions

A list of functions providing fit indices for the model. Refer to fit.functions parameter in dominanceAnalysis function.

null.model

Applicable only for linear mixed models. It refers to the null model against which to test the submodels, i.e., only random effects, without any fixed effects.

...

Additional arguments provided to lm or lmer (not implemented yet).

Details

Use summary() to obtain a nicely formatted data.frame.

Value

An object of class bootDominanceAnalysis containing:

boot

The results of the bootstrap analysis.

preds

The predictors analyzed.

fit.functions

The fit functions used in the analysis.

c.names

A vector where each value represents the name of a specific dominance analysis result. Names are prefixed with the type of dominance (complete, conditional, or general), and the fit function used, followed by the names of the first and second predictors involved in the comparison.

m.names

Names of each one the predictor pairs.

terms

The terms analyzed.

R

The number of bootstrap resamples.

Examples

lm.1 <- lm(Employed ~ ., longley)
da.boot <- bootDominanceAnalysis(lm.1, R = 1000)
summary(da.boot)