Se vc ler a documentacao do comando sample(), notara' que existe um argumento que vc pode configurar para que a amostragem seja feita com reposicao. b Em 26 de abril de 2013 11:11, Matheus Monteiro <matheusmont@gmail.com> escreveu:
Andre, obrigado mas preciso instalar o pacote gtools mesmo. Vou ter que encontrar uma maneira
Tem que implementar um AASc isto é com reposicao. A funcao sample faz por default um AASs. Revise isto
Abs
Matheus
Em 26 de abril de 2013 09:33, <andrebvs@bol.com.br> escreveu:
Veja se é isso:
# População:
pop <- c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761, 1555039,1416014,1410430,1333767,1244518,1098630,1039610,1016128,953393, 867067,830231,817590,805397,801746,774886,742428,680496,668877,654786, 643603,619746,619536,613815,600692,587701,568099,561329,526338,525225, 515707,491807,483821,474596,474596,472300,460062,446911,442701,433158, 425169,424948,422569)
Npop <- length(pop)
# Estratos:
EstI <- c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761, 1555039,1416014)
EstII <- c(1410430,1333767,1244518,1098630,1039610,1016128,953393,867067,
830231,817590,805397,801746,774886,742428,680496,668877,654786,643603,619746,
619536,613815,600692,587701,568099,561329,526338,525225,515707,491807,483821, Â 474596,474596,472300,460062,446911,442701,433158,425169,424948,422569)
nEstI<-length(EstI) nEstII<-length(EstII)
# Amostras:
set.seed(2013) A1 <- sample(EstI,5) A2 <- sample(EstII,5)
# Estimativas:
Médias <- c(mean(A1),mean(A2)) Variâncias <- c(var(A1),var(A2)) Desv_Pad <- sqrt(Variâncias)
Result <- data.frame(Médias,Variâncias,Desv_Pad) rownames(Result) <- c("EstI","EstII") Result
Att. André
________________________________ Em 26/04/2013 01:56, Matheus Monteiro < matheusmont@gmail.com > escreveu: Estou trabalhando com os 50 municipios mais populosos do Brasil, foi pedido para dividir a população em dois estratos, onde no primeiro estrato esteja os 10 mais populosos e no segundo o restante. Preciso obter uma amostra de 5 municípios de cada estrato e calcular a estimativa da media amostral e variancia. Sei que preciso instalar o pacote gtools mas não estou conseguindo. Não estou conseguindo, alguém pode ajudar? Por enquanto está assim:
pop = c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761,1555039,1416014,1410430,1333767,1244518,1098630,1039610,1016128,953393,867067,830231,817590,805397,801746,774886,742428,680496,668877,654786,643603,619746,619536,613815,600692,587701,568099,561329,526338,525225,515707,491807,483821,474596,474596,472300,460062,446911,442701,433158,425169,424948,422569)
Npop=length(pop)
Npop
EstI=c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761,1555039,1416014)
EstII=c(1410430,1333767,1244518,1098630,1039610,1016128,953393,867067,830231,817590,805397,801746,774886,742428,680496,668877,654786,643603,619746,619536,613815,600692,587701,568099,561329,526338,525225,515707,491807,483821,474596,474596,472300,460062,446911,442701,433158,425169,424948,422569)
nEstI=length(EstI)
nEstII=length(EstII)
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.