tests/testthat/test-gorica_and_sem.R

In cjvanlissa/gorica: Evaluation of Inequality Constrained Hypotheses Using GORICA

if(FALSE){


library(gorica)
library(lavaan)

#  Specify and fit the confirmatory factor model
model1 <- '
    A =~ Ab + Al + Af + An + Ar + Ac
    B =~ Bb + Bl + Bf + Bn + Br + Bc
'

# Use the lavaan sem function to execute the confirmatory factor analysis
fit1 <- sem(model1, data = sesamesim, std.lv = TRUE)

# Inspect the parameter names (such that hypotheses can be specified using these)
coef(fit1)

# Formulate hypotheses
hypotheses1 <- "
A=~Ab > .6 & A=~Al > .6 & A=~Af > .6 & A=~An > .6 & A=~Ar > .6 & A=~Ac > .6 &
B=~Bb > .6 & B=~Bl > .6 & B=~Bf > .6 & B=~Bn > .6 & B=~Br > .6 & B=~Bc > .6
"

# Call gorica
set.seed(100)
out_c <- gorica(fit1, hypotheses1, comparison = "complement", standardize=TRUE) # TO DO doet laatste iets??
out_c
# TO DO Krijg nu een negatieve PT voor het complemet (zou iets tussen 0 en 12 moeten zijn, ws dichtbij 12)
# NB In restriktor gaat dit wel goed - zie regressie voorbeeld hieronder.
out_u <- gorica(fit1, hypotheses1, standardize=TRUE) # TO DO doet laatste iets??
out_u


}