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 <[email protected]> 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, <[email protected]> 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 < [email protected] > 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 >> [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. > > > > _______________________________________________ > 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. _______________________________________________ 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.
