Matheus, você quer fazer especificamente obter uma amostra levando em conta a alocação proporcional não apenas do segundo grupo, mas em ralação aos 50 municípios, é isso velho? Se for, o exemplo que eu te passei não satisfaz, pois alí a alocação proporcional leva em conta a população apenas dos 40 municípios restantes.
Uma saída seria gerar o sample a partir dos 50 e então verificar se o município pertence ao primeiro ou ao segundo grupo. prop <- pop/sum(pop) sub1 <- pop[(1:10)] sub2 <- pop[(11:50)] set.seed(51) sub3 <- sample(pop, 10, replace=TRUE, prob=prop) # identifica match(sub1,sub3) match(sub2,sub3) 2013/4/29 Matheus Monteiro <[email protected]> > Alguém? > > > Em 29 de abril de 2013 01:26, Matheus Monteiro > <[email protected]>escreveu: > > Apenas uma dúvida, em sub3 nao seria: >> >> sub3 <- sample(sub1+sub2, 10, replace=TRUE, prob=prop) >> >> pois preciso trabalhar com os dados contidos nos dois estratos >> >> desculpe se estiver falando besteira >> >> >> Em 29 de abril de 2013 01:21, Matheus Monteiro >> <[email protected]>escreveu: >> >> Obrigado amigo. >>> >>> >>> Em 29 de abril de 2013 01:14, Daniel Marcelino >>> <[email protected]>escreveu: >>> >>> Isso tá com cara de dever de casa. >>>> >>>> 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) >>>> >>>> # order >>>> pop <- sort(pop, T) >>>> >>>> # Alocação >>>> sub1 <- pop[(1:10)] >>>> sub2 <- pop[(11:50)] >>>> >>>> # Alocação proporcional >>>> prop <- sub2/sum(sub2) >>>> >>>> set.seed(51) >>>> sub3 <- sample(sub2, 10, replace=TRUE, prob=prop) >>>> >>>> # Stats >>>> m <- c(mean(sub1),mean(sub3)) >>>> var <- c(var(sub1),var(sub3)) >>>> stddev <- sqrt(var) >>>> >>>> tab <- data.frame(m, var, stddev) >>>> rownames(tab) <- c("+Pop","-Pop") >>>> tab >>>> >>>> Daniel >>>> >>>> >>>> 2013/4/28 Matheus Monteiro <[email protected]> >>>> >>>>> Ainda trabalhando com os 50 municipios mais populosos do Brasil, >>>>> divido a população em dois estratos, onde no primeiro estrato esteja os 10 >>>>> mais populosos e no segundo os 40 restantes. >>>>> Preciso obter uma amostra com reposição de 10 municípios usando o >>>>> alocação proporcional e calcular a estimativa da média e variancia. >>>>> Alguém pode ajudar? >>>>> >>>>> 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) >>>>> >>>>> >>>>> >>>>> >>>>> _______________________________________________ >>>>> 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. >>>>> >>>> >>>> >>>> >>>> -- >>>> \begin{signature} >>>> Daniel Marcelino >>>> ☁ [email protected] >>>> ☎ (514) 343 6111 #3799 >>>> Skype: d.marcelino >>>> ✎ 3200 Jean Brillant, Office C5071 >>>> Montreal, QC; H3T 1N8 >>>> Canada >>>> \end{signature} >>>> >>>> "Small steps toward a much better world" >>>> >>>> _______________________________________________ >>>> 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. > -- \begin{signature} Daniel Marcelino ☁ [email protected] ☎ (514) 343 6111 #3799 Skype: d.marcelino ✎ 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada \end{signature} "Small steps toward a much better world"
_______________________________________________ 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.
