Alguém?
Em 29 de abril de 2013 01:26, Matheus Monteiro <[email protected]>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 > <[email protected]>escreveu: > > Obrigado amigo. >> >> >> Em 29 de abril de 2013 01:14, Daniel Marcelino <[email protected]>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 <[email protected]> >>> >>>> 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 >>>> [email protected] >>>> 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 >>> ☁ [email protected] >>> ☎ (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 >>> [email protected] >>> 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 [email protected] 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.
