[R-br] Simulação Gama e weibull

Bernardo Rangel Tura tura em centroin.com.br
Quinta Agosto 2 06:04:14 BRT 2012


On 08/01/2012 07:03 PM, Simone toni ruiz corrêa wrote:
> Boa noite,
> Gostaria de, primeiramente, parabenizar pela iniciativa desta lista!!!
> Sou "novata" no R e estou trabalhando com simulação de temperatura média
> do ar e radiação solar usando as *distribuições gama e weibull*.
> As simulações estão ok, mas para descrição em Material e métodos de
> minha Tese, necessito saber *qual a metodologia usada pelo R nestas
> simulações*.
> Por exemplo, se o R utiliza a simulação inversa ou direta e se há artigo
> referencia neste assunto.
> Neste mesmo assunto, quando verifiquei a *aderencia* (teste de
> kolmogorov-smirnov)**da dist. weibull aos dados de radiação solar (série
> diária), muitos valores apresentaram a mensagem *optimization failled*.
> Alguém saberia me explicar o que acontece?
> Agradeço antecipadamente,
> Simone
>
Bom dia Simone

Primeiro seja bem vinda!

Segundo um grade macete é usar o help do R, por exemplo se digitar o 
comando ?rgamma você irá encontrar

     ‘rgamma’ for ‘shape >= 1’ uses

      Ahrens, J. H. and Dieter, U. (1982).  Generating gamma variates by
      a modified rejection technique.  _Communications of the ACM_,
      *25*, 47-54,

      and for ‘0 < shape < 1’ uses

      Ahrens, J. H. and Dieter, U. (1974).  Computer methods for
      sampling from gamma, beta, Poisson and binomial distributions.
      _Computing_, *12*, 223-246.

do mesmo jeito ?rweibull

Source:

      ‘[dpq]weibull’ are calculated directly from the definitions.
      ‘rweibull’ uses inversion.

References:

      Johnson, N. L., Kotz, S. and Balakrishnan, N. (1995) _Continuous
      Univariate Distributions_, volume 1, chapter 21.  Wiley, New York.

Terceiro:

Permita um pequeno comentário, na língua portuguesa o termo cero é 
observância e não aderência, porém no caso em questão provavelmente 
estamos discutindo ajuste.

Quanto ao erro fica difícil ajudá-la sem ver o comando e sem uma pequena 
amostra de seus dados. É o chamando exemplo mínimo ...

[]s
Tura




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