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@leg.ufpr.br> Para: r-br@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