Question

我创建了这个简单的资产增长模型并得到225分。唯一的问题是我想重复这10次，基本上存储前225个点并再次重复这个过程，这次只改变标准的正常随机变量，这里是做了什么

St <- 10
u <- 0.15
sigma <- 0.1
h <- seq(0,1,by=1/225)
t <- h[-1]

for (j in 1:10) {
  z <- rnorm(225)  
  for (i in 1:225){
     ST <- St*exp((u-0.5*(sigma^2))*t - sigma*z[j]*sqrt(t))
  }
}

Answer 1

您能否确认我的计算是否正确：

St <- 10
u <- 0.15
sigma <- 0.1
z <- rnorm(225) 
h <- seq(0,1,by=1/225)
t <- h[-1]

List = list()
for (i in 1:225){
  S <- St * exp( (u - 0.5 * (sigma^2)) * t[i] - sigma * z[i] * sqrt(t[i]) )
  List[[i]] = S
} 

Matrix = do.call(rbind, List)
Matrix

如果是，那么你想做同样的10次吗？由于你没有回应，我将继续并假设这是你想要的。这是代码：

St <- 10
u <- 0.15
sigma <- 0.1
h <- seq(0,1,by=1/225)
t <- h[-1]

BigList = list()
for (g in 1:10){

  z <- rnorm(225) 
  List = list()
  for (i in 1:225){
    S <- St * exp( (u - 0.5 * (sigma^2)) * t[i] - sigma * z[i] * sqrt(t[i]) )
    List[[i]] = S
  } 

  BigList[[g]] = do.call(rbind, List)
}
All = do.call(cbind, BigList)
All

如何创建用于存储数据的向量

1 个答案: