Question

我正在使用R中的Sim.DiffProc包来模拟Stratonovich随机积分。使用以下代码，我可以模拟从t = 0到t = 5的随机积分的5条路径：

fun=expression(w) strat=st.int(fun, type="str", M=5, lower=0, upper=5)

如果st.int()函数没有给出各种t中的值作为输出，我怎样才能得到t = 5中随机积分的值？

Answer 1

我不确定你的意思是t = 5。 $ X矩阵是时间序列的一部分：

> str(strat)
List of 8
 $ X           : mts [1:1001, 1:5] 0.0187 0.0177 0.0506 0.0357 0.0357 ...
  ..- attr(*, "dimnames")=List of 2
  .. ..$ : NULL
  .. ..$ : chr [1:5] "X1" "X2" "X3" "X4" ...
  ..- attr(*, "tsp")= num [1:3] 0 5 200
  ..- attr(*, "class")= chr [1:3] "mts" "ts" "matrix"
 $ fun         : symbol w
 $ type        : chr "str"
 $ subdivisions: int 1000
 $ M           : num 5
 $ Dt          : num 0.005
 $ t0          : num 0
 $ T           : num 5
 - attr(*, "class")= chr "st.int"

如果它是值矩阵的第五行是你的意思，它将是：

> (strat$X[5 , ])
          X1           X2           X3           X4           X5 
0.0031517578 0.0161278426 0.0003616453 0.0097594992 0.0012617410

模拟随机积分

1 个答案: