Alocação proporcional no R
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)
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 <matheusmont@gmail.com>
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 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.
-- \begin{signature} Daniel Marcelino ☁ dm.silva@umontreal.ca ☎ (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"
Obrigado amigo. Em 29 de abril de 2013 01:14, Daniel Marcelino <dmarcelino@live.com>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 <matheusmont@gmail.com>
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 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.
-- \begin{signature} Daniel Marcelino ☁ dm.silva@umontreal.ca ☎ (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 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.
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 <matheusmont@gmail.com>escreveu:
Obrigado amigo.
Em 29 de abril de 2013 01:14, Daniel Marcelino <dmarcelino@live.com>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 <matheusmont@gmail.com>
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 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.
-- \begin{signature} Daniel Marcelino ☁ dm.silva@umontreal.ca ☎ (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 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.
Alguém? Em 29 de abril de 2013 01:26, Matheus Monteiro <matheusmont@gmail.com>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 <matheusmont@gmail.com>escreveu:
Obrigado amigo.
Em 29 de abril de 2013 01:14, Daniel Marcelino <dmarcelino@live.com>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 <matheusmont@gmail.com>
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 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.
-- \begin{signature} Daniel Marcelino ☁ dm.silva@umontreal.ca ☎ (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 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.
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 <matheusmont@gmail.com>
Alguém?
Em 29 de abril de 2013 01:26, Matheus Monteiro <matheusmont@gmail.com>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 <matheusmont@gmail.com>escreveu:
Obrigado amigo.
Em 29 de abril de 2013 01:14, Daniel Marcelino <dmarcelino@live.com>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 <matheusmont@gmail.com>
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 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.
-- \begin{signature} Daniel Marcelino ☁ dm.silva@umontreal.ca ☎ (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 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.
-- \begin{signature} Daniel Marcelino ☁ dm.silva@umontreal.ca ☎ (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"
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 <dmarcelino@live.com>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 <matheusmont@gmail.com>
Alguém?
Em 29 de abril de 2013 01:26, Matheus Monteiro <matheusmont@gmail.com>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 <matheusmont@gmail.com>escreveu:
Obrigado amigo.
Em 29 de abril de 2013 01:14, Daniel Marcelino <dmarcelino@live.com>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 <matheusmont@gmail.com>
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 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.
-- \begin{signature} Daniel Marcelino ☁ dm.silva@umontreal.ca ☎ (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 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.
-- \begin{signature} Daniel Marcelino ☁ dm.silva@umontreal.ca ☎ (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 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.
participantes (2)
-
Daniel Marcelino -
Matheus Monteiro