童话牧场 发表于 2013-3-6 13:37
楼住,现在是否有蒙特卡罗模拟的SAS程序?
蒙特卡罗模拟法运用SAS软件的具体程序:
data eg1;
seed=12345;
array b(10000);
array x(16);
do j=1 to 10000;
do i=1 to 16;
y=normal(seed); /*产生服从标准正态分布的随机数*/
x(i)=y;
end;
array s(17);
s(1)=6.3482;
do i=1 to 16;
s(i+1)=s(i)+s(i)*(-0.000157/16+0.000881/4*x(i)); /*对s(i)进行迭代*/
end;
b(j)=s17; /*保留迭代结果*/
end;
do i=1 to 9999;
do j=i+1 to 10000;
if b(j)<b(i) then do;
temp=b(j);b(j)=b(i);b(i)=temp; /*升序排列*/
end;
end;
end;
proc print;
var b500; /*输出5%分位处的值*/
run;
以上是我很久以前写的,也不知道对错