davharris · February 13, 2015 21:21
diff --git a/low_power.R b/low_power.R
 n = 10
 delta = 0.1
 sd = 1

 # Calculate power with n and delta
 power.t.test(n = n, delta = delta, sd = 1, type = "one.sample") # ~5% power

 # Run 10,000 t.tests
 p = replicate(1E4, t.test(rnorm(n, mean = delta, sd = 1))$p.value)

 mean(p < .025 | p > .975) # reject ~5.2% of the time
	n = 10
	delta = 0.1
	sd = 1

	# Calculate power with n and delta
	power.t.test(n = n, delta = delta, sd = 1, type = "one.sample") # ~5% power

	# Run 10,000 t.tests
	p = replicate(1E4, t.test(rnorm(n, mean = delta, sd = 1))$p.value)

	mean(p < .025 \| p > .975) # reject ~5.2% of the time