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')
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
