<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Perfeito Walmes,<br>
<br>
Era isto mesmo que eu buscava,<br>
<br>
Abraço,<br>
<br>
Alexandre<br>
<br>
<br>
<div class="moz-cite-prefix">Em 20/03/2013 17:31, Walmes Zeviani
escreveu:<br>
</div>
<blockquote
cite="mid:CAFU=Ekb9-T2UZC_v3n-2_aoz=fhu6ZwF06a2a2cjYsexJWm-2w@mail.gmail.com"
type="cite"><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
moz-do-not-send="true" 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
moz-do-not-send="true"
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>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
======================================================================
Alexandre dos Santos
Proteção Florestal
Coordenador do curso Técnico em Florestas
Vice Coordenador do curso de Engenharia Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
======================================================================
</pre>
</body>
</html>