Alguém?

Em 29 de abril de 2013 01:26, Matheus Monteiro <matheusmont@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@gmail.com> escreveu:

Obrigado amigo.

Em 29 de abril de 2013 01:14, Daniel Marcelino <dmarcelino@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@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@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@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@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@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.