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