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
|
%macro ... ;
.
.
.
%let ind=%eval(&nbr_suj-1);
%do %until (&ind = 1) ;
/*initialisation*/
proc iml;
rvcc=shape(0,%eval(&nbr_suj-2),&nbr_sim);print rvcc ;quit;
%do sim=1 %to &nbr_sim;
.
.
obtention des tables config et config_k
.
.
.
.
proc iml;
use config;
read all var _num_ into mat;print mat ;
use config_k;
read all var _num_ into mat2; print mat2;
w=mat*T(mat);
w2=mat2*T(mat2);s=w*w2;
m=trace(s);
m2=sqrt(trace(w**2)*trace(w2**2));
k=%eval(%eval(&nbr_suj)-%eval(&ind));
rvCC[k,&sim]=m/m2;
quit;
%end;
%let ind=%eval(&ind-1);
%end;
.
.
.
.
.
%mend; |
Partager