[R-br] filliben

[Elias T Krainski]

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')