<font face="trebuchet ms,sans-serif">Alexandre,<br><br>Segue um CMR simples:<br><br><span style="font-family:courier new,monospace"># tamanho de amostra para medida de interesse que<br># é o erro absoluto ao estimar a média da população<br>
# erro = abs(média_estimada-média_verdadeira)<br><br>pop.total <- rnorm(10000, mean=20, sd=15)    # população<br>m <- mean(20)                                # média verdadeira<br>i <- unique(floor(2^seq(1, 12, by=0.5))); i  # tamanhos de amostra<br>
nr <- 100                                    # número de execuções de reamostragem<br>fun <- function(size, pop, m){<br>  # função que calcula o erro a partir de uma amostra<br>  x <- sample(pop, size, replace=TRUE)<br>
  # retorna a diferença entre estimado e real<br>  abs(mean(x)-m)<br>}<br><br>results1 <- sapply(i,<br>                   function(j){<br>                     replicate(nr, fun(size=j, pop=pop.total, m=m))<br>                   })<br>
str(results1) # matriz: colunas correspondem à i<br>results2 <- colMeans(results1) # tira a média para cada i<br>results2<br><br>plot(results2~i)<br>plot(results2~log2(i))<br><br># qual o tamanho de amostra para ter um erro no ponto<br>
# médio entre o máximo (quando n=1) e o mínimo (n->infinito, erro=0)<br><br>n0 <- nls(results2~A*B^(-log2(i)+1), start=list(A=1, B=1))<br>coef(n0)<br><br>perc <- 0.3                           # erro tolerado de 30% do máximo (usuário define)<br>
xi <- -log(perc)/log(coef(n0)["B"])+1 # tamanho de amostra correspondente<br><br>plot(results2~log2(i),<br>     xlim=c(1, max(log2(i))), ylim=c(0,max(results2, coef(n0)["A"])),<br>     xlab="Tamanho da amostra (log2 unidades)",<br>
     ylab="Erro absoluto médio ao estimar a média")<br>abline(v=1, lty=2)<br>with(as.list(coef(n0)), curve(A*B^(-x+1), add=TRUE))<br>with(as.list(coef(n0)), abline(h=A*c(1,perc), v=xi, lty=2, col=2))<br><br>ceiling(2^xi) # tamanho de amostra à ser considerado<br>
</span><br>À disposição.<br>Walmes.<br><br clear="all"></font><div><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a><br>
skype: walmeszeviani<br style="font-family:trebuchet ms,sans-serif"></span><span style="font-family:trebuchet ms,sans-serif">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>
</div><br>