Formula em LaTeX: F_{T}(t;\alpha,\beta) = \Phi(\frac{1}{\alpha}(\frac{t}{\beta}^(1/2) - \frac{\beta}{t}^(1/2))) Pessoal com gerar números aleatórios com distribuição Birbaum-Saunder em C? Tenho a distribuição acumulada da Birbaum-Saunder que é escrita como acumulada de uma normal e tenho a densidade da Birbaum. A formula da acumulada está escrita acima para t e beta maiores que zero. Pensei primeiramente em pegar números aleatórios normais pela biblioteca GSL e fazer ser igual à t na formula acumulada acima. Daria certo? Se for pela transformada inversa eu teria que gerar números uniformes[0,1] e fazer F(t) = U. Depois eu tinha que ver qual o t em F me daria U, sendo U o número uniforme gerado. Outra complicação é que tenho uma acumulada da normal que é uma integral de -inf a t. Existe outra forma mais fácil? -- Saudações, Pedro Rafael Diniz Marinho.
O número de opções disponíveis para simular de uma distribuição depende do que você tem dela: 1) expressão da função de distribuição acumulada => transformada inversa 2) expressão da função densidade de probabilidade => amostragem por aceitação rejeição Veja um exemplo inocente de gerar v.a. com a função densidade sin(x)*I(0<x<pi) http://www.leg.ufpr.br/~walmes/ensino/ce083-2012-01/ce083-2012-01-aula12.R As opções vão além das que eu mencionei, depende do que você tem. No R felizmente já temos a distribuição implementada. Veja http://rss.acs.unt.edu/Rdoc/library/VGAM/html/bisaUC.html http://rss.acs.unt.edu/Rdoc/library/VGAM/html/bisa.html À 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 (2)
-
Pedro Rafael -
Walmes Zeviani