[R-br] Gerar multinomial

Tiago Vieira tiagovieira7 em yahoo.com.br
Quarta Janeiro 11 11:30:00 BRST 2012


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
> 
> 
> 
>
_______________________________________________
R-br mailing list
R-br em listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120111/0a1e4b5d/attachment.html>


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