#==============================================================================
# Determinando o menor tamanho de amostra
# Exemplo pag. 6-8 da aula_09.pdf
#==============================================================================
n   <- 2:30
p1  <- 2*pnorm(sqrt(n)/5)-1
p2  <- pchisq(1.44*n,n-1)-pchisq(0.64*n,n-1)

out <- cbind(n,p1,p2,p1*p2)
colnames(out) <- c("n","p1","p2","p1*p2")
print(out)

       n        p1        p2      p1*p2
 [1,]  2 0.2227026 0.1682130 0.03746147
 [2,]  3 0.2709655 0.2675678 0.07250164
 [3,]  4 0.3108435 0.3406564 0.10589083
 [4,]  5 0.3452792 0.3992418 0.13784988
 [5,]  6 0.3757939 0.4483605 0.16849112
 [6,]  7 0.4032988 0.4906854 0.19789284
 [7,]  8 0.4283924 0.5278313 0.22611891
 [8,]  9 0.4514938 0.5608617 0.25322554
 [9,] 10 0.4729107 0.5905216 0.27926402
[10,] 11 0.4928775 0.6173583 0.30428204
[11,] 12 0.5115777 0.6417880 0.32832441
[12,] 13 0.5291583 0.6641368 0.35143349
[13,] 14 0.5457398 0.6846659 0.37364942
[14,] 15 0.5614220 0.7035891 0.39501036
[15,] 16 0.5762892 0.7210835 0.41555264
[16,] 17 0.5904133 0.7372985 0.43531085
[17,] 18 0.6038561 0.7523613 0.45431798
[18,] 19 0.6166715 0.7663813 0.47260550
[19,] 20 0.6289066 0.7794535 0.49020345
[20,] 21 0.6406032 0.7916608 0.50714049
[21,] 22 0.6517983 0.8030766 0.52344400
[22,] 23 0.6625250 0.8137657 0.53914014
[23,] 24 0.6728131 0.8237858 0.55425387
[24,] 25 0.6826895 0.8331886 0.56880908
[25,] 26 0.6921785 0.8420206 0.58282858
[26,] 27 0.7013024 0.8503238 0.59633418
[27,] 28 0.7100815 0.8581362 0.60934671
[28,] 29 0.7185345 0.8654924 0.62188612
[29,] 30 0.7266783 0.8724238 0.63397146