I'll throw a few ideas out there. After looking at your example syntax, I briefly played around with the power specification in line 203 [configural <- semPower.aPriori(effect = c(.08, .095), effect.measure = ‘RMSEA’, alpha = .05, power = .90, df = c(1347, 1365))]. When I enter “power = .10” there, it still returns n = 10. Maybe something is missing in the syntax?
So far, I’ve only worked with Monte Carlo simulations, using the great tool by Wang and Rhemtulla:
Wang, Y. Andre, and Mijke Rhemtulla. 2021. “Power Analysis for Parameter Estimation in Structural Equation Modeling: A Discussion and Tutorial.” Advances in Methods and Practices in Psychological Science 4 (1): 2515245920918253.
https://doi.org/10.1177/2515245920918253 This allows the power to be calculated separately for each parameter.
local power (power of determining specific effects within the model). See also here:
https://jihongzhang.org/posts/2022-04-29-power-analysis-for-sem/This has worked best for me so far.
Best regards, Marcus