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...)
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.
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;
Merci de votre aide
Partager