全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
24531 43
2011-12-15
请教:R程序下怎样用蒙特卡洛模拟股市大盘指数,望高人指教
二维码

扫码加我 拉你入群

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

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

全部回复
2011-12-15 18:51:47
顶起来
二维码

扫码加我 拉你入群

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

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

2011-12-15 20:59:08
底下程序取自matlab,
楼主有R基础,应该很容易转换.
Numerical methods in finance And Economics :
               A Matlab-based introduction
AssetPaths.m

% usage: SPaths = AssetPath(S0,mu,sigma,T,NSteps,NReps)
% description: returns NReps rows whose entries are a random path
% for asset with drift mu, volatility sigma and initial price S0,
% starting at t=0 ending with t=T in steps of dt = T/NSteps.


S0=50;      %initial price S0
mu=0.04;    %drift mu
sigma=0.3;  % volatility sigma
T=1;        %1 year
NSteps=252;
NReps=1000;

% usage: SPaths = AssetPath(S0,mu,sigma,T,NSteps,NReps)
% description: returns NReps rows whose entries are a random path
% for asset with drift mu, volatility sigma and initial price S0,
% starting at t=0 ending with t=T in steps of dt = T/NSteps.
SPaths = zeros(NReps,NSteps+1);
SPaths(:,1) = S0;
dt = T/NSteps;
nudt = (mu-0.5*sigma^2)*dt;
sidt = sigma*sqrt(dt);
for i = 1:NReps
  for j = 1:NSteps
    SPaths(i,j+1) = SPaths(i,j)*exp(nudt + sidt*randn);
  end
end
二维码

扫码加我 拉你入群

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

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

2011-12-15 21:11:24
估计要看一些统计计算方面的书,推荐《Introduction to Probability Simulation and Gibbs Sampling with R》,你可以搜搜相关的毕业论文,《统计模拟及其R实现》也可以作参考。
二维码

扫码加我 拉你入群

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

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

2011-12-16 09:07:36
epoh老师,您好!
Error: Missing MATLAB operator.这个错误是什么意思?如何更正?
>> % usage: SPaths = AssetPath(S0,mu,sigma,T,NSteps,NReps)
% description: returns NReps rows whose entries are a random path
% for asset with drift mu, volatility sigma and initial price S0,
% starting at t=0 ending with t=T in steps of dt = T/NSteps.
S0=50;      %initial price S0
mu=0.04;    %drift mu
sigma=0.3;  % volatility sigma
T=1;        %1 year
NSteps=252;
NReps=1000;
% usage: SPaths = AssetPath(S0,mu,sigma,T,NSteps,NReps)
% description: returns NReps rows whose entries are a random path
% for asset with drift mu, volatility sigma and initial price S0,
% starting at t=0 ending with t=T in steps of dt = T/NSteps.
SPaths = zeros(NReps,NSteps+1);
SPaths(:,1) = S0;
dt = T/NSteps;
nudt = (mu-0.5*sigma^2)*dt;
sidt = sigma*sqrt(dt);
for i = 1:NReps
  for j = 1:NSteps
    SPaths(i,j+1) = SPaths(i,j)*exp(nudt + sidt*randn);
  end
end

>>
以上程序可以行动,我想问一下:如何把模拟出的结果做出路径趋势图呢?这样看起来更直观。

我曾经看过您写的brown motion的R模拟程序,您并做出了图,非常直观。
二维码

扫码加我 拉你入群

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

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

2011-12-16 10:01:28
zhangtao 发表于 2011-12-16 09:07
epoh老师,您好!
Error: Missing MATLAB operator.这个错误是什么意思?如何更正?
>> % usage: SPaths  ...
zhangtao兄说的可是底下
股票价格模拟路径图.
   AssetPath.jpg

二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群