[R-br] simulacao gama

[Eliardo Costa]

Natalia,

Tente também o pacote "Runuran".

Att, Eliardo.

Em 25 de outubro de 2012 14:29, Walmes Zeviani
<walmeszeviani em gmail.com>escreveu:

> Inicialmente posso sugerir 2 abordagens
>
> # simulando da gamma e truncando entre 3 e 10
> curve(dgamma(x, shape=3, scale=2), 0, 20)
> x <- rgamma(100, shape=3, scale=2)
>
> # opção 1, truncar após simular
> x <- x[x>=3 & x<=10]
> length(x) # terá que simular mais que 100 para ter n útil > 100
>
> # opção 2, simular após truncar (transformação integral da probabilidade)
> curve(pgamma(x, shape=3, scale=2), 0, 20)
> xtrun <- c(3,10)
> ytrun <- pgamma(xtrun, shape=3, scale=2)
> abline(v=xtrun, h=ytrun)
>
> u <- runif(100, ytrun[1], ytrun[2])
> x <- qgamma(u, shape=3, scale=2)
>
> curve(pgamma(x, shape=3, scale=2), 0, 20)
> rug(x)
> rug(u, side=2)
>
> curve(dgamma(x, shape=3, scale=2), 0, 20)
> rug(x)
> abline(v=xtrun, col=2)
>
> À 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 em ufpr.br
> skype: walmeszeviani
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
> _______________________________________________
> 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/20121025/a57200d5/attachment-0001.html>