[R-br] Estratificação no R

[Benilton Carvalho]

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 em 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 em 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 em 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 em 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 em 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.