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

[Matheus Monteiro]

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 em 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 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"
>
> _______________________________________________
> 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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130429/78cbbf41/attachment.html>