Distribuição Empírica

domingo, 28 de outubro de 2007

Esse é um exemplo encontrado no help do software R.
Comando: help (ecdf)



CÓDIGOS PARA UTILIZAÇÃO NO SOFTWARE R
Fn <- ecdf(rnorm(12))
tt <- seq(-2,2, by = 0.1)
y <- round(rnorm(12),1); y[3] <- y[1]
Fn12 <- ecdf(y)
require(graphics)
op <- par(mfrow=c(3,1), mgp=c(1.5, 0.8,0), mar= .1+c(3,3,2,1))
F10 <- ecdf(rnorm(10))
plot(F10)
plot(F10, verticals= TRUE, do.points = FALSE)
plot(Fn12 , lwd = 2) ; mtext("lwd = 2", adj=1)

0 Observações: