pdf("ex0303.pdf")


# Power Calculations for Example 3.3
#  Note that what r calls delta is mu1-mu2, not (mu1-mu2)/sigma
#  sd is sigma, sig.level is alpha, power is 1-beta, alt is alternative (2-sided here)

power.t.test(delta=15,sd=34.4,sig.level=0.05,power=0.8,alt="two.sided")



dev.off()