Pwr function

Calculates power based on a model fit