Example SEM implementation for McDonald's ω coefficient

[mihaiconstantin]

# Install the `remotes` package if not available.
install.packages("remotes")

# Install the `powerly` package containing the `SEM` example.
remotes::install_github("mihaiconstantin/powerly", ref = "example/sem-omega")

# Load the package.
library(powerly)

# Run the method for the example.
results <- powerly(
    #  The lower and upper bounds of the sample size search range.
    range_lower = 150,
    range_upper = 450,

    # How many sample sizes to be in the candidate grid.
    samples = 50,

    # The number of Monte Carlo replications for each sample in the candidate grid.
    replications = 30,

    # This is the explicit estimand implementation name (i.e., redundant in new API).
    measure = "omega",

    # This is what we call reliability in version `2.0.0`.
    statistic = "power",
    
    # Check if the research question was answered within ε = 0.05 accuracy.
    measure_value = 1 - 0.05,

    # This is the value for the reliability (i.e., not McDonald's ω).
    statistic_value = .8,

    # This is the `SemOmega` model implemented for this example.
    model = "sem-omega",

    # Number of cores for the parallel processing.
    cores = 9,

    # Verbose output.
    verbose = TRUE
)

# Plot the sample size recommendation.
plot(results, step = 3)

# Validate.
validation <- validate(results, replications = 1500, cores = 9)

# Plot the validation.
plot(validation)