I am using mixedpower following the notebook of Kumle et al here
My code is as follows:
powerfull1 <- mixedpower(model=full1, data=alldum, fixed_effects= c("core.non"), simvar="ID", steps=c(250, 275, 300, 325, 350), critical_value=2, n_sim = 1000)
The simulations stop after the first step, with the following error message:
[1] "Estimating power for step:"
[1] 250
[1] "Simulations for step 250 are based on 714.285714285714 successful single runs"
Error in [<-.data.frame
(*tmp*
, column_name, value = c(0.069, 0.027, :
replacement has 5 rows, data has 7
I have explored the solutions to this kind of error here.
However it's not clear where in mixedpower a dataframe is created, so I can't see how to apply the solution.
Any help? Thanks!