[R-br] Gerar multinomial

Mauro Sznelwar sznelwar em uol.com.br
Quarta Janeiro 11 22:39:25 BRST 2012


Muito obrigado! Agora deu certo!

  Complementando, assim dará certo:


  table(sample(1:3, size=100, replace = T, prob=c(0.1, 0.6, 0.3)))




------------------------------------------------------------------------------
  De: Paulo Justiniano <paulojus em leg.ufpr.br>
  Para: r-br em listas.c3sl.ufpr.br 
  Enviadas: Quarta-feira, 11 de Janeiro de 2012 8:32
  Assunto: Re: [R-br] Gerar multinomial


  troque

  > table(sample(1:3, n=100, prob=c(0.1, 0.6, 0.3)))

  por

  table(sample(1:3, size=100, prob=c(0.1, 0.6, 0.3)))


  On Tue, 10 Jan 2012, Mauro Sznelwar wrote:

  > Muito obrigado pelo retorno professor! Mas não consegui rodar o que em enviou, veja o que aconteceu no meu exemplo:
  >  
  > table(sample(1:3, n=100, prob=c(0.1, 0.6, 0.3)))
  > Erro em sample(1:3, n = 100, prob = c(0.1, 0.6, 0.3)) :
  >   argumento(s) não utilizado(s) (n = 100)
  > 
  > Suponha que o vetor (X1, X2, X3) sao contagens de uma distr. multinomial
  > com probabilidades (p1, p2, p3) e que voce queira fixar o total
  > N = X1+X2+X3
  > 
  > neste caso pode usar
  > 
  > table(sample(1:3, n=N, prob=c(p1, p2, p3)))
  > 
  > substituindo N, p1, p2, p3 pelo valores adequados
  > 
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120111/ed3013ea/attachment.html>


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