[R-br] Estratificação no R

Rodrigo Sant'Ana rodrigo.gringo em gmail.com
Sexta Abril 26 13:04:06 BRT 2013


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 em 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 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.
> _______________________________________________
> 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/a2a583cf/attachment.html>


Mais detalhes sobre a lista de discussão R-br