proc iml;
mu={0,0};
sig={1 0.9, 0.9 1};
a=half(sig)`;
seed=2343543;
p=2;
n=30;
x=rannor(j(p,n,seed));
y=a*x+mu*j(1,n,1);
ybar=y[,:];
yc=y-ybar*j(1,n,1);
s=yc*yc`/(n-1);
val=eigval(s);
minval=min(val);
print minval;
run;

|
楼主: guzhou88
|
1658
0
[原创博文] SAS中如何实现循环语句? |
加好友,备注cda京ICP备16021002号-2 京B2-20170662号
京公网安备 11010802022788号
论坛法律顾问:王进律师
知识产权保护声明
免责及隐私声明


