Discussion :

[ANOVA]

Bonjour,

je dois traduire un algorithme MATLAB en R. Seulement, je ne connais absolument pas Matlab et je bloque sur cette instruction :

V(:,iter+1)=max(aux,2)

Que fait elle ? Pouvez vous m'expliquer plus particulièrement max(aux,2), car j'ai trouvé ce que signifie V(:,iter+1) (on prend le iter+1 ième vecteur colonne de la matrice V).

V est une matrice
aux est une matrice
iter est une variable qui compte les itérations d'une boucle for

Le commentaire à côté de cette ligne de code est :

% take the max value over vij ignoring the missing values

Merci de votre aide.

PS : pour info, cet algorithme vient de la page 11 de ce document : matlab

[Caro-Line]

La version de MATLAB utilisée est-elle indiquée ?

Car c'est étrange dans ma version, max(aux,2) renverra une matrice de la taille de aux, avec les valeurs inférieures à 2 remplacées par 2 (en gros pour chque élément il prend le max entre aux et 2).
Donc là tu essayerais de faire rentrer une matrice dans un vecteur colonne (puisque ton interprétation est exacte).
Donc pour moi ce code ne devrait pas fonctionner...

Et ce ne serait pas plus simple d'avoir l'algo en algo plutôt que de passer par du MATLAB ?

[phryte]

Salut.

max(aux,2) renverra une matrice de la taille de aux, avec les valeurs inférieures à 2 remplacées par 2 (en gros pour chque élément il prend le max entre aux et 2).

Caro a raison. Exemple :

Code :

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

1
2
3
4
V=zeros(5,5);
iter=1;
aux=[0 1 2 3 4]';
V(:,iter+1)=max(aux,2)

Résultats :

Code :

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

1
2
3
4
5
6
7
V =
 
     0     2     0     0     0
     0     2     0     0     0
     0     2     0     0     0
     0     3     0     0     0
     0     4     0     0     0

[ANOVA]

Et ce ne serait pas plus simple d'avoir l'algo en algo plutôt que de passer par du MATLAB ?

Si je pouvais avoir l'algo plutôt que du code MATLAB, ça m'arrangerait, seulement, le seul document dont je dispose est celui ci :

matlab

et l'algo d'où vient la ligne de code qui pose problème est page 11.

La version de MATLAB utilisée est-elle indiquée ?

Non, la version de Matlab n'est pas indiquée, mais le document date de 2003...en plus, j'ai même pas matlab, donc je ne peux pas faire de tests moi même...

[ANOVA]

Je me permet de faire remonter ce sujet pour demander un petit service à quelqu'un qui aurait matlab. Je ne dispose pas de Matlab (malheureusement). Je ne peux donc pas faire les tests seul. Est ce que quelqu'un pourrait exécuter ce code et me dire s'il fonctionne et ce qu'il fait ?

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
28
29
30
31
32
33
34
35
36
37
38
39
clear
MaxIter=30;
beta=0.9;
K=0:0.01:1;
[rk,dimK]=size(K);
V=zero(dimK,MaxIter);
 
for iter=1:MaxIter
    iter
    aux=zeros(dimK,dimK)+NaN;
 
    for ik=1: dimK
        for ik2=1:(ik-1)
            aux(ik,ik2)=log(K(ik)-K(ik2))+beta*V(ik2,iter);
        end
    end
 
V(:,iter+1)=max(aux,2);
 
end
[Val,Ind]=max(aux,2);
 
optK=K(Ind);
optK=optK+Val*0;
optC=K'-optK';
 
figure(1)
plot(K,V);
xlabel('Size of Cake');
ylabel('Value Function');
 
 
% plot optimal consumption rule as a function of cake size
figure(2)
plot(K,[optC K'],'LineWidth',2)
xlabel('Size of Cake');
ylabel('Optimal Consumption');
text(0.4,0.65,'45 deg.line','FontSize',18)
text(0.4,0.13,'Optimal Consumption','FontSize',18)

Merci.

[Caro-Line]

Ben...il plante

Ligne 6 : il faut remplacer zero par zeros :

Code :

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

V=zeros(dimK,MaxIter);

Et (comme je l'avais prédit, je vais me reconvertir dans la voyance ), là :

Code :

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

V(:,iter+1)=max(aux,2);

donne :

??? Subscripted assignment dimension mismatch.

Pour les explications voir mon post #2.

Je suppose (avec un degré d'incertitude maximum cependant) que cela vient d'une mauvaise utilisation de la fonction MAX.
Il faudrait voir la documentation de cette fonction dans la version de MATLAB utilisée par l'auteur (2003 si on regarde ici :Quelles sont les différentes versions de MATLAB ?, on peut supposer que c'est la r13..., ayant seulement accès à la documentation de la R13SP2, ben l'aide de MAX est la même)

Ou alors comprendre ce qu'il a voulu faire avec son max.

[ANOVA]

Effectivement, cela vient peut être d'un problème de version. En allant sur ce site (http://www.mathworks.com/access/help...x+matlab&meta=), on trouve :

C = max(A,[],dim) returns the largest elements along the dimension of A specified by scalar dim. For example, max(A,[],1) produces the maximum values along the first dimension (the rows) of A.

Donc je pense que ce qu'il veut faire avec cette instruction, c'est récupérer la valeur maximale par rapport aux colonnes. Du coup, avec la dernière version de matlab, la syntaxe serait la suivante :

Code :

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

V(:,iter+1)=max(aux,[],2);

Est ce que vous pouvez tester le code avec cette modification et me dire si ça marche svp ?

Merci beaucoup à caro95470 pour sa rapidité.