有懂iml的么 论文急用啊 求各位sas高高手看看这个程序为什么有错？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

根据《金融计算与建模:理论、算法与SAS程序》数里面对于First to Default的一篮子信用违约互换定价的例子
例子里面是十个债券
我给缩减到了3个  同时变换了些数字  整体程序应该没问题
但是始终跑不出来结果
求帮忙啊  
毕业论文急用。。

程序如下：

proc iml;
cov={1.0000    0.4015    0.4193
      0.4015    1.0000    0.3639
      0.4193    0.3639    1.0000};
h={0.07,0.25,0.12};
T={[10000]0};
V={[10000]0};
r=0.03;
value=0;
rv=rannor(repeat(1000,10000,3));
Y=rv*(root(cov));
i=1;
j=1;
do while(i<=10000);
    do while(j<=3);
      x=cdf('Normal',Y[i,j]);
      Y[i,j])=-(log(1-x))/h[j];
      j=j+1;
      end;
T[i]=min(Y[i,]);
if T[i]<6 then v[i]=50000*exp(-r*T[i]);
    value=value+v[i];
    j=1;
    i=i+1;
end;
value=value/10000;
print value;
run;
quit;

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：论文急用 高高手 金融计算与建模 Default repeat 程序 论文

呜呜~~

proc iml;
cov={1.0000    0.4015    0.4193,   /*协方差矩阵，用逗号换行*/
      0.4015    1.0000    0.3639,
      0.4193    0.3639    1.0000};
h={0.07,0.25,0.12};
T={[10000]0};
V={[10000]0};
r=0.03;
value=0;
rv=rannor(repeat(1000,10000,3));
Y=rv*(root(cov));
i=1;
j=1;
do while(i<=10000);
    do while(j<=3);
      x=cdf('Normal',Y[i,j]);
      Y[i,j]=-(log(1-x))/h[j];  /*多了一个括号*/
      j=j+1;
      end;
T=min(Y[i,]);
if T<6 then v=50000*exp(-r*T);
    value=value+v;
    j=1;
    i=i+1;
end;
value=value/10000;
print value;
quit;

ziyenano 发表于 2013-1-8 21:11
proc iml;
cov={1.0000    0.4015    0.4193,   /*协方差矩阵，用逗号换行*/
      0.4015    1.0000     ...

天呐！！！
太感谢了！！！
您吧地址给我  我给您寄两斤铁棍山药表示感谢！！！
真心的啊！！

urmine5683 发表于 2013-1-8 21:19
天呐！！！
太感谢了！！！
您吧地址给我  我给您寄两斤铁棍山药表示感谢！！！

呵呵，不用了哇~心领了！

ziyenano 发表于 2013-1-8 21:21
呵呵，不用了哇~心领了！

这怎么行！知恩图报！！

私信我吧  或者旺旺找我：铁棍儿山药

我家人自己种的 保证正宗的 哈哈哈

ziyenano 发表于 2013-1-8 21:21
呵呵，不用了哇~心领了！

ziyenano你不要客气了，给她地址吧，她已经给我寄了两斤了。

Imasasor 发表于 2013-1-8 21:29
ziyenano你不要客气了，给她地址吧，她已经给我寄了两斤了。

版主是来捧场的么~包邮不

ziyenano 发表于 2013-1-8 21:33
版主是来捧场的么~包邮不

0元包邮，她家种了几十亩山药

Imasasor 发表于 2013-1-8 21:35
0元包邮，她家种了几十亩山药

这样啊~那我就不客气了，比论坛币实在，谢谢楼主了