Вы можете использовать replicate
, который вернет вам 10 столбцов, для которых вы можете использовать colMeans
library(SMFI5)
simulateloop <- replicate(10, sim.vasicek(0.001987455,0.001614103,0.000186756,-0.001,30,1))
Чтобы использовать for
l oop, вы можете сохранить результат в список
simulateloop <- vector('list', length = 10)
for(i in 1:10){
simulateloop[[i]]<- sim.vasicek(0.001987455,0.001614103,0.000186756,-0.001,30,1)
}
, а затем используйте sapply
, чтобы взять mean
sapply(simulateloop, mean)