Caros Alexandre e Walmes. Estou usando este teste de decisão de número de amostras para ver se minha amostragem está ok.
Mas fiquei na dúvida sobre o valor 10000 da linha pop.total <- rnorm(10000, mean=20, sd=15) Vou exemplificar com dados de previsão de safra em café, no caso, são relativos a altura da planta de café. 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. Estes 3 pontos amostrais estão dentro de uma área de 1.42 ha com 1775 plantas. sd=0.2587516 mean= 3.013333 assim acho que o valor 10000 seria 1775, ok? Posso determinar o numero de amostras para para uma área de 20 ha com 37085 plantas? Grato Hélio Em 27 de março de 2013 12:30, Hélio Gallo Rocha <[email protected]>escreveu: > Caros Alexandre e Walmes. > > Estou usando este teste de decisão de número de amostras para ver se minha > amostragem está ok. > > Mas fiquei na dúvida sobre o valor 10000 da linha > > pop.total <- rnorm(10000, mean=20, sd=15) > > > > Vou exemplificar com dados de previsão de safra em café, no caso, são > relativos a altura da planta de café. > 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. > Estes 3 pontos amostrais estão dentro de uma área de 1.42 ha com 1775 > plantas. > > sd=0.2587516 > mean= 3.013333 > > assim acho que o valor 10000 seria 1775, ok? > > Posso determinar o numero de amostras para para uma área de 20 ha > com 37085 plantas? > > Grato > > Hélio > > > > Em 20 de março de 2013 18:49, ASANTOS [via R-br] < > [email protected]> escreveu: > > Perfeito Walmes, >> >> Era isto mesmo que eu buscava, >> >> Abraço, >> >> Alexandre >> >> >> Em 20/03/2013 17:31, Walmes Zeviani escreveu: >> >> Alexandre, >> >> Segue um CMR simples: >> >> # tamanho de amostra para medida de interesse que >> # é o erro absoluto ao estimar a média da população >> # erro = abs(média_estimada-média_verdadeira) >> >> pop.total <- rnorm(10000, mean=20, sd=15) # população >> m <- mean(20) # média verdadeira >> i <- unique(floor(2^seq(1, 12, by=0.5))); i # tamanhos de amostra >> nr <- 100 # número de execuções de >> reamostragem >> fun <- function(size, pop, m){ >> # função que calcula o erro a partir de uma amostra >> x <- sample(pop, size, replace=TRUE) >> # retorna a diferença entre estimado e real >> abs(mean(x)-m) >> } >> >> results1 <- sapply(i, >> function(j){ >> replicate(nr, fun(size=j, pop=pop.total, m=m)) >> }) >> str(results1) # matriz: colunas correspondem à i >> results2 <- colMeans(results1) # tira a média para cada i >> results2 >> >> plot(results2~i) >> plot(results2~log2(i)) >> >> # qual o tamanho de amostra para ter um erro no ponto >> # médio entre o máximo (quando n=1) e o mínimo (n->infinito, erro=0) >> >> n0 <- nls(results2~A*B^(-log2(i)+1), start=list(A=1, B=1)) >> coef(n0) >> >> perc <- 0.3 # erro tolerado de 30% do máximo >> (usuário define) >> xi <- -log(perc)/log(coef(n0)["B"])+1 # tamanho de amostra correspondente >> >> plot(results2~log2(i), >> xlim=c(1, max(log2(i))), ylim=c(0,max(results2, coef(n0)["A"])), >> xlab="Tamanho da amostra (log2 unidades)", >> ylab="Erro absoluto médio ao estimar a média") >> abline(v=1, lty=2) >> with(as.list(coef(n0)), curve(A*B^(-x+1), add=TRUE)) >> with(as.list(coef(n0)), abline(h=A*c(1,perc), v=xi, lty=2, col=2)) >> >> ceiling(2^xi) # tamanho de amostra à ser considerado >> >> À disposição. >> Walmes. >> >> >> ========================================================================== >> Walmes Marques Zeviani >> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) >> Departamento de Estatística - Universidade Federal do Paraná >> fone: (+55) 41 3361 3573 >> VoIP: (3361 3600) 1053 1173 >> e-mail: [hidden email]<http://user/SendEmail.jtp?type=node&node=4658689&i=0> >> skype: walmeszeviani >> twitter: @walmeszeviani >> homepage: http://www.leg.ufpr.br/~walmes >> linux user number: 531218 >> ========================================================================== >> >> >> >> _______________________________________________ >> R-br mailing list[hidden email] >> <http://user/SendEmail.jtp?type=node&node=4658689&i=1>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código >> mínimo reproduzível. >> >> >> -- >> ====================================================================== >> 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)[hidden email] >> <http://user/SendEmail.jtp?type=node&node=4658689&i=2> >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658689&i=3> >> ====================================================================== >> >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4658689&i=4> >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >> código mínimo reproduzível. >> >> ------------------------------ >> If you reply to this email, your message will be added to the >> discussion below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Teste-para-decisao-do-numero-otimo-de-amostra-tp4658684p4658689.html >> To unsubscribe from R-br, click >> here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4> >> . >> NAML<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> >> > > > > -- > Hélio Gallo Rocha > IFSULDEMINAS - Câmpus Muzambinho > -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
