[R-br] Estratificação no R

[Matheus Monteiro]

Estou trabalhando com os 50 municipios mais populosos do Brasil, foi pedido
para dividir a população em dois estratos, onde no primeiro estrato esteja
os 10 mais populosos e no segundo o restante.
Preciso obter uma amostra de 5 municípios de cada estrato e calcular a
estimativa da media amostral e variancia.
Sei que preciso instalar o pacote gtools mas não estou conseguindo.
Não estou conseguindo, alguém pode ajudar?
Por enquanto está assim:

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)

Npop=length(pop)

Npop

EstI=c(11376685,6390290,2710968,2648532,2500532,2395785,1861838,1776761,1555039,1416014)

EstII=c(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)

nEstI=length(EstI)

nEstII=length(EstII)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130426/348c729f/attachment.html>