I know how to compute power in equal sample size anova, but I have a difficulty when the data is unequal sample size......

This is how to compute power in equal sample size anova, the "n=" means it can only be used in equal sample size.

pwr.anova.test(k=,n=,f)

I've already compute f, and this is my code:

f<-sqrt(sum((means-mean(dta1$Y))^2/(3*10207)))

Do you have any solution to solve the "unequal" problem?

Thank you in advance. I'm very appreciate!