Gostaria de saber onde posso achar comandos do R para geração de números aleatórios da distribuição multinomial. Não posso usar o comando rmultinom que existe, preciso implementar algo para três variáveis, teta1, teta2 e teta3. Desde já agradeço.
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 On Mon, 9 Jan 2012, Mauro Sznelwar wrote:
Gostaria de saber onde posso achar comandos do R para geração de números aleatórios da distribuição multinomial. Não posso usar o comando rmultinom que existe, preciso implementar algo para três variáveis, teta1, teta2 e teta3. Desde já agradeço.
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
Mauro, É boa prática, faz parte do aprendizado e quase obrigação de todo usuário consultar sempre a documentação das funções. Ela existe para isso. Às vezes na digitação de códigos direto na mensagem (código não avaliado no R, enviado de supetão) acontece de errarmos a digitação ou coisas semelhantes. Pense por exemplo que pelo horário o Professor Paulo só vai te responder amanhã. E aí, não é mais rápido olhar a documentação? À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
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
olá R-br, Gostaria de saber se algum de vocês possui o artigo intitulado "A Homogeneity Test Applied to Precipitation Data" de H. Alexandersson (1986) para me fornecer. Obrigado pela ajuda. Até mais.
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
_______________________________________________ R-br mailing list R-br@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.
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@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
_______________________________________________ R-br mailing list R-br@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.
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
As probabilidades dentro prob=c() tem que somar 1? Complementando, assim dará certo: table(sample(1:3, size=100, replace = T, prob=c(0.1, 0.6, 0.3)))
nao, sao pesos a documentacao da funcao explica isto ?sample vai te informar com mais detalhes On Thu, 12 Jan 2012, Mauro Sznelwar wrote:
As probabilidades dentro prob=c() tem que somar 1? Complementando, assim dará certo: table(sample(1:3, size=100, replace = T, prob=c(0.1, 0.6, 0.3)))
Mauro, Isso está descrito na sessão *Details* da documentação da função sample(). Pro favor, consulte-a. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (5)
-
Lucas Venezian Povoa -
Mauro Sznelwar -
Paulo Justiniano -
Tiago Vieira -
Walmes Zeviani