在数据集R可以找到x1和X2，n1,n2的相关系数

这位高人真是大好人，太感谢了！再次谢谢您！

data crackman;
input x1 x2@;
%num(&num.-1);
cards;
0 0
0 0
0 0.194524181
0 0.193628445
0 0.324881475
0 0.62184472
0 0.620721755
0 0.951667341
0 0.950444341
0 1.125868305
0 0.947999226
0.463629947 0.94677711
0.662146589 1.122072013
1.337229367 1.887472204
1.334961048 3.191768201
1.33269345 4.152501155
1.833266353 5.722478385
1.830842011 8.047885729
3.234993863 12.75743694
4.828810466 16.64131955
6.230338604 19.63580102
8.540719884 22.39520766
10.65538505 25.34358262
12.94672109 27.95276074
15.42105955 30.14475545
16.46129195 30.70560581
17.53262754 32.43165634
19.20164089 34.82412985
19.77129671 35.43639784
20.34921228 37.32067156
22.14289494 38.61075375
21.53021689 39.92961874
22.13351014 36.6747002
22.74547163 37.95048576
22.74073388 39.25470498
23.99583954 41.26761358
23.99100983 40.58462724
23.98618031 41.26070298
23.34707735 41.25724776
24.61984585 40.57433489
23.97169283 41.93735249
23.33272929 41.24688242
23.96203542 41.93039321
;
run;
proc sql;
select count(*) as num into:num from crackman;
quit;
proc corr data=crackman out=a(keep=x1 where=(x1<=1));
var _all_;
run;
proc sql ;
select x1 into:nu from a where _name_='x2';
quit;
data b;
set a;
n=&nu;
m=n-x1;
if m<0 then group+0 ;
else group+1;
obs=_n_-2;
where group=1;
run;
proc print data=b;
run;
%macro num(num);
%do i=1 %to &num.;
%let j=%eval(&i.+2);
x&j.=lag&i.(x2);
%end;
output;
%mend num;

结果里面的X1为R值
OBS为错位次数
是5

5# crackman

方法虽然不尽相同但总是大同小异.