We're conducting a power analysis using wp.poisson to get the relative risk for different levels of a normally distributed predict (mean = 0, sd = 1). Our dependent variable is binary, but it was suggested that we use a poisson rather than logistic since the prevalence of the behavior we're observing is over 10% (22-40% has been observed in previous studies). The observed intercept in simulated data is -1.91, based on the 22% occurrence, so the IRR is 0.148. I used the relative risk for a small, medium and large effect as 1.22, 1.86, and 3.00, as suggested by Olivier et al (2017). I get an error when I try to run the medium and large effects of "non-finite function value." Here's my code:
medium <- wp.poisson(n = seq(40, 100, 10), exp0 = 0.15, exp1 = 1.86, alpha = 0.05, power = NULL, family = "normal", parameter = c(0, 1))
Any suggestions for trouble shooting?