logvraisemblance sous R

Bonjour
quelqu'un aurait l'amabilité de m'expliquer un peu le raisonnement effectué pour ce code svp,

je cherche à calculer la logvraisemblance de la fonction (échantillon distribué selon une loi normale) :

le code est :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22

LVS<‐function(x,n){ x<‐as.matrix(x,n,1) mu<‐seq(1,5,length=n) sigma<‐seq(1,5,length=n) T=nrow(x); lnL=matrix(0,n,n); for (i in 1:n){ for (j in 1:n){ mm2=sum((x‐mu[i])^2) lnL[i,j]=‐T*log(sigma[j])‐1/2*(mm2/(sigma[j]^2));}} persp(mu mu,sigma,lnL, theta = 70, phi = 30, expand = 0.5, col = "green"); max=min(lnL); for (i in 1:n){ for (j in 1:n){ if(max<lnL[i,j]){max=lnL[i,j]; indexi=i;indexj=j}}} mumax=mu[indexi];sigmamax=sigma[indexj]; return(list(mu=mu,sigma=sigma,lnL=lnL,sol=cbind(max,mumax,sigmamax)))}

---

j'ai mis en gras les lignes que je ne comprends vraiment pas, je suis prenante pour toute explication , merci beaucoup !!!

Bonjour,

Code:

---

1 2	mu<‐seq(1,5,length=n) sigma<‐seq(1,5,length=n)

---

C'est une suite entre 1 et 5 avec des sauts de n.

Code:

---

1 2 3 4

seq(1,5,1) 1 2 3 4 5 seq(1,5,2) 1 3 5

---

Code:

---

lnL=matrix(0,n,n);

---

Crée une matrice avec que des zéros.

Code:

---

persp(mu mu,sigma,lnL, theta = 70, phi = 30, expand = 0.5, col = "green");

---

Desine un graphique en 3 dimensions avec des paramètres pour l’orientation de la vues theta, phi et expand.

Cordialement.