Obrigado Daniel, com o auxilio do seu exemplo anterior, eu usei os seguintes comandos # Alocação sub1 <- pop[(1:10)] sub2 <- pop[(11:50)]
# Alocação proporcional aux1 <- sample(sub2, 8 ) aux2 <- sample(sub1, 2) set.seed(51) sub3 <- c(aux1,aux2) vou compara-los. abraço Em 29 de abril de 2013 17:19, Daniel Marcelino <[email protected]>escreveu: > 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. >
_______________________________________________ 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.
