neg.rmsea: Equivalence Tests for RMSEA

View source: R/neg.rmsea.R

neg.rmseaR Documentation

Equivalence Tests for RMSEA

Description

Function performs one of four equivalence tests for RMSEA fit index.

Usage

neg.rmsea(
  mod,
  alpha = 0.05,
  eq.bound,
  modif.eq.bound = FALSE,
  ci.method = "not.close",
  nbootpd = 50L,
  nboot = 250L,
  round = 3,
  plot = TRUE,
  saveplot = FALSE
)

## S3 method for class 'neg.rmsea'
print(x, ...)

Arguments

mod

lavaan model object

alpha

desired alpha level (default = .05)

eq.bound

upper end of equivalence interval for comparison; must be .01, .05, .08 or .10 if modif.eq.bound = TRUE

modif.eq.bound

should the upper end of the equivalence interval be modified? (default = FALSE)

ci.method

method used to calculate confidence interval; options are "not.close" or "yhy.boot"; "not.close" corresponds to (1-2alpha) percent CI, "yhy.boot" corresponds to (1-2alpha) percent boot CI (default = "not.close")

nbootpd

number of bootstrap samples by "yhy.boot" for pd function

nboot

number of bootstrap samples if "yhy.boot" is selected as ci.method (default = 250L)

round

number of digits to round equivalence bound and confidence interval bounds (default = 3)

plot

logical, plotting the results (default = TRUE)

saveplot

saving plots (default = FALSE)

x

object of class neg.rmsea

...

extra arguments

Details

The user specifies the lavaan fitted model object, the desired equivalence bound, and method of confidence interval computation. By default, the function does not modify the equivalence bounds according to Yuan et al. (2016). The user can also choose to instead run an equivalence test using a modified equivalence bound if the equivalence bound to be modified is .01, .05, .08, or .10. Alpha level can also be modified.

For information on modified equivalence bounds see Yuan, K. H., Chan, W., Marcoulides, G. A., & Bentler, P. M. (2016). Assessing structural equation models by equivalence testing with adjusted fit indexes. Structural Equation Modeling: A Multidisciplinary Journal, 23(3), 319-330. doi: https://doi.org/10.1080/10705511.2015.1065414.

The proportional distance quantifies the proportional distance from 0 to the nearest negligible effect (equivalence) interval (here, eiU). As values get farther from 0 the relationship becomes more substantial, with values greater than 1 indicating that the effect falls outside of the negligible effect (equivalence) interval.

Value

returns a list including the following:

  • title1 The title of the RMSEA equivalence test. The appropriate title of the test will be displayed depending on the ci.method chosen and whether modif.eq.bound is TRUE or FALSE.

  • rmsea_index The RMSEA index.

  • ci.method The method for confidence interval calculation (direct computation or bootstrap).

  • rmsea_eq The upper end of the 1-2*alpha confidence interval for the RMSEA index.

  • eq.bound The equivalence bound.

  • PD Proportional distance (PD).

  • cilpd Lower bound of the 1-alpha CI for the PD.

  • ciupd Upper bound of the 1-alpha CI for the PD.

Author(s)

Rob Cribbie cribbie@yorku.ca and Nataly Beribisky natalyb1@yorku.ca

@export neg.rmsea

Examples

d <- lavaan::HolzingerSwineford1939
hs.mod <- 'visual =~ x1 + x2 + x3
textual =~ x4 + x5 + x6
speed =~ x7 + x8 + x9'
fit1 <- lavaan::cfa(hs.mod, data = d)
neg.rmsea(alpha = .05, mod = fit1, eq.bound = .05, ci.method = "not.close", modif.eq.bound = FALSE,
round = 5, nboot = 50)

negligible documentation built on Sept. 11, 2024, 9:24 p.m.