Discussion :

[Gharbi1776]

Bonjour,
Pouvez-vous m'aider s'il vous plaît, ça fait des jours que j'essaie de faire marcher cet algorithme mais rien n'y fait. J'ai toujours des messages d'erreur, j'ai bien installé le package Nloptr , j'ai confirmé avec un exemple, il a marché mais j'ai toujours ces deux messages, je n'arrive pas à savoir d'où vient le problème.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 likelihood<-function (alphaopt,betaopt)
 { 
  return( - ( ( alphaopt ^ n ) * ( betaopt ^ n ) * prod(betaopt) * exp( ( - alphaopt ) * sum(betaopt) ) ) )
 }
 
  constraints<-function (alphaopt,betaopt) 
 {
  return( c( ( n / alphaopt ) - sum(betaopt), ( n / betaopt ) + sumlog(betaopt) - ( alphaopt * sumprod(betaopt) ) ) )
 }
 res<- c()
 x0<- c(0, 0)
 
 opts <- list("algorithm"="NLOPT_LD_LBFGS", "xtol_rel"=3)
 
 nloptr( x0=x0, eval_f=likelihood, eval_grad_f=constraints, opts=opts)

Error in .checkfunargs(eval_f, arglist, "eval_f") :
eval_f requires argument 'betaopt' but this has not been passed to the 'nloptr' function.

et les fonctions sum(betaopt), prod(betaopt) et sumlog(betaopt) ont bien été définies au début de l'algorithme.
Pouvez-vous m'aider svp ?

[dev_ggy]

Bonsoir,

Il est possible que les fonctions ne prennent pas plus d'un paramètre.

Dans ce cas au lieu de donner alphaopt,betaopt séparément tu prends un paramètre X = c(alphaopt,betaopt).

Je ne suis pas sûr que cela marche, faute de temps je ne peux essayer pour toi.

Bon courage.

[Gharbi1776]

Bonjour,
Tout d'abord, je vous remercie d'avoir répondu à ma demande, j'avais tenté cette possibilité mais ça n'avait pas marché mais suite à votre réponse j'ai refait une tentative avec quelques changements de variables, je suis parvenu à obtenir l'output ci-dessous, seulement, je ne comprends pas pourquoi le nombre d'itérations est de zéro.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
Call:
nloptr(x0 = opt0, eval_f = likelihood, eval_grad_f = constraints,     opts = opts)
 
 
Minimization using NLopt version 2.4.0 
 
NLopt solver status: -1 ( NLOPT_FAILURE: Generic failure code. )
 
Number of Iterations....: 1 
Termination conditions:  xtol_rel: 1e-04 
Number of inequality constraints:  0 
Number of equality constraints:    0 
Current value of objective function:  0 
Current value of controls: 0 0