Discussion :

[regall]

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 :

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
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 !!!

[dev_ggy]

Bonjour,

Code :

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

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 :

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

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

Code :

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

lnL=matrix(0,n,n);

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

Code :

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

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.