[R-br] Teste para decisão do número ótimo de amostra
ASANTOS
alexandresantosbr em yahoo.com.br
Quarta Março 20 14:01:56 BRT 2013
Boa tarde pessoal,
Gostaria de utilizar um teste para decidir o número ótimo de
amostral, baseado na abordagem que estou empregando, se faço:
pop.total<-rnorm(10000,mean=20,sd=15) ## População total
i<-c(5000,50)## Quantidade de amostras
n.sim<-999 ## Número de simulações
RES=NULL
for(Nsim in 1:n.sim){
for(Nam in 1:length(i)){
amostras<-sample(pop.total,i[Nam]) ## Amostragem com reposição
m.amostras<-mean(amostras) ## Media encontrada nas amostras
err.amost<-100-((mean(amostras)/mean(pop.total))*100) ### Erro amostral
RES=rbind(RES,c(Nsim,i[Nam],err.amost))
}}
####### 50 amostras
am.50<-RES[RES[,2]==50,]
shapiro.test(am.50[,3]-mean(am.50[,3])) ## Distribuição normal,então
inferências sobre os parâmetros média e desvio estão corretos
mean(am.50[,3])## Média do erro amostral
sd(am.50[,3])/sqrt(length(am.50[,3]))## Erro padrão do erro amostral
sd(am.50[,3])## Desvio padrão
#
####### 5000 amostras
am.5000<-RES[RES[,2]==5000,]
shapiro.test(am.5000[,3]-mean(am.5000[,3])) ## Distribuição normal,então
inferências sobre os parâmetros média e desvio estão corretos
mean(am.5000[,3])## Média do erro amostral
sd(am.5000[,3])/sqrt(length(am.50[,3]))## Erro padrão do erro amostral
sd(am.5000[,3])## Desvio padrão
#
Percebo que a média do erro amostral e o erro padrão desta
média são baixos (próximos de zero) tanto quando minha amostra tem
tamanho 50 quanto 5000, porém o desvio padrão, obviamente, aumenta
quando reduzo o número de amostras. Gostaria de saber que teste poderia
empregar para escolha de um número ótimo de amostras, uma vez que eu não
posso usar como critério o erro amostral, quando simulo este cenário 999
vezes e emprego inúmeros tamanhos amostrais?
Obrigado,
--
======================================================================
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)
e-mails:alexandresantosbr em yahoo.com.br
alexandre.santos em cas.ifmt.edu.br
======================================================================
Mais detalhes sobre a lista de discussão R-br