power_analysis.R