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
>