プロット複数の幾何ブラウン運動

Question

私は、次のコードを使用してRに幾何ブラウン運動のさまざまなシミュレーションを作成しました：

m <- 10 
n <- 1000 
mu <- 0.15 
sigma <- 0.3 
s0 <- 10 
T1 <- 5 
set.seed(70967993) 
MotiBr <- matrix(rnorm(m*n, 0, sqrt(T1/n)), n, m) 
MotiBr <- apply(MotiBr, 2, diffinv) 
temp <- (0:n)/n*T1 
MotiBr <- (mu - sigma^2/2)*temp + sigma*MotiBr 
MotiBr <- t(s0*exp(MotiBr)) 

は今、私はこのようなグラフに結果をプロットしたい：

これはGBMではないことに注意してください。私は得たい結果を公開するために投稿しました。誰かが私のシミュレーションを使ってそのようなプロットを構築するコードをMotiBrに投稿できますか？

PS：私はT1に値5を与えますが、私が意味するのは[0,5]です。ちょうどT1 <- 5でコード化するのは正しいですか？

Answer 1

私はあなたがmatplot機能を探していると思われます。

matplot(temp, t(MotiBr), type = "l", col = "black") 

は、ここで以下のプロット

Answer 2

を生成し、あなたのデータのプロットを作成する方法です。

まず、マトリックスを長いフォーマットのデータフレームに変換します。

データをプロットする1つの方法は、ggplot2パッケージです。

library(ggplot2) 
ggplot(dat, aes(x = time, y = val, group = idx)) + 
    geom_line()