Algorithme de Metropolis Hastings

NaceraTizi

Bonjour tous le monde
Pour estimer le paramètre beta dans un modèle statistique (Delt, alp, delt et sig sont les autres paramètres de notre modèle et fvraisRCA est la fonction de vraisemblance), j'ai utilisé l'algorithme de Metropolis Hastings. Je vous donne la partie de mon programme Matlab ou il affiche l'erreur ainsi  que le message d'erreur qu'il affiche. Si vous pouvez m'aider. 

for i=2:N    

betastar=beta(i-1)+randn(1)*sig1;

alpha1=min(1,fvraisRCA(betastar,alp(i-1),Delt(i-1),delt(i-1),sig(i-1))/fvraisRCA(beta(i-1),alp(i-1),Delt(i-1),delt(i-1),sig(i-1)));

u1=rand(1);

 if u1<alpha1

     beta(i)=betastar;

 else

     beta(i)=beta(i-1);

 end

Out of memory. The likely cause is an infinite

recursion within the program.

Error in fvraisRCA (line 43)

alpha1=min(1,fvraisRCA(betastar,alp(i-1),Delt(i-1),delt(i-1),sig(i-1))/fvraisRCA(beta(i-1),alp(i-1),Delt(i-1),delt(i-1),sig(i-1)));

 Merci

gerard0

Bonjour.

La deuxième erreur te dit qu'il se passe quelque chose dans fvraisRCA, donc c'est dans cette routine qu'il faut aller voir.

Cordialement.