################################################### ### chunk number 1: ################################################### n = 15 pop.sd = 21 SD = pop.sd/sqrt(n) SD mu = 75 region=qnorm(c(.025,.975),mu,SD) region ################################################### ### chunk number 2: ################################################### mu = 80 pnorm(region[1],mu,SD) pnorm(region[2],mu,SD,lower.tail=FALSE) ################################################### ### chunk number 3: ################################################### mus = seq(55,95,by=5) cbind(mu=mus,power=pnorm(region[1],mus,SD)+ pnorm(region[2],mus,SD,lower.tail=FALSE)) ################################################### ### chunk number 4: ################################################### mus = seq(50,100,by=1) plot(mus,pnorm(region[1],mus,SD)+ pnorm(region[2],mus,SD,lower.tail=FALSE),type="l", xlab="Population mean (mu)",ylab="Power") ################################################### ### chunk number 5: ################################################### n.draws = 200 mu = 40 pop.sd = 10 n = 18 draws = matrix( rnorm( n.draws * n, mu, pop.sd ), 18 ) get.p.value = function( x ) t.test( x, mu = 40 )$p.value pvalues = apply( draws, 2, get.p.value ) hist(pvalues,breaks=20,prob=TRUE) abline(h=1,lwd=2,lty=2) ################################################### ### chunk number 6: ################################################### mu = 37 draws = matrix( rnorm( n.draws * n, mu, pop.sd ), 18 ) pvalues = apply( draws, 2, get.p.value ) hist(pvalues,breaks=20,prob=TRUE) abline(h=1,lwd=2,lty=2)