[R-br] Gerar multinomial

[Mauro Sznelwar]

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/20120110/cb51c24a/attachment.html>