The within argument was ill defined.
power_curve2 <- powerCurve(model, within="g", breaks=c(30,40,50,100), nsim = 100, test = fixed("x", method = "z") ) plot(power_curve2)
gives you something like: