标绘多个随机游动

Question

我相信新的R和观看YouTube视频用下面的代码：

require(yuima) 
m1 <- setModel(drift='mu*s', diffusion='sigma*s', 
    state.var='s', time.var='t', solve.var='s', xinit=100) 
x <- simulate(m1, true.param = list(mu=.1, sigma=.2)) 
plot(x) 

以上是亩乘以S，和sigma乘以秒。 我想将这些模拟中的1000个绘制成一条多色线条。任何投入将不胜感激。

Answer 1

有可能是一个更好的方法来做到这一点，但我最终挖掘到模拟对象，并拉出我需要的部分。 （一般来说，这是很危险的/不好的做法;如果包的作者给你访问的是得到你想要的比特的方法，你应该使用它们。）

设置：

require(yuima) 
m1 <- setModel(drift='mu*s', diffusion='sigma*s', 
    state.var='s', time.var='t', solve.var='s', xinit=100) 

运行一个函数仿真和返回数据（simulate确实有一个​​说法，但它并没有在这种情况下做什么...）

simhack <- function() { 
    s <- simulate(m1, true.param = list(mu=.1, sigma=.2)) 
    return(s@data@zoo.data[[1]]) 
} 

运行1000次 - 自动组合结果为1000列的矩阵：

r <- replicate(1000,simhack()) 

简介：

par(las=1,bty="l") ## cosmetic 
matplot(r,type="l",lty=1)