<p>hull and white（1987）提出的期权定价模型如下：</p><p>dlogs=(mu-0.5*sigma^2)*dt+sigma*dw1&nbsp; --------(1)</p><p>d(sigma^2)=(a+b*sigma^2)dt+k*sigma*dw2------(2)</p><p>cov(dw1,dw2)=p---------------------------------------(3)</p><p>其中dt=1/52(数据为周数据)，换成离散形式后：</p><p>dlogs=(mu-0.5*sigma^2)*dt+sigma*sqr(dt)*e1&nbsp; --------(1)'</p><p>d(sigma^2)=(a+b*sigma^2)dt+k*sigma*sqr(dt)*e2------(2)'</p><p>cov(dw1,dw2)=p---------------------------------------(3)'</p><p>其中(1)'为量测方程，（2）'为状态方程，e1,e2服从N(0,1)的扰动项</p><p>数据dlogs在附件中，或任取一股票的价格，进行下dlog运算即可。</p><p>&nbsp;下面是本人的程序，能运行但肯定有错误，敬请高手指点！将附上必要的酬劳</p><p>&nbsp;!tt= 1/52</p><p>!sqrtt=sqr(!tt)</p><p><br/>'c(1)=mu c(2)=a c(3)=b c(4)=cox st=sigma^2<br/>'specify quasi-likelihood state-space<br/>sspace sv<br/>sv.append @signal dlogs =(c(1)-0.5*st) *!tt+&nbsp; [ename=e1]<br/>sv.append @state st =st(-1)+(c(2)+c(3)*st(-1))*!tt+ [ename=e2,var=1]<br/>@evar var(e1)=sqr(st)*!sqrtt*[var = 1]<br/>'<a href="mailto:'@evar">@evar</a> var(e2)=c(4)*sqr(st)*!sqrtt*[var = 1]<br/>@evar cov(e1,e2)=c(5)</p><p>'set starting values<br/>c(1) = 0.01<br/>c(2) = 0.06<br/>c(3) = 0.05<br/>c(4)=0.03<br/>c(5)=0.05</p><p>'estimate<br/>show sv.ml(showopts,m=500,c=1e-7)</p><p>&nbsp;</p><br/>

[此贴子已经被作者于2008-7-30 22:23:24编辑过]