The purpose of this code is to create 500 simulated datasets and find the number of p-values lower than 0.05 (power) as well as the mean parameter estimate calculated with f