Exemplo de comparação de eficiência de estimadores

## X ~ N(0,1) N=1000 amostras de tamanho 20
sim1 <- matrix(rnorm(20*1000), nc=20)
medias1 <- apply(sim1,1,mean)
medianas1 <- apply(sim1,1,median)
var(medias1)
var(medianas1)
rm(sim1)
 
##
## X ~ N(3.1, 0,4/3) N=5000, amostras de tamanho 50
sim2 <- matrix(rnorm(50*5000), nc=20)
medias2 <- apply(sim2,1,mean)
medianas2 <- apply(sim2,1,median)
var(medias2)
var(medianas2)