Discussion :

[Geniuz77]

Bonjour, j'aimerais savoir comment effectuer l'opération suivante :

Je dois effectuer 1000 fois l'estimation de coefficients de régression (entre autres), et les placer dans une matrice vide, dans le but ultime de calculer la moyenne de ceux-ci. (Il fallait auparavant générer les variables en question...)

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
40
41
42
43
44
45
46
47
%Création des variables
N=50;  
 
x1=randn(N,1);
 
x2=randn(N,1);
 
Et=randn(N,1);
 
B0=0.5;
B1=1.2;
B2=0.8;
 
y=B0+(B1*x1)+(B2*x2)+Et;
 
T=N;
K=2;
 
%Création de la matrice X
X = [ones(size(x1)),x1,x2];
 
%Calcul des coefficient
b=inv(X'*X)*X'*y;
B=b';
 
%Calcul des résidus
e = y-X*b;
 
%Calcul de la variance des résidus
sigma2 = (e'*e)/(T-K);
 
%Calcul de la variance des coefficients
etype = sqrt(diag(var(b)));
 
 
 
%Créaion de la des matrices de résultats
 
Betas=NaN(size(b),m); 
 
m=1000
for i=(1:m);
    for j=(1:3)
 
 
end;
end;

Comme vous le constatez mon problème se situe au niveau de la boucle, j'aimerais avoir des informations à propos de la syntaxe a utiliser. Ma matrice vide sera de 3x1000 étant donné que ma matrice des coeficients et de 3x1.

Merci de votre aide

[Invité]

Bonsoir, et bienvenue sur le forum

Code :

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

1
2
3
Betas=NaN(length(b),m); 

m=1000

Attention de définir ta variable m avant de l'utiliser.

Je te conseille dans un premier temps la lecture de Qu'est-ce que la préallocation de mémoire ?
Préallocation qui est bien sûr extensible en 2D (NaN est aussi une solution, au même titre que zeros ou ones, à toi de voir selon le cas).

[tachmou]

Bonjour,

Où se situe ton problème exactement ? Tu as une erreur ?

Sinon, pour la programmation des boucles, je ne mettrais pas de point-virgule à la fin des lignes "for" et "end".

Cordialement,

[Cheetor]

Ayant vu ton programme, je me suis permise d'y effectuer quelques modifications :

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
40
41
42
43
 
% Listes des constantes
N=50;   % Taille des vecteurs d'entrées
B0=0.5; % Valeur du coeff 0
B1=1.2; % Valeur du coeff 1
B2=0.8; % Valeur du coeff 2
T=N;    % Affectation de la valeur N à la valeur T (utile ?)
K=2;    % utilisé dans une formule 
NbrCoeff = 3; % Correspond aux nombres de coeffs, à modifier s'il y en a plus ou moins
 
%Création de la matrice de résultats 
iter=1000;  % Nombre d'itérations souhaité
Betas=NaN(NbrCoeff,iter);   % Initialisation de la matrice Betas
 
for ind_iter=(1:iter)   
    % Vecteurs aléatoires (signaux d'entrées) modifié à chaque itération
    % Il est possible de les sortir de la boucle mais du coup, les colonnes
    % de Betas restent identiques
    x1=randn(N,1);
    x2=randn(N,1);
    Et=randn(N,1);
 
    % Formule liant les entrées à la sortie
    y=B0+(B1*x1)+(B2*x2)+Et;
 
    %Création de la matrice X
    X = [ones(size(x1)),x1,x2];
 
    %Calcul des coefficients
    b=inv(X'*X)*X'*y; % Obtention des coeff
    B=b'; % Passage en ligne (inutile si Betas est en 3x1000)
    Betas(:,ind_iter) = b; % Copie du résultat dans la colonne souhaitée
 
    %Calcul des résidus
    e = y-X*b;
 
    %Calcul de la variance des résidus
    sigma2 = (e'*e)/(T-K);
 
    %Calcul de la variance des coefficients
    etype = sqrt(diag(var(b)));
 
end

Tout ce qui correspond à des constantes ont été mises en début de programme.
Ensuite, si tu veux estimer 1000 fois des coefficients, il est préférable d'avoir des vecteurs qui se modifient à chaque itération. J'ai donc passé x1, x2 et Et dans la boucle afin qu'ils soient modifiés à chaque passage dans la boucle.

Comme tu as une modification à chaque itération, tu dois recalculer Y, ta matrice X,... donc ils sont placés dans la boucle.
Du coup, tu as tes 3 coefficients dans B et tu peux les copier dans la matrice Betas.

Après, j'ai peut-être mal compris le problème/demande.