6 Dez
2013
6 Dez
'13
19:06
apenas para ilustrar o poder desse tipo de teste: filliben = function(x) { y = sort(x) n = length(y) p = (1:n-0.3175)/(n+0.365) p[1] = 1-0.5^(1/n) p[n] = 0.5^(1/n) p = qnorm(p) cor(y,p) } e = replicate(1000, filliben(rexp(100, 1))) u = replicate(1000, filliben(runif(100))) n = replicate(1000, filliben(rnorm(100))) hist(n, xlim=range(e,u,n)) hist(u, add=TRUE, border='blue') hist(e, add=TRUE, border='red')