[R-br] Gerar multinomial

Paulo Justiniano paulojus em leg.ufpr.br
Quarta Janeiro 11 10:30:21 BRST 2012


Certamente tTiago, faltou o argumento replace!!
obrigado por apontar a correção

On Wed, 11 Jan 2012, Tiago Vieira wrote:

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


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