[R-br] Estratificação no R

[Matheus Monteiro]

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.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130426/db0847cd/attachment.html>