C orroborando o que o Benilton e o André colocaram, basta você utilizar o argumento replace = TRUE previsto na função sample() do R. Isso resolverá seu problema de implementação de uma amostragem aleatória simples com reposição.
Sugiro que dê uma estudada no help da função. ?sample Abs ___________________________________________________ *Rodrigo Sant'Ana* -- Mestre em Ciência e Tecnologia Ambiental - MCTA/UNIVALI Graduado em Oceanografia - CTTMar/UNIVALI Universidade do Vale do Itajaí - UNIVALI -- Em 26 de abril de 2013 11:39, Benilton Carvalho <[email protected]>escreveu: > 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. >
_______________________________________________ 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.
