Learn R Programming

cSEM (version 0.5.0)

testOMF: Test for overall model fit

Description

lifecycle::badge("maturing")

Usage

testOMF(
 .object                = NULL, 
 .alpha                 = 0.05,
 .fit_measures          = FALSE,
 .handle_inadmissibles  = c("drop", "ignore", "replace"), 
 .R                     = 499, 
 .saturated             = FALSE,
 .seed                  = NULL,
 ...
)

Value

A list of class cSEMTestOMF containing the following list elements:

$Test_statistic

The value of the test statistics.

$Critical_value

The corresponding critical values obtained by the bootstrap.

$Decision

The test decision. One of: FALSE (Reject) or TRUE (Do not reject).

$Information

The .R bootstrap values; The number of admissible results; The seed used and the number of total runs.

Arguments

.object

An R object of class cSEMResults resulting from a call to csem().

.alpha

An integer or a numeric vector of significance levels. Defaults to 0.05.

.fit_measures

Logical. (EXPERIMENTAL) Should additional fit measures be included? Defaults to FALSE.

.handle_inadmissibles

Character string. How should inadmissible results be treated? One of "drop", "ignore", or "replace". If "drop", all replications/resamples yielding an inadmissible result will be dropped (i.e. the number of results returned will potentially be less than .R). For "ignore" all results are returned even if all or some of the replications yielded inadmissible results (i.e. number of results returned is equal to .R). For "replace" resampling continues until there are exactly .R admissible solutions. Depending on the frequency of inadmissible solutions this may significantly increase computing time. Defaults to "drop".

.R

Integer. The number of bootstrap replications. Defaults to 499.

.saturated

Logical. Should a saturated structural model be used? Defaults to FALSE.

.seed

Integer or NULL. The random seed to use. Defaults to NULL in which case an arbitrary seed is chosen. Note that the scope of the seed is limited to the body of the function it is used in. Hence, the global seed will not be altered!

...

Can be used to determine the fitting function used in the calculateGFI function.

Details

Bootstrap-based test for overall model fit originally proposed by Beran1985;textualcSEM. See also Dijkstra2015;textualcSEM who first suggested the test in the context of PLS-PM.

By default, testOMF() tests the null hypothesis that the population indicator correlation matrix equals the population model-implied indicator correlation matrix. Several discrepancy measures may be used. By default, testOMF() uses four distance measures to assess the distance between the sample indicator correlation matrix and the estimated model-implied indicator correlation matrix, namely the geodesic distance, the squared Euclidean distance, the standardized root mean square residual (SRMR), and the distance based on the maximum likelihood fit function. The reference distribution for each test statistic is obtained by the bootstrap as proposed by Beran1985;textualcSEM.

It is possible to perform the bootstrap-based test using fit measures such as the CFI, RMSEA or the GFI if .fit_measures = TRUE. This is experimental. To the best of our knowledge the applicability and usefulness of the fit measures for model fit assessment have not been formally (statistically) assessed yet. Theoretically, the logic of the test applies to these fit indices as well. Hence, their applicability is theoretically justified. Only use if you know what you are doing.

If .saturated = TRUE the original structural model is ignored and replaced by a saturated model, i.e., a model in which all constructs are allowed to correlate freely. This is useful to test misspecification of the measurement model in isolation.

References

See Also

csem(), calculateSRMR(), calculateDG(), calculateDL(), cSEMResults, testMICOM(), testMGD(), exportToExcel()

Examples

Run this code
# ===========================================================================
# Basic usage
# ===========================================================================
model <- "
# Structural model
eta2 ~ eta1
eta3 ~ eta1 + eta2

# (Reflective) measurement model
eta1 =~ y11 + y12 + y13
eta2 =~ y21 + y22 + y23
eta3 =~ y31 + y32 + y33
"

## Estimate
out <- csem(threecommonfactors, model, .approach_weights = "PLS-PM")

## Test
testOMF(out, .R = 50, .seed = 320)

Run the code above in your browser using DataLab