
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)

O pacote que não instala deve ser problema no R do meu Ubuntu, mas em relação ao que se pede como devo proceder? Em 26 de abril de 2013 01:56, Matheus Monteiro <matheusmont@gmail.com>escreveu:
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)

Isso parece mais um dever de casa do que dúvida. Pra fazer a amostragem por amostra estratificada tu deve separar em estratos (o que tu ja fez), calcular o tamanho amostral em cada estrato (que tu estipulou sendo 5) e fazer uma amostra aleatoria simples dentro de cada estrato (com o sample()). Para análise, recomendo o pacote survey 2013/4/26 Matheus Monteiro <matheusmont@gmail.com>
O pacote que não instala deve ser problema no R do meu Ubuntu, mas em relação ao que se pede como devo proceder?
Em 26 de abril de 2013 01:56, Matheus Monteiro <matheusmont@gmail.com>escreveu:
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)
_______________________________________________ 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.

Andre, obrigado mas preciso instalar o pacote gtools mesmo. Vou ter que encontrar uma maneira Tem que implementar um AASc isto é com reposicao. A funcao sample faz por default um AASs. Revise isto Abs Matheus Em 26 de abril de 2013 09:33, <andrebvs@bol.com.br> escreveu:
Veja se é isso:
*# População:*
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)
*# Estratos:*
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)
*# Amostras:*
set.seed(2013) A1 <- sample(EstI,5) A2 <- sample(EstII,5)
*# Estimativas:*
Médias <- c(mean(A1),mean(A2)) Variâncias <- c(var(A1),var(A2)) Desv_Pad <- sqrt(Variâncias)
Result <- data.frame(Médias,Variâncias,Desv_Pad) rownames(Result) <- c("EstI","EstII") Result
*Att.* *André*
------------------------------ Em 26/04/2013 01:56, *Matheus Monteiro < matheusmont@gmail.com >*escreveu: 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)
_______________________________________________ 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.

Se vc ler a documentacao do comando sample(), notara' que existe um argumento que vc pode configurar para que a amostragem seja feita com reposicao. b Em 26 de abril de 2013 11:11, Matheus Monteiro <matheusmont@gmail.com> escreveu:
Andre, obrigado mas preciso instalar o pacote gtools mesmo. Vou ter que encontrar uma maneira
Tem que implementar um AASc isto é com reposicao. A funcao sample faz por default um AASs. Revise isto
Abs
Matheus
Em 26 de abril de 2013 09:33, <andrebvs@bol.com.br> escreveu:
Veja se é isso:
# População:
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)
# Estratos:
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)
# Amostras:
set.seed(2013) A1 <- sample(EstI,5) A2 <- sample(EstII,5)
# Estimativas:
Médias <- c(mean(A1),mean(A2)) Variâncias <- c(var(A1),var(A2)) Desv_Pad <- sqrt(Variâncias)
Result <- data.frame(Médias,Variâncias,Desv_Pad) rownames(Result) <- c("EstI","EstII") Result
Att. André
________________________________ Em 26/04/2013 01:56, Matheus Monteiro < matheusmont@gmail.com > escreveu: 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)
_______________________________________________ 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.

C orroborando o que o Benilton e o André colocaram, basta você utilizar o argumento replace = TRUE previsto na função sample() do R. Isso resolverá seu problema de implementação de uma amostragem aleatória simples com reposição. Sugiro que dê uma estudada no help da função. ?sample Abs ___________________________________________________ *Rodrigo Sant'Ana* -- Mestre em Ciência e Tecnologia Ambiental - MCTA/UNIVALI Graduado em Oceanografia - CTTMar/UNIVALI Universidade do Vale do Itajaí - UNIVALI -- Em 26 de abril de 2013 11:39, Benilton Carvalho <beniltoncarvalho@gmail.com>escreveu:
Se vc ler a documentacao do comando sample(), notara' que existe um argumento que vc pode configurar para que a amostragem seja feita com reposicao. b
Em 26 de abril de 2013 11:11, Matheus Monteiro <matheusmont@gmail.com> escreveu:
Andre, obrigado mas preciso instalar o pacote gtools mesmo. Vou ter que encontrar uma maneira
Tem que implementar um AASc isto é com reposicao. A funcao sample faz por default um AASs. Revise isto
Abs
Matheus
Em 26 de abril de 2013 09:33, <andrebvs@bol.com.br> escreveu:
Veja se é isso:
# População:
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)
# Estratos:
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)
# Amostras:
set.seed(2013) A1 <- sample(EstI,5) A2 <- sample(EstII,5)
# Estimativas:
Médias <- c(mean(A1),mean(A2)) Variâncias <- c(var(A1),var(A2)) Desv_Pad <- sqrt(Variâncias)
Result <- data.frame(Médias,Variâncias,Desv_Pad) rownames(Result) <- c("EstI","EstII") Result
Att. André
________________________________ Em 26/04/2013 01:56, Matheus Monteiro < matheusmont@gmail.com > escreveu: 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)
_______________________________________________ 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.
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.

Muito obrigado pessoal, agradeço. Em 26 de abril de 2013 14:48, <andrebvs@bol.com.br> escreveu:
Utilize o seguinte:
A1 <- sample(EstI,5,replace=T) # com reposição A2 <- sample(EstII,5,replace=T) # com reposição
*Att.* *André*
------------------------------ Em 26/04/2013 13:04, *Rodrigo Sant'Ana < rodrigo.gringo@gmail.com >*escreveu: C orroborando o que o Benilton e o André colocaram, basta você utilizar o argumento replace = TRUE previsto na função sample() do R. Isso resolverá seu problema de implementação de uma amostragem aleatória simples com reposição.
Sugiro que dê uma estudada no help da função. ?sample
Abs
___________________________________________________ *Rodrigo Sant'Ana* -- Mestre em Ciência e Tecnologia Ambiental - MCTA/UNIVALI Graduado em Oceanografia - CTTMar/UNIVALI Universidade do Vale do Itajaí - UNIVALI --
Em 26 de abril de 2013 11:39, Benilton Carvalho < beniltoncarvalho@gmail.com<http://../../../undefined/compose?to=beniltoncarvalho@gmail.com>
escreveu:
Se vc ler a documentacao do comando sample(), notara' que existe um argumento que vc pode configurar para que a amostragem seja feita com reposicao. b
Em 26 de abril de 2013 11:11, Matheus Monteiro <matheusmont@gmail.com<http://../../../undefined/compose?to=matheusmont@gmail.com>> escreveu:
Andre, obrigado mas preciso instalar o pacote gtools mesmo. Vou ter que encontrar uma maneira
Tem que implementar um AASc isto é com reposicao. A funcao sample faz por default um AASs. Revise isto
Abs
Matheus
Em 26 de abril de 2013 09:33, <andrebvs@bol.com.br<http://../../../undefined/compose?to=andrebvs@bol.com.br>> escreveu:
Veja se é isso:
# População:
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,52633 8,525225,
515707,491807,483821,474596,474596,472300,460062,446911,442701,433158,
425169,424948,422569)
Npop >> # Estratos:
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)
# Amostras:
&g t;> set.seed(2013)
A1 >> A2 >> # Estimativas:
Médias >> Variâncias >> Desv_Pad >> Result >> rownames(Result) >> Result
Att. André
________________________________ Em 26/04/2013 01:56, Matheus Monteiro < matheusmont@gmail.com<http://../../../undefined/compose?to=matheusmont@gmail.com>> escreveu: 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 inst alar 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,613 815,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)
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br<http://../../../undefined/compose?to=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<http://../undefined/compose?to=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<http://../../../undefined/compose?to=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.
participantes (5)
-
andrebvs@bol.com.br
-
Benilton Carvalho
-
Matheus Monteiro
-
Rodrigo Coster
-
Rodrigo Sant'Ana