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

[Matheus Monteiro]

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