|
0
|
1 function x_dw=mz_6bis_durbin_watson(x_tr_mz)
|
|
|
2
|
|
|
3 // pas de calcul matriciel car pas assez de memoire
|
|
|
4 // x_dw doit varier entre 0 et 2
|
|
|
5 // 0=ce qu on recherche
|
|
|
6
|
|
|
7 [n,q]=size(x_tr_mz.d);
|
|
|
8
|
|
|
9 // calcul de la norme des differences et des vect-colonne de x_tr_mz
|
|
|
10 x_diff=x_tr_mz.d(2:n,:)-x_tr_mz.d(1:n-1,:);
|
|
|
11 rmsce=zeros(q,1);
|
|
|
12 rmsce2=zeros(q,1);
|
|
|
13 for i=1:q;
|
|
|
14 rmsce(i)=sqrt((x_diff(:,i)'*x_diff(:,i))/(n-1));
|
|
|
15 rmsce2(i)=sqrt((x_tr_mz.d(:,i)'*x_tr_mz.d(:,i))/(n));
|
|
|
16 end
|
|
|
17
|
|
|
18 x_dw=rmsce./rmsce2;
|
|
|
19
|
|
|
20
|
|
|
21 endfunction
|
