A função pgamma faz justamente o que você mencionou, retorna o valor da probabilidade para o evento Pr(X<x) sendo X~Gamma(shape, scale). Nível de confiança não vem ao caso porque isso **não é uma estimativa** de probabilidade. Agora se você a partir de uma amostra aleatória da variável aleatória X que você considera ter distribuição Gamma você **estimou os parâmetros** shape e scale (por máxima verossimilhança, por exemplo), então probabilidades calculadas a partir deles podem ter incerteza associada, que pode ser representada por um intervalo/região de confiança por exemplo. Para isso você terá que obter os extremos dos valores de probabilidade para todos os valores de shape e scale internos à uma região de confiança. Procure um estatístico local para orientações mais aprofundadas. À 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================