Em 27 de março de 2013 12:30, Hélio Gallo Rocha <heliogallorocha@gmail.com> 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] <ml-node+s2285057n4658689h27@n4.nabble.com> 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]
skype: walmeszeviani
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
_______________________________________________
R-br mailing list
[hidden email]
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] 
        [hidden email]                    
======================================================================
_______________________________________________
R-br mailing list
[hidden email]
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.
NAML
--
Hélio Gallo Rocha
IFSULDEMINAS - Câmpus Muzambinho