2017-06-19 22 views
1

我已经创建了这个简单的资产增长模型并获得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)) 
    } 
} 
+0

圣< - 10 Ú< - 0.15 西格玛< - 0.1 ħ< - SEQ(0,1,由= 1/225) 吨< - H [-1] (i = 1:225)ST <-St * exp((u-0.5 *(sigma^2))* t(对于(j,以1:10为单位) - sigma * z [j] * sqrt(t)) } } –

回答

0

,你能否证实,如果我的计算是正确的:

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 
+0

Hi @ AK88,非常感谢。它的辉煌,它的作品 –

+0

随时接受我的回答:)) – AK88

+0

嘿,希望我能。但它表示“记录的声望低于15的人的选票被记录下来,但不会改变显示的分数后的宣传。” –