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: