ajustement à une loi gamma dans le cadre de GLM

Postez ici vos questions, réponses, commentaires ou suggestions - Les sujets seront ultérieurement répartis dans les archives par les modérateurs

Modérateur : Groupe des modérateurs

Marianne Caris
Messages : 6
Enregistré le : 13 Fév 2013, 13:43

ajustement à une loi gamma dans le cadre de GLM

Messagepar Marianne Caris » 16 Déc 2013, 21:05

Bonjour,
novice sous R je suis perplexe suite à ces messages d'erreur:
> fitdistr(coutMoyenRC,densfun="gamma")
Erreur dans optim(x = c(5.65, 5.9817, 0, 0, 16.66, 4.545, 11.2, 15.6075, : la valeur initiale dans 'vmin' n'est pas finie
> fitdistr(coutMoyenRC,densfun="gamma",lower=0.0001)
Erreur dans optim(x = c(5.65, 5.9817, 0, 0, 16.66, 4.545, 11.2, 15.6075, : L-BFGS-B nécessite des valeurs finies de 'fn'
> fitdistr(freq,densfun="gamma",lower=0.0001)
Erreur dans optim(x = c(7L, 7L, 0L, 0L, 20L, 5L, 13L, 18L, 13L, 0L, 0L, 11L, : L-BFGS-B nécessite des valeurs finies de 'fn'
> fitdistr(coutMoyenRC,densfun="gamma",list(shape = 1, rate = 0.1),lower=0.0001)
Erreur dans optim(x = c(5.65, 5.9817, 0, 0, 16.66, 4.545, 11.2, 15.6075, : non-finite finite-difference value [1]

quelle erreur commeté-je?

Merci pour votre aide

Marianne

dicko ahmadou
Messages : 444
Enregistré le : 21 Nov 2009, 20:15

Messagepar dicko ahmadou » 16 Déc 2013, 23:01

Salut

Sans tes données, on ne peut que rester perplexe.
Le code ci dessous qui vient de l'aide de fitdistrmarche bien.

Code : Tout sélectionner

require(MASS)
set.seed(123)
x <- rgamma(100, shape = 5, rate = 0.1)
fitdistr(x, "gamma")


Donc pour pouvoir en dire plus il faudrait que tu fournisses un example minimal que l'on peut reproduire (comme le code plus haut).
The best thing about being a statistician is that you get to play in everyone's backyard.
John Tukey

Myriam barry
Messages : 1
Enregistré le : 28 Juil 2017, 07:48

Re: ajustement à une loi gamma dans le cadre de GLM

Messagepar Myriam barry » 28 Juil 2017, 09:45

Bonjour,je suis confrontée au même problème et j'aimerai savoir comment l'avez résolu. Merci


Retourner vers « Questions en cours »

Qui est en ligne

Utilisateurs parcourant ce forum : Google [Bot] et 1 invité