没头没脑的问题，好在哪里，没人回答狠正常，谁知道到底要什么帮助，再说这里都是统计和计量的，这种问题要么到金融版去问，要么去matlab版

顺便，陈毅恒的书在这里有review，http://www.riskbook.com/link/chan_and_wong_(2006).htm

这年头，搞统计的，搞计量的，搞经济的，搞会计的，都敢说是自己是搞金融的

[此贴子已经被作者于2007-2-4 8:40:34编辑过]

确实也发现,这里有个很不好的现象,下载的东西是很多下载的人也很多,不知道真正去看能懂的有几个...论坛上的问题都没人解决,下载那么多东西有什么用

回楼主,我只会用R/SPLUS算VaR,GARCH模型出来后预测1-step,multi-step volatility,分别算出他们的mean和variance,然后就简单了,具体有公式现在不记得了

如何用matlab做蒙特卡罗模拟计算VaR，能不能请高手指教一下~~~~~

山东师范大学 。郭繁的硕士论文

第一步：产生你选择分布的随机数
第二步：根据求出的方差和e=h^0.5*n求出扰动项，其中，e是扰动项，h为方差，n为随机数
第三步：利用均值方程和扰动项求出收益
第四步：将收益按照从小到大的顺序排列，求出给定置信水平下的损失，即为所求VaR

谢谢楼主分享！！！

这可是大牛glasserman写的， 而且附带variance reduction

我也在做这个  算套期保值  痛苦啊

哎  同志们   如果什么都这个都让别人告诉你 ，文章就写的差不多。有些时候还得靠自己

嗯，说得深入浅出。。。受益匪浅了

I CAN USE R to do it and guess matlab is pretty similar

看过好多关于用Eviews做蒙特卡洛模拟计算VaR的，但是具体程序，也不知道啊，希望高手解答。

请问楼上R来做需要具体程序，请指教

哇哈哈~ 大侠出现了，首先monte carlo模拟算VaR的方法很多， 看你的标题貌似你要先拟合garch再用monte carlo， 那不就是garch simulation了么，直接模拟了取quantile呗。
我擦，才注意到时08年的老帖了......

谢谢以上各楼，受益匪浅

1）设定一个rolling window
2）对这个window做GARCH拟合（用GARCH前，我会先检验一下数据）
3）monte Carlos 模拟
4）取quantile

tong wen

为什么不直接以excel模拟?

忘记了，

MATLAB做还真不太会 不过蒙特卡洛可以用水晶球来做

第一步：历史模拟发

function [VaR CVaR] = fun1( x , alpha ,M)
n=length(x);
%初始化Var，CVaR，rep使用了遗传算法工具箱中的函数
% VaR=rep([0],[1 n-M]);
% CVaR=rep([0],[1 n-M]);
VaR=repmat(0,[1 n-M]);
CVaR=repmat(0,[1 n-M]);
for i=1:n-M%用历史模拟，所以需要有前M=365天的数据辅助来计算第366天的VaR
  R=x(i:i+M);%366天的数据
  VaR(i) = prctile(R,alpha*100);%366天的数据按分位数求得VaR
  r=[];
%r=R(R<=VaR)

%找出数据中小于分位数的数据,给r;
for j=1:length(R)
  if R(j)<=VaR(i)
      r=[r R(j)];
  end
end
%对r求均值，即为CVaR
  CVaR(i)=(mean(r));
end

end



第二步：蒙特卡洛模拟，其中用到了历史模拟
function [VaR CVaR] = fun3(x , alpha,M)
n=length(x);
% v=var(x);%错了，
v=std(x);
u=mean(x);
xnew=normrnd(u,v,[1 n]);
[VaR CVaR]=fun1(xnew,alpha,M);
end

主函数调用

%get the data from the excel file
%Data.xls保存的是没有"%"的涨幅,2005.06.01-2011.06.24；
%X为什么只取了2006到2011年的呢？那就是因为滑动窗口大小为366，预测每个第366的数据。
clc;
clear all;
Data=xlsread('Data.xls');
Data=Data/100;
X=xlsread('NUM.xls');
%
alpha=1-0.95;
year=365;
X=X(1:length(Data)-year);
[VaR1 CVaR1]=fun(Data,alpha,year,1);
% [VaR2 CVaR2]=fun(Data,alpha,year,2);
[VaR3 CVaR3]=fun(Data,alpha,year,3);
% [VaR4 CVaR4]=fun(Data,alpha,year,4);
h=figure(1);
set(h,'color','w');
plot(X,VaR1,'b-*');
hold on;
% plot(X,VaR2,'r-*');
plot(X,VaR3,'m-*');
% plot(X,VaR4,'y-*');
xlabel('06/06/2006--24/06/2011');
legend('VaR-HS','VaR-NORM','VaR-CN','VaR-CF',2);%2表示在左上角
hold off;

h=figure(2);
set(h,'color','w')
plot(X,CVaR1,'b-*');
hold on;
% plot(X,CVaR2,'r-*');
plot(X,CVaR3,'m-*');
% plot(X,CVaR4,'y-*');
% xlabel('06/06/2006--24/06/2011');
% xlabel('06/06/2006--24/06/2011');
legend('CVaR-HS','CVaR-NORM','CVaR-CN','CVaR-CF',2);
hold off;

GARCH
MATLAB 用GARCH（参数），求得条件方差方程，迭代算出来后，带入就好了