[R-br] Alocação proporcional no R

[Daniel Marcelino]

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 em gmail.com>

> Alguém?
>
>
> Em 29 de abril de 2013 01:26, Matheus Monteiro <matheusmont em 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 em gmail.com>escreveu:
>>
>> Obrigado amigo.
>>>
>>>
>>> Em 29 de abril de 2013 01:14, Daniel Marcelino <dmarcelino em 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 em 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 em 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 em 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 em 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 em 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 em 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"
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130429/22738309/attachment.html>