A script demonstrating the use of power analysis with R