Caros Alexandre e Walmes.<div><br></div><div>Estou usando este teste de decisão de número de amostras para ver se minha amostragem está ok.</div><div><br></div><div>Mas fiquei na dúvida sobre o valor 10000 da linha</div><div>
<div class="im" style="color:rgb(80,0,80);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)"><blockquote type="cite" style="border-left-width:2px;border-left-style:solid;border-left-color:rgb(204,204,204);padding:0px 1em">
<font face="trebuchet ms,sans-serif"><span style="font-family:'courier new',monospace">pop.total <- rnorm(10000, mean=20, sd=15)<br><br></span></font></blockquote><div><br></div><div><br></div></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
Vou exemplificar com dados de previsão de safra em café, no caso, são relativos a altura da planta de café.</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
Estas medidas referem a coleta em 3 pontos da lavoura, cada ponto é composto por 5 plantas, então sd e mean são de 15 plantas.</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
Estes 3 pontos amostrais estão dentro de uma área de 1.42 ha com 1775 plantas.</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)"><br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
sd=0.2587516</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">mean= 3.013333</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">assim acho que o valor 10000 seria 1775, ok?</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Posso determinar o numero de amostras para para uma área de 20 ha com 37085 plantas?</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Grato</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<br></div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">Hélio</div><div style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:13px;background-color:rgb(255,255,255)">
<br></div></div><br><div class="gmail_quote">Em 27 de março de 2013 12:30, Hélio Gallo Rocha <span dir="ltr"><<a href="mailto:heliogallorocha@gmail.com" target="_blank">heliogallorocha@gmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Caros Alexandre e Walmes.<div><br></div><div>Estou usando este teste de decisão de número de amostras para ver se minha amostragem está ok.</div>
<div><br></div><div>Mas fiquei na dúvida sobre o valor 10000 da linha</div><div><div class="im">
<blockquote type="cite" style="border-left-width:2px;border-left-style:solid;border-left-color:rgb(204,204,204);padding:0px 1em"><font face="trebuchet ms,sans-serif"><span style="font-family:'courier new',monospace">pop.total <- rnorm(10000, mean=20, sd=15)<br>
<br></span></font></blockquote><div><br></div><div><br></div></div><div>Vou exemplificar com dados de previsão de safra em café, no caso, são relativos a altura da planta de café.</div><div>Estas medidas referem a coleta em 3 pontos da lavoura, cada ponto é composto por 5 plantas, então sd e mean são de 15 plantas.</div>
<div>Estes 3 pontos amostrais estão dentro de uma área de 1.42 ha com 1775 plantas.</div><div><br></div><div>sd=0.2587516</div><div>mean= 3.013333</div><div><br></div><div>assim acho que o valor 10000 seria 1775, ok?</div>
<div><br></div><div>Posso determinar o numero de amostras para para uma área de 20 ha com 37085 plantas?</div><div><br></div><div>Grato</div><div><br></div><div>Hélio</div><div><br></div><div><br></div><br><div class="gmail_quote">
Em 20 de março de 2013 18:49, ASANTOS [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4658689h27@n4.nabble.com" target="_blank">ml-node+s2285057n4658689h27@n4.nabble.com</a>></span> escreveu:<div><div class="h5">
<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Perfeito Walmes,<br>
<br>
Era isto mesmo que eu buscava,<br>
<br>
Abraço,<br>
<br>
Alexandre<br>
<br>
<br>
<div>Em 20/03/2013 17:31, Walmes Zeviani
escreveu:<br>
</div>
<blockquote style="border-left:2px solid #cccccc;padding:0 1em" 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 href="http://user/SendEmail.jtp?type=node&node=4658689&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</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" rel="nofollow" link="external" 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></fieldset>
<br>
<pre>_______________________________________________
R-br mailing list
<a href="http://user/SendEmail.jtp?type=node&node=4658689&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<pre 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: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)
<a href="http://user/SendEmail.jtp?type=node&node=4658689&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<a href="http://user/SendEmail.jtp?type=node&node=4658689&i=3" rel="nofollow" link="external" target="_blank">[hidden email]</a>
======================================================================
</pre>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4658689&i=4" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr noshade size="1" color="#cccccc">
<div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
<a href="http://r-br.2285057.n4.nabble.com/R-br-Teste-para-decisao-do-numero-otimo-de-amostra-tp4658684p4658689.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-Teste-para-decisao-do-numero-otimo-de-amostra-tp4658684p4658689.html</a>
</div>
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4" target="_blank">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></blockquote></div></div></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>