Calcul avec la fonction optim sous R

Bonjour,

nous essayons minimiser la logvraisemblance d'une fonction gamma (translatée). Cependant nous n'arrivons pas a obtenir un résultat cohérent. Force est de constater que les résultats ne sont pas valides

(projet3 nous renvoie un vecteur de 22265 données toutes > a -7.92)

fc=function(par,x)
log(gamma(par[2]))-par[1]*log(par[2])-(par[1]-1)/length(x)*sum(log(x-7,92377))+par[2]/length(x)*sum(x-7,92377)

Estimmv=function(){
x=projet3()
emv=optim(c(100,100), fc, x=x, lower=c(0.001,0.001), upper=c(100000,100))
return(emv)
}

Merci d'avance d'essayer de nous aider