[R-br] Gerar multinomial

[Paulo Justiniano]

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
> 
> 
> 
>