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? Estou com dúvida aqui, pois preciso encontrar uma maneira de alocar proporcionalmente a população inteira e não apenas o estrato 2 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 DÚVIDA prop <- sub2/sum(sub2) set.seed(51) sub3 <- sample(sub2, 10, replace=TRUE, prob=prop)