tests/testthat/test-obs.resid.R

In philchalmers/faoutlier: Influential Case Detection Methods for Factor Analysis and Structural Equation Models

context('obs.resid')

test_that('obs.resid run', {

    #Exploratory
    nfact <- 3
    ORresult <- obs.resid(holzinger, nfact)
    ORresult.outlier <- obs.resid(holzinger.outlier, nfact)
    expect_equal(ORresult$std_res[1:3], c(-0.01285074, -0.75147106, -1.42004924),
                 tolerance = 1e-5)
    expect_is(ORresult, 'obs.resid')
    expect_is(ORresult.outlier, 'obs.resid')
    expect_is(plot(ORresult), 'trellis')
    expect_is(plot(ORresult.outlier), 'trellis')

    #-------------------------------------------------------------------
    suppressMessages(model <- sem::specifyModel(file='sem-model/sem-model.txt', quiet=TRUE))
    ORresult <- suppressWarnings(obs.resid(holzinger, model))
    ORresult.outlier <- suppressWarnings(obs.resid(holzinger.outlier, model))
    expect_equal(ORresult$std_res[1:3], c(0.2548177, -0.5300287, -1.8518586),
                 tolerance = 1e-5)
    expect_is(ORresult, 'obs.resid')
    expect_is(ORresult.outlier, 'obs.resid')
    expect_is(plot(ORresult), 'trellis')
    expect_is(plot(ORresult.outlier), 'trellis')

    #-------------------------------------------------------------------
    model <- 'F1 =~  Remndrs + SntComp + WrdMean
    F2 =~ MissNum + MxdArit + OddWrds
    F3 =~ Boots + Gloves + Hatchts'

    obs.resid2 <- obs.resid(holzinger, model, orthogonal=TRUE)
    obs.resid2.outlier <- obs.resid(holzinger.outlier, model, orthogonal=TRUE)
    expect_equal(ORresult$std_res[1:3], c(0.2548177, -0.5300287, -1.8518586),
                 tolerance = 1e-5)
    expect_is(obs.resid2, "obs.resid")
    expect_is(obs.resid2.outlier, "obs.resid")
    expect_is(plot(obs.resid2), 'trellis')
    expect_is(plot(obs.resid2.outlier), 'trellis')
})