[R-br] filliben
Elias T Krainski
eliaskrainski em yahoo.com.br
Sexta Dezembro 6 17:06:43 BRST 2013
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')
Mais detalhes sobre a lista de discussão R-br