R语言写蒙特卡洛模拟

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

用R写股票定价公式之后（重复100次，每一次默认一年253天）想用plot把每一年的路径画出来，我自己写的函数如下
for (i in 1:100){
plot(1:253,data[i,],type="l")
}

但是只能显示出来一条线，但是我已经写了一个for循环了
找朋友帮忙写了一个

for (i in 1:100){
  if(i == 1) plot(x,data[i,],"l",ylim=c(0,120))
#先画出第一条模拟路径
  if (i != 1) lines(x,data[i,],ylim=c(0,120))


就可以产生100条线 ，请问一下R问什么要这么写？我看matlab直接用一个for循环，写个plot就能完成

本人小白，还请各位大神指点

扫码加我 拉你入群

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

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

分享0 收藏5 回帖

关键词：蒙特卡洛模拟 蒙特卡洛 蒙特卡 R语言 MATLAB

plot是高级绘图，添加线的话用lines。R绘图最好的是ggplot2

1. require(ggplot2)
   
2. require(reshape2)
   
3. 
   
4. T <- 1
   
5. n <- 1000
   
6. delta <- T/n
   
7. n_sample <- 5
   
8. 
   
9. set.seed(123)
   
10. rand_var <- matrix(rnorm(n*n_sample), ncol = n_sample, nrow = n)
    
11. browian <- apply(rand_var, 2, cumsum)
    
12. 
    
13. stock <- data.frame(id=as.factor(rep(1:n_sample, each=n)),
    
14.                     time=rep(1:n, n_sample),
    
15.                     ret=as.vector(browian))
    
16. 
    
17. p <- ggplot(stock) + geom_line(aes(y=ret, x= time,  colour=id))
    
18. 
    
19. ggsave("browian1.pdf", p)

复制代码

deem 发表于 2017-8-10 00:13
plot是高级绘图，添加线的话用lines。R绘图最好的是ggplot2

谢谢了      我试试  麻烦了

deem 发表于 2017-8-10 02:10

还有想问一下 plot难到只能话一条线吗？ 这个好不智能啊  那那个循环就白写了

世界和平X 发表于 2017-8-9 23:27
用R写股票定价公式之后（重复100次，每一次默认一年253天）想用plot把每一年的路径画出来，我自己写的函数如 ...

用matplot画多条线

世界和平X 发表于 2017-8-9 23:27
用R写股票定价公式之后（重复100次，每一次默认一年253天）想用plot把每一年的路径画出来，我自己写的函数如 ...

谢谢 我去试试

wangfs111222 发表于 2017-8-10 13:21
用matplot画多条线

谢谢 我去试试

deem 发表于 2017-8-10 02:10

大神 你这个代码写的太棒了  谢谢了  我好好学习学习