Discussion :

[Invité]

Pouvez-vous en fait me corriger le code suivant :

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
16
17
18
19
20
21
22
23
24
25
26
27
x=c(1,2,3,4,5,6)
fx=c(1486,694,195,37,10,1)
moy <- sum(x*fx)/sum(fx)
 
 
theta=moy
i=0
repeat
 
{
n=sum(fx)
deriveesec=function(theta)
{
-sum(fx)/theta**2+(n*exp(-theta))/(1-exp(-theta))**2
}
i=i+1
thetaold=theta
derivee=-n+sum(fx)/theta-(n*exp(-theta))/(1-exp(-theta))
 
 
theta=theta-derivee/deriveesec
 
if (sqrt(sum((theta-thetaold)**2))<0.000001)
{break}
 
print(i)
}

l'erreur retournée est la suivante :

Error in derivee/deriveesec : non-numeric argument to binary operator

[mamounMob]

Code :

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

1
2
 
theta=theta-derivee/deriveesec

tu divises un numérique par une fonction !
pas étonnant qu'il t'affiche argument non numérique.

[Invité]

J'ai corrigé le code précédent par le suivant, les dérivées sont numériques :

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
16
17
18
19
20
21
22
23
theta=moy
i=0
repeat
 
{
n=sum(fx)
deriveesec={-sum(fx)/theta**2+(n*exp(-theta))/(1-exp(-theta))**2
}
i=i+1
thetaold=theta
derivee=
{
f=-n+sum(fx)/theta-(n*exp(-theta))/(1-exp(-theta))
}
 
 
theta=theta-derivee/deriveesec
 
if (sqrt(sum((theta-thetaold)**2))<0.001)
{break}
 
print(i)
}

Cependant, il est impossible d'afficher theta et le message d'erreur est le suivant :

Error in if (sqrt(sum((theta - thetaold)^2)) < 0.001) { :
missing value where TRUE/FALSE needed

[mamounMob]

Code :

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

1
2
3
4
5
6
 
# i = 8
> theta
[1] -892.6866
> exp(-theta)
[1] Inf

Exécute ton code à la main, itération par itération si il le faut

[Sengar]

Salut

Je dirais que tu t'es trompée dans le calcul de la dérivée.
A vue de nez j'aurais dis:

Code :

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

1
2
3
4
  deriveesec={-sum(fx)/theta**2 -
                (n*exp(-theta))/(1-exp(-theta))**2 +
                (n*exp(-theta))/(1-exp(-theta))
}

Là ça converge.

Mais c'est plus trop un problème de R

[Invité]

Merci beaucoup, le programme fonctionne, il restera maintenant à le comparer à la fonction implémentée sous R.
C'est curieux, pourtant je crois avoir recopié la formule que le prof nous a donnée. Je voulais donc m'assurer que tu as bien utilisé la formule de la dérivée première du dessous pour calculer la dérivée seconde à la main.

Sinon quelqu'un d'entre vous s'y connaît en SAS ? Car je n'ai pas de réponse sur forum pour le moment